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Etude des intégrales a’ une équation aux dérivées 
partielles du quatriéme ordre dans un domaine 
illimité et quelques propriétes de leurs dérivées 


J. WoLSKA-BOCHENEK 


Commumicated by R. FINN 


1. Introduction 
Dans le travail [7] nous avons étudié |’équation aux dérivées partielles 


yA Ay(A, t) es =o) (1) 


; ‘ oe? een ae 
y étant un paramétre, A= i an END désignant l’opérateur de Laplace, connue 


dans l’hydrodynamique d’un liquide visqueux. En s’appuyant sur la solution 
fondamentale bien connue 


O (Ao 8; i no dq (2) 


UAB 


2)» (t—7) 
nous avons étudié les intégrales de l’équation (1) analogues aux potentiels de 
charge plane et de simple couche relativement a l’équation de chaleur dans un 
domaine borné. 

Dans le travail [7] nous avons profité des résultats récents de W. PoGORZELSKI 
[Z], [2], [3], [4], et du travail de J. SCHAUDER [6]. 

Dans ce travail nous nous proposons d’étudier les propriétés du potentiel 
de charge plane relativement a l’équation (1) dans un domaine non-borné, les 
propriétés du potentiel de double couche relativement a l’équation (1), de méme 
que les propriétés des dérivées de ces potentiels, et quelques nouvelles propriétés 
du potentiel de simple couche. 

Comme résultat de nos études, nous allons obtenir quelques nouvelles propriétés 
du potentiel de simple et de double couche relativement a ]’équation de chaleur, 
dont les plus importantes seront les propriétés de la dérivée normale et tangentielle 
du potentiel de double couche. 

Dans nos considérations nous allons supposer toujours que la courbe fermée 
C limitant le domaine non-borné D posséde une tangente continue en tout point 
PEC, et que l’angle que font les tangentes aux deux points P et F, de la courbe 
C satisfasse a la condition de Hélder suivante: 


|dpp,| SconstrSp, O<aS1. (3) 
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2. Etude des intégrales de 1’équation (1) 


Potentiel de charge plane 
On appelle potentiel de charge plane relativement a l’équation (1) l’intégrale 
suivante 


VA, = chic Ader el Baca nade (4) 


ou D est un domaine non-borné, limité par la courbe fermée C vérifiant la con- 
dition (3). 

La fonction @(B, t), dite densité de charge, est determinée dans le domaine 
[BED, 0XtST]. Si la densité @ est une fonction bornée dans le domaine 
[BED, 0<1tSXT] et intégrable dans toute partie bornée et mesurable de cette 
région, alors le potentiel (4) est determiné dans tout le plan pour OS?<T, 
et posséde la limitation: 

V(A,)| < Ky +H sup|o! | Oe Ax 2o\pour ize KR 6) 
AS20968 1747; 
et admet en tout point les dérivées premiéres par rapport aux coordonnées du 
point A, continues, données par la formule: 


V;(A,t) = =f fou A,t; B,t) 0(B, t) dogdt. (6) 


En outre nous concluons immédiatemment que les dérivées (6) admettent la 
limitation 
0 AaA) pour ify n ey 


IV (4,0)| < Ky A+ sup || a 


A> pour 2c he 


Il en résulte de méme que le potentiel (4) et les dérivées (6) satisfont 4 la con- 
dition de Lipschitz-Hélder 


|V(A, #) — V(A,, | < Hy é*** sup|o| 744, 
V4| (A, #) — Va (Ay, 8) |< H,8"** supl|o|744, 0<d0<1, (9 


— 
(ee) 
Sey (LS 


Ky, Ho, K,, Hy, Ry étant des constantes positives. 

Selon la méme supposition concernant la densité 9, le potentiel (4) admet 
les dérivées secondes par rapport aux coordonnées du point A et la dérivée 
premiere par rapport au temps: 


t 
V4 (A, 2) at JJ ora (A, t; B,t) 0(B, t) dopdt, 

(10) 
als if w, (A, t; B, t) 0(B, t) dogdt 


0 


en tout point extérieur au domaine D et pour O<f<T. 


Ces dérivées vérifient l’équation donnée (1), puisque la fonction AV vérifie 


’équation de chaleur: 
yA(Av) — SAY) — 


en tout point A extérieur au domaine D pour 0<#<T. 
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Théoréme 1. Sz la densité 0(B, t), bornée et continue dans la région [BE D, 
~~ OSTST), vérifie la condition de Holder: 


|o(B, t) — e (By, t)| < const rep, Ose t; (11) 


dans une région [BE D’, 0<tST], ou D’ est la partie du domaine D, alors le 
potentiel (4) admet les dérivées secondes par rapport au point A et la dérivée premiére 
par rapport au temps en tout point A intérieur du domaine D’' pour O<tST. 


Démonstration. Nous nous appuierons sur le théoréme 2 de notre travail [7]. 
Considérons donc un domaine borné et mesurable D” contenant a l’intérieur le 
point arbitraire A€D et situé dans D’. Nous écrivons le potentiel (4) sous la 
forme d’une somme des intégrales 


VAL eV eatA, te Vom aa, 2) (12) 


étendues au domaine D” et au domaine D—D"”. Or, d’aprés le travail [7] les 
termes de la somme (12) admettent les secondes dérivées par rapport au point 
A et de méme la premiére dérivée par rapport a ¢. Il en résulte de méme que 
e laplacien du potentiel (7) s’exprime par la formule 


AV (A, t) =| ff Aol A,t; B,t) 0(B, t) dogdt (13) 


dans le domaine A €D’, 0<iST. 


Corollaire. La fonction (12) est identique au potentiel de charge plane rela- 
tivement a l’équation de chaleur 


ae 


AV(A, 1) = $f fee "t—*) 9(B, 1) dag dt 


t—T 


et a lintérieur du domaine [A €D’, 0<i<T], a supposé que (11) satisfasse, 
d’aprés les études du travail [7] a l’équation de Poisson 
(AV) 


A(AV) a b= =Aar 0. (14) 


Aussi, d’aprés la remarque précédente, en s’appuyant sur les résultats du travail 
[3], nous pouvons constater que le laplacien du potentiel de charge plane posséde 
les dérivées premiéres par rapport aux coordonnées du point A, de la forme 
suivante: 


(AM =f If ol 0B, 2) doy dr (15) 


dans tout plan, et pour OSt<T. 


D’aprés les études du travail [3] nous tirons les limitations suivantes des 


intégrales uniformément et absolument convergentes: 
5 | <Ks f~’ sup|o Orda St rg = Tig (16) 
|AV|<K,¢~“*sup|o Ree) Sl hae ine (17) 
(Ge), | < Ket sul e ee A 1) SE at ipa ys (18) 
|(AV)4|< K,#~*sup|o VAS] Meee hy: (19) 
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aye OV oV 
K, et Ks étant des constantes positives. En outre les fonctions Be er ay 


satisfont aux conditions: 
@V(A,t)  @V(Ay,4) 
at at 
| éV (A, t) | E V (A, 1a 
A 


(20) 


~*sup|o| 7, 


< H,f~*sup|o| 74.7 (21) 


i 2 


ot at 
ot 0<0<1, H, et H, étant des constantes positives. 


Remarque. D’aprés les résultats récents de W. PoGoRzELskI [9] on peut 
démontrer les propriétés précédentes du potentiel de charge plane, méme dans 
le cas ou la fonction o(B, t) n’est pas bornée. 


Potentiel de simple couche 


On appelle potentiel de simple couche relativement a l’équation (1) l’intégrale 


=f folds: Q, t) u(Q, t) dlo dt. (22) 


Les propriétés de ce potentiel ont été discutées dans notre travail [7]. 


suivante 


Nous allons étudier quelques propriétés des dérivées du potentiel (22). 
Evidémment, si nous supposons que la fonction w(Q, tT) est bornée et intégrable 
dans le domaine [Q0€C, 0XtXT] nous pouvons démontrer que le potentiel 
(22) Rees la dérivée par rapport a la variable ¢, donnée par la formule 


t r4Q 
ie was meer Sy 
A d= = f fos 0,1) u(Q,1)dlodt 2d J es dlgdt (23) 
0 


en tout oe, ALE IMAC, OSS I. 
La dérivée (23) a la forme du potentiel de simple couche pour l’équation de 


chaleur et d’aprés les études du travail [/] vérifie la condition de Hélder de la 


f , , e z if 
ae |W (A, #) — Wy (A, &)| < Ag t*74.4, sup| | + hg sup|u| |t—4,|"”, (24) 


0<0<1, 0< <1, et posséde la limitation 


|W |< ky? sup| ul, 


(25) 
k, et hy étant des constantes positives. 


Nous allons démontrer que sous les autres suppositions concernant la densité 
Ht, la dérivée (23) peut satisfaire a la condition de Hélder par rapport A la variable 
t avec l’exposant arbitrairement inférieur a l’unité. 


Théoréme 2. Sz la fonction u(Q, t) est continue par rapport aux coordonnées 
du point QEC pour 0<tSXT et satisfait a la condition de Hélder 
|H(Q, 7) —#(0,8|S (26) 


0<t<t, 0<f’<1, alors la dérivée (23) du Sore de simple couche (22) satisfatt 
a la condition de Hélder suivante: 


Wi (4,4) — Wi (A,4)] S MeL, + -* sup [ul] |t— 4)” (27) 
Oty >t, Sm Ag, O29 


} 
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Démonstration. Etudions la différence 


peli 
W; (A, t) —W; (A, 4) == - f {ae AMG (Q, T) dl, lo 
) “OC 
Te ole tee TAQ 
4v(t,—T) 
gD fife ee PI y(0.3) ah ae 
0c 


rAQ 


=-3 i fee PO 0.9 — 40.0) digae + 


a5 : ffs e "69TH (Q,7)—K(Q,4)]dlpdt+ (28) 
ee & 
4 : 1 Ae 
Av T) 
m3 Af rye [u(Q,4) —“(Q,#)]dlgdrt+ 
0. ¢ 
a ees A 4¥ (t;—T) = 
1 sl fae L(Q, t) dlg dt 
0c 


a, ae 
4v(t—T) 
= f{ [ee 1(Q,t) dlp dr. 
0 Cc 


Les trois premiéres intégrales dans l’expression (28) satisfont a la condition de 
Holder avec l’exposant 0<#<1, ce qui résulte immédiatement des résultats du 
travail de W. PoGoRZELSKI [/]. Considérons donc les deux derniéres intégrales 
(28). Evidemment nous aurons 


= {fae rae) 10(Q,8)dlgdt — > ffmae' u(Q,t)dlodt 


a ft (A,t; Q,0) —(4, hh; Q,0)]u(Q, A) dlp (29) 
C 
1 ae ee: 
| oe Ate (¢ —t,) u(Q, 8) dg, 
CG * 


ou ?¢, est choisi arbitrairement dans l’intervalle (¢,, 7). D’aprés la limitation 
connue de la fonction sous-intégrale nous obtenons 


essere 


(30) 


ou <A,x< 1, et arrivons a la thése de notre théoréme. Enfin, en supposant que 
la fonction (Q, t) satisfasse a la condition de Hélder 


|(Q, 7) —H(Q, 2) <7 100, (34) 
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on arrive aux inégalités suivantes: 
|W (A, 4)| < hi (sup|e| + H,) ¢, (32) 
["7s(A, t) — Wa (Ay, 8) <M (sup |u| + H,) O74, (33) 
ki, et hi étant des constantes positives dépendant de la courbe C, 0<A<1. 


De méme, d’aprés les résultats de W. PoGorzELSKI [3], sous la supposition 
(31) pour la fonction 7, nous arrivons aux inégalités 


Sp (Pt) — 3, WB] < hes oup| ye] O* + H,) Pe, (34) 
d / —j) Ay 

a7 | < Aa (sup |u| A + HP), (35) 

|(AW)4| < kg (sup|u|e-*+ A, e7*), (36) 


kg, ky, hg 6tant des constantes positives, 0<Ay<1, O<A,<1. 


Potentiel de double couche 
On appelle potentiel de double couche relativement a |’équation (1) l’intégrale 
suivante: 


YAQ 


=f fico algae f f HIE 7 2(0,2)dlgdt (37 
0C 


oul yy désigne |’angle que fait la normale intérieure au point Q de la courbe C 


avec le vecteur QA. 
En outre, nous avons démontré dans le travail [7] que le potentiel (37) 
satisfait aux inégalités 
|U(P,t) — UB, | < mA? sup |e] 7h, (38) 


[U(P, |< kee 


supposition faite que la fonction ¢(Q, tT) est continue et bornée dans la région 
[QEC, 0OStST], h, et k, étant des constantes positives. 


Dérivée normale du potentiel de double couche 
Ecrivons le sg (37) sous la forme suivante: 


U(A, #) ~ f fsme re rea nee 


YAQ 


~ ff eee t liga + ff IL 1) dlgdt 


YAQ 
COS ay fet 
— 040) 4v(t—r) = 
a Se <4 ,t)dlgdt- 
i 10 | \f@0 odt+ 
oc 
t 


spl thara (0,1) —C(P, D)]dilpdt +22 [aie 1) dt. 


0 
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I] en résulte que la dérivée du potentiel (37) dans la direction d’une normale 
au point P de la courbe C (le point A est situé sur cette normale) s’exprime par 
les intégrales 
dU(A,?) cos y; fs 
— anp = {fal eye Zz =< C(Q, t) dlgdt + 


YAQ 


2 


rAQ 
COS Yo COS Yp .~ Gy(t—z) 
ih 2 (t—T) s C(Q, t) dlodt + (41) 


Cos 0 
LS wet aa EWG, 2) EAP) hee 

En s’appuyant sur les études de J. ScHAUDER [6], concernant la dérivée normale 
du potentiel logarithmique de double couche, nous pouvons démontrer les théo- 
rémes suivants: 

Théoréme 3. Sz la fonction €(Q, t) posséde la dérivée par rapport a V’arc 
vérifiant la condition de Hélder, 
d d 
az $ (2: uke azo (G17) 


continue par rapport a la variable t, et si la courbe C satisfait a la condition (3), 
ow O<B Sa, alors la dérivée normale du potentiel (37) tend uniformément vers la 
valeur linute 


Oi xt 
oe ee ffd eg le Ss ~1]£(0,) dy dt + 
dnp 


as i, ro OQ. (42) 


FLED AAS adup =a 
cos y°, cos 9 rPQ 
if Byes Se ENON Valo dee (43) 
ads) . _ au(P,n 
7 ffl ¥PQ )E@. 2) ag Wd T) | dlo dt a, ; 


et cette valeur limite existe au sens de la valeur principale de Cauchy, y% dant 
—> 
angle entre la normale intérieure au point Q avec QP, et y% étant l’angle entre la 
— 
normale intérieure au point P avec PQ. 


Démonstration. Nous allons montrer d’abord que la valeur limite oe 
Pe 


existe au sens de la valeur principale de Cauchy. Pour la troisiéme intégrale 
de l’expression (43) cette propriété résulte directement du théoréme de J. ScHau- 
DER [Joc. cit.] concernant la dérivée normale du potentiel logarithmique de double 
couche. Quant a la seconde intégrale (43) — elle est absolument convergente 
a la seule supposition que la fonction €(Q, t) soit continue et bornée et que la 
courbe C satisfasse a la condition (3). Evidemment, d’aprés la limitation de la 


fonction sous-intégrale 
0 pats ue 
COS YQ COS vey 49 (t—7) = const PSA 3 (44) 
2v(t—t) (t—2)4 rE gr? 
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ow Ai a Ls 
nous pouvons constater, en choisissant A<1, a+A>%, que Vintégrale 


t i 
* COS 7/9 COS Y'p Ay ( P, = dl ve 
13 2v(t—o) é "o( ) dla 
0€ 
est absolument convergente. 
La premiére intégrale (43) posséde la méme propriété. En effet, écrivons 


t a eS 7 St oe 
[fe my _¢ W-91C(0,2)dlydt 
J J dup \ rpo 
0 3 
“ - 2c0sy% cosy ae 
aa cosh cosy Trte  -r0=9 (9,4) digdu + sm 
. YP Q 2% (t—T) 


2 
Vv p* 


9 a al’ 
+ {fH YP*O 2. see=*) £(Q; T) digdt =J,+Je, 


YP QO 29( (¢—7) 


ot. P* est le point arbitraire de l’arc PQ. La distance 7psg étant comparable 
avec pg — nous pouvons constater que la fonction sous- -intégrale dans ]’intégrale 
J, posséde la limitation (44) d’ot il résulte que la fonction jf, est absolument 
convergente. La fonction sous-intégrale dans l’intégrale J, posséde la limitation 

rp+Q 
cost = 7pxy 2 a0) 
¥pqg 2v(t—T) 


a (46) 
PQ 


Nous en pouvons déduire, en choisissant 3<A<1, que l’intégrale J, est absolu- 
ment convergente. 

Maintenant, de la méthode bien connue dans la théorie classique du potentiel, 
on peut démontrer que 


eee dU(P, t) 


ae a | <e Sl'A P| << ey 


c’est-a-dire la thése du théoréme 3. 
En s’appuyant sur l’expression (43), nous pouvons démontrer la limitation 
suivante: 
dU(P = 
| U(P, t) | = (?- *sup|C| a ¢sup | at 7 | +e] (47) 


dnp 
k, étant une constante positive. 
Théoréme 4. Sz la fonction €(Q, t) satisfait a la condition (42), la courbe C 


satisfait a la condition (3), alors la valeur limite (43) de la dérivée normale du 
potentiel (46) satisfait a la condition de Holder suivante: 


aU(P,t) .. dU(B,?) 
{ dnp dnp | = hy th * sup \c| So sup | =|+ tH) 4, VPP, (48) 


ou h, est une constante positive. 


Demonstration, La valeur limite (43) se compose de trois intégrales, dont 
es deux premiéres possédent la singularité faible. Donc, d’aprés les méthodes 
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bien connues dans la théorie du potentiel, il est facile 4 démontrer, que’elles 
satisfont a la condition de Hélder avec l’exposant #> 0, arbitrairement inférieur 


a lunité. Quant a la troisiéme intégrale (43), la démonstration du théoréme 
énoncé sera la méme que dans le travail [7]. 


Dérivée tangentielle du potentiel de double couche 
Théoréme 5. Sz la fonction £(Q, t), vérifiant la condition de Hélder 


I5(Q, 7) —0(Q1, 0)| < Ay rp, (49) 


est continue par rapport a la variable tv et si o+B,>1, « étant l’exposant de la 
supposition 3, alors le potentiel (37) en tout point P de la courbe C elle-méme, 
possede la dérivée par rapport al arc, vérifiant la condition de Holder, 


NONIZ. CROCE ai ; =i a+ By — 
v S: <h,[f*sup|¢| tH] 7ep" 4, (50) 


ou la constante positive Ing dépend de la courbe C. 


Démonstration. D’aprés l’expression (40), nous pouvons écrire la dérivée par 
rapport a l’arc du potentiel (37), en tout point P de la courbe C sous la forme 
suivante: 


oe Cos YO pect 
=f fae pate 1)£(Q, 2) dlgde + 
COs } sin 7 sthe 
0 P , 4v(t—z) | 
if 2v(t—t)  ~ C(Q, t) dlg dt + (51) 


+ Lf dl 22) 800.9 2 nage. 


\ YPO> 
oc 


Donc, d’aprés les considérations analogues, que dans le théoréme 3, nous pouvons 
constater que les deux premieres intégrales de l’expression (51) aux singularités 
faibles, satisfont a la condition de Hélder avec l’exposant arbitrairement inférieur 
al’unité. Donc, la troisiéme intégrale, d’apreés les études de J. SCHAUDER [Joc. cit.], 
satisfait 4 la condition de Holder avec l’exposant «+, —1, c.q.f.d. 


De méme nous obtenons facilement les limitations: 


|U4 (A oe es ie (53) 


Dérivée par rapport a la variable t du potentiel de double couche 


D’aprés le résultat du travail [7] la dérivée par rapport a la variable ¢ du 
potentiel de double couche s’exprime par l’intégrale: 


2 
v2, 


t AQ 
oo a r40 COS YQ Sonar 
Ui (A.N= ff afar 6(Q, 1) dlgdt (54) 
0c 
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et vérifie les inégalités: 


|U/ (P, t) — U/ (P,,4)| <hg sup || |r? ee P24 tafe] (55) 


IU) (P,d)| <Aysup|é|¢ 2 * (56) 


ol 0<9<1,0< 0 <1, hg, kg étant des constantes positives. Par la méme méthode 
que dans le théoréme 2, nous pouvons démontrer le théoréme suivant: 

Théoréme 6. Si la fonction £(Q, 1), définie dans le domaine [QEC, OS tS T], 
est continue par rapport aux coordonnées du point Q et vérifie la condition 


Ducat, (57) 


alors la dérivée (54) satisfait a la condition de Holder suivante: 


”, (58) 


|U; (A, t) — U/ (A, 4)| < hg (Hz + t-*** sup|C|) |t—4, 


ou t>t, 1—ga<Ayy<1, O<PS1. 


Théoréme 7. Sz la fonction €(Q, t) vérifie la condition (49) et la condition 
(57), st la courbe C satisfait a le condition (3), si a+f,>1 et si a+28’>1, alors 
la dérivée du potentiel (37) par rapport a la variable t, ayant la forme (54) posséde 
la dérivée par rapport al arc en tout point P de la courbe C, vérifiant la condition 
de Holder 


a7 ( (Poti = 7U (| <h, (#- *H,+4¢-*sup|C| + H,) rep. (59) 


Démonstration. Ecrivons donc la fonction (54) sous la forme suivante: 


I= ff COsiyg, 7400, 40a [f(Q, t) —€(Q,#)] dlgdt+ 


YAQ 4y( ie 
(60) 
YA0 


[Te eM etena fe ze.905 


I] en résulte que la dérivée par rapport a l’are de la fonction (54) en tout point 
P de la courbe C elle-méme s’exprime par les intégrales aux singularités faibles 
suivantes: 


ai LU.) -f fa 
(61) 


“fa[gBe% ajfesaee f4(eAjet0n-cooia 


YP OQ 


Ypg COS 7h 


ee ZG ital [f(Q, t) —C(Q, é)] dlgdt + 


En effet pour la premiére intégrale cette propriété est évidente d’aprés la limi- 
tation de la fonction sous intégrale 


d | Ypo.cos BL i } 
fae [eager 7 OO ah eta ae © 
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: d COS Yo SIN yp—sin } Be 
remarque faite que ai (cos y) = vo = es et a+2f’>1. En choisissant 
3—a—2p’ we solar 
— “5 B <A<1 nous pouvons constater que la premiére intégrale (61) est uni- 
formément convergente. La seconde intégrale (61) se compose de deux membres 

PQ ; PQ 
d [cosy | — at __ cosy sin y ~ 4"? 
fa| rpg \° See ae if Laine pel Cia la 
C C 
2cos y% sin yP—sind 7 2820 9) 
| [= YO YP PO 2 4vt e(Q, t) dq, 
e YPO 2vt 
6 
_ bo 
oti la différence f mee =) est transformée de la méme facon que dans |’ex- 


pression (45). 
En tenant compte que 7p+¢ est comparable avec rpg, nous arrivons aux limi- 
tations suivantes: 
1PQ 
0 (Oe es 
COS sin Wiphes 
eee iy fe 
2vt 


const 
rea cee (64) 
| 2cos y sin yp—sin 0 rp«o 2 ave 


| YPQ Qvt 
4— 


En choisissant 4}, de l’intervalle ( ZS 1), nous pouvons constater que l’intégrale 


(63) est uniformément convergente. La derniére intégrale (61) d’aprés les études 
du travail [7] posséde la méme propriété. 

On peut démontrer, par une méthode analogue a celle employée dans la 
théorie du potentiel, que les deux premiéres intégrales satisfont a la condition 
de Holder avec l’exposant 0<?< 1 et la derniére, d’aprés les études de J. SCHAU- 
DER [loc. cit.] 4 la condition de Hélder avec l’exposant «+ 6,—1. 

-Il en résulte donc la thése du théoréme 7. Evidemment on arrive a l’inégalité 
suivante: eethee Ros tee as Ahm Ae ig 
5, (UP, A)]| < Re? * He + -* sup || + H,). (65) 

Théoréme 8. Si la fonction €(Q, t) vérifie les conditions (42), (57), supposition 
faite que x<’ <\, la courbe C vérifie la condition (3), supposition faite que 0<P Sz, 
alors la dérivée du potentiel (37) par rapport a la variable t, ayant la forme (54), 
posséde la dérivée normale an [U, (A, t)], qui tend uniformément vers la valeur 
limite 


d iy witas j 
iD cans AO) = aap 8 
! Nenre Seis 
s ee ee Or 49(t—T) [o(Q, t) —€(Q, 8) ] alg dt + 
Jed) dnp. | 4(t—)2 i 
i ‘ (66) 
pel ate COS VO ema 
| | = (« 4 1)]£(0.9 alot 
_ ff 4_[cos yd Lee (P.t) dl 
+f i (%2 |e(@4 C(P, t)] dl, 


et cette valeur limite existe au sens de la valeur principale de Cauchy. 
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Démonstration. Nous allons démontrer d’abord que la valeur limite (66) 
existe au sens de la valeur principale de Cauchy. Cette propriété est évidente 
pour la derniére intégrale (voir [7]). Les autres intégrales sont absolument con- 
vergentes. Notamment, d’aprés la méthode appliquée dans les theoremes 3 et 6, 
nous pouvons constater que la fonction sous intégrales dans la premiere inte- 
grale (66) posséde la limitation 


2 
’PQ 


d | 1PQ SVQ, av(t—#) | 9 220m SOUEE ia ioe , 
: : a Beare Fale A <i =e T SE yy, 
dnp 4v(t t)? [f(Q ) ( ) ( t)? | rpg” 2p’—20 roo? 2B 


En choisissant 2 —f’<A<1, nous pouvons constater que la premiére intégrale 
(66) est uniformément convergente. 


La fonction sous-intégrale dans la seconde intégrale posséde la limitation 


SA ae gM & 
_d_ |cosyQ (- arty 
anup YPOQ ‘ 


En choisissant +<A,<1, nous pouvons constater que la seconde intégrale (66) 
est aussi uniformément convergente. Maintenant, d’aprés les méthodes bien 
connues dans la théorie du potentiel, nous pouvons démontrer, que 


d d 
A—>PEC dnp dnp 


const 


ee 


1 ; 1 


9 3 ps 
YPO YPO 


[U4] = 3 [U(P, 2] 


c.q.f.d. 


Théoréme 9. Sz les conditions supposées dans le théoreme 8 sont satisfattes, 
alors la dérivée normale de la fonction (54) satisfatt a la condition de Holder suivante: 
LATED 4 See iee e Je Hi+ j-sup|t|+H;] rn. (67) 

i i aa ila Gp Se ; Peas : ; 
Démonstration. La preuve du théoréme résulte immédiatemment des con- 
sidérations précédentes, concernant les intégrales aux singularités faibles, et des 
études de J. SCHAUDER [Joc. cit.]. Nous aurons de méme les limitations évidentes 

d y —h , 1 cy ope 
Laie CNP A)]] < Ry [PAH + sup |e] +H, (68) 
, , , aa | Lae, 1 ore *. sf 

|[U'(A, Ta | < Rs G “Hy + 4. sup|¢| +H; + H,|. (69) 


Théoréme 10. Selon les mémes suppositions que dans le théoréme 7 la dérivée 
par rapport a Varc de la fonction U/(A, t) satisfait 4 la condition de Holder par 
rapport ala variable t 


d 9 


1, , A aw) me , 702 = oO 
ay UP.) — 5 (OP, 4))| < hs [te HE + sup |e] + He]|4¢—4l4%, (70) 


ow e€=min[a+2p’—1, A(a+9), 2(a+f6,—1)], ee Ei ae SE hi est une con- 
stante positive. ¥ 


Démonstration. Nous allons démontrer le théoréme cité suivant la méthode 
indiquée par W. PoGorzELski dans le travail [3]. D’aprés l’expression (61) nous 


pouvons décomposer la dérivée a [U;(A, t)] en trois parties que nous désignons 
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par f,, J2., &. Etudions la premiére différence 


mA? pre) 
- [fa 


[fe eee owns ee [5(Q, t) —C(0,4)] alo dt 


EO COs vO 


4y(t—a)? ~ me gee t) —€(Q,t)] dlgdt + 


- tke (74) 
¥pQ COS y? i) 
-[fi pore wi "|ee )—C(Q,0] digdr + 
{fale Din oo t) —C(Q,2)|dlgdt + 
ffé e929 WEN] (0,2) —£(0,A)]dlgde— Jit het 
dl 4y(t,—)? ail 94@tT=hfitheths: 
D’aprés la limitation (62) nous aurons 
9 LPP Sak a+2f’—-1 ele Sil 
Wat oe te constHst ok ae SG) 


Pour étudier la seconde vote (67) considérons le cercle JJ de centre P et 
dont le rayon rj n’est pas fixé pour le moment. Décomposons les intégrales 
Jie et Jz3 en deux parties: 


Stat he= he + hE + et hs) (73) 


prises suivant la portion C; de la courbe C intérieure au cercle /7 et la portion 
C —C,, extérieure au cercle J/. 

Il suffit d’étudier le cas oti vz est le rayon suffisamment petit 77< 6, ot 6 
est le rayon du cercle K assez petit pour que la partie de la courbe C intérieure 
au cercle / puisse étre coupée par chaque paralléle 4 la normale en P en un 
point au plus. Evidemment nous aurons 


Woes const Het yet th tas 
he ee CONSE Te yen err 


See 28 


(74) 


ou 


Transformons maintenant la somme (J, + J;3)°-“@ de la facon suivante: 


Srot+ Sis) OH iy fal 


0 C—Cy 


ty rea a rhe 
a / i d | rq cos 7h 2 a= d | YpQ COS YQ e Ph-dlly (75) 
Li a Fa Ped her yp sera il eas 

C-—Cy 


x (6(Q, 7) —6(Q, A)] diode = fie + fig. 


PQ | 
gee 7) (0,4) —£(0,0) dlgadv + 
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La premiére intégrale (75) satisfera a linégalité: 


\Jis|<conste “i, |7— Ale —— : (76) 
VII 


. . , ” , / s oN 
En appliquant a la fonction sous-intégrale dans lintégrale Jiz le théoréme des 
accroissements nous arrivons a l’inégalité suivante: 


ty bes ’ 
| al <const |¢— 4] He f poe al = dlgdt 


i C=en 


(77) 


< const f~? H,|t —4, |" 


ot la constante ne dépend pas de 777. 


Choisissons maintenant le rayon 777 en posant 


71 = const |t — t/? <6. (78) 
Evidemment nous avons 


16 S21 >) 5146, 76224 —a9= 6 


Il en résulte d’aprés (72), (74), (77), (76) 


a+2p’—1 . &+28’—-1 


|A(P, 4) — f,(P,4)|< const’ ¢ He |e—4 (79) 


Appliquant la méthode précédente aux intégrales J, (P, t) et J,(P, t) nous arrivons 
aux inégalités 


| J2(P, t) — h(P. 4)| < const [¢-* sup |¢| + H, #7] |t— hapieeee (80) 


a+p,—1 a+f,—1 
Ja.) — (P,4)| <const [Pew] |t—4P 4. (81) 


En rapprochant les résultats (79), (80), (81) nous arrivons a la thése du théoréme 10. 


La démonstration analogue appliquée a la dérivée normale du potentiel 
U/(P, t) nous améne au théoréme suivant. 


Théoréme 11. Selon les suppositions admises dans le théoréme 8 et la sup- 
position 
dt(P,t)  aC(P,t) 


GH dl 


|< Hele — a, 0<BS1; (82) 


la dérivée normale de la fonction U/(P, t) satisfait a la condition de Hélder suivante: 


d 


F ad 5 s! ey) = fond ioe bi & 
gq UP, (UP, h))| <h[f * Ae+-e-* sup || + H+ Hi] |¢—4 [© *, (83) 


ow €;=min (2p —1, 2B, 4B) e £<A,y<1, he étant une constante positive. 
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1. Introduction 


Dans la théorie du mouvement plan d’un liquide visqueux incompressible, 
on déduit |’équation: 


oAy op ocAy op OAy , 
= : ) 1 
Des Oy ox ap ap (1) 


ou p(x, y, ¢), dite fonction du courant, est définie par les égalités: 


7) é 
ji ae ja : (1’) 
: oy : Ox 
> = = = 
ou v, et v, sont les composantes du vecteur de vitesse v, p(x, y, t)=rot, F, ot F 
est le vecteur des forces extérieures; v est le coéfficient de viscosité cinématique. 


Dans ce travail nous nous proposons de trouver une fonction w(x, y, f), qui 
a l’intérieur du domaine non-borné D, limité par l'ensemble fini des courbes 


fermées disjointes C,, Cy,..., C, satisfasse a l’équation (1) et en tout point P de la 
frontiére C=C, +C,+---+C, satisfasse aux conditions limites suivantes: 
fim «ee im, ee oe 


ASPecow at ayPpec dn 


ou la premiére des conditions (2) est équivalente a la condition: 


lim w(A,t)=g(d); 


A->PE€EC 
les fonctions g(t) et G(P, t) étant données. 


En outre nous demandons que la fonction w(x, y, ¢) satisfasse a la condition 
initiale 


lim y(A,#)=f(A); AED (3) 


t—>0 


ou la fonction /(A) est donnée. 
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La résolution du probléme mathématique (1), (2), (3) permet d’obtenir la 
résolution du probléme fondamental de l’hydrodynamique, c’est-a-dire, de dé- 
terminer le mouvement non-stationnaire d’un liquide visqueux dans un domaine 
non-borné. 


En effet, au point de vue de I’hydrodynamique, les conditions (2) expriment 
que le liquide visqueux a la vitesse donnée tangente au bord du domaine. 


Evidemment, la premiére condition (2) exprime qu’en tout point P de la 
frontiére C la composante normale du vecteur de la vitesse est nulle, c’est-a-dire 
que le vecteur de la vitesse ¥ est tangent au bord du domaine. 


La seconde condition (2) exprime que la composante tangentielle du vecteur 
de la vitesse du fluide a une valeur donné en tout point P de la frontiére C. 
Dans le cas G(P, t)=0 les conditions (2) expriment que le liquide visqueux 
adhére au bord du domaine. Quant a la condition initiale nous pouvons remarquer 
que la fonction f(A) est déterminée a une constante prés par les composantes 
vy et v} du vecteur de la vitesse, données a priori au moment initial. En effet 
la fonction /(A) s’exprime par l’intégrale curviligne 


f(A) =f vj dx —v?dy + const, (3’) 


A,A 


A, étant un point arbitraire du domaine D, étant supposé que les intégrales 


fujdx—vdy=0, aa ae ee 
Cy 


La condition initiale (3) étant satisfaite, nous en pouvons déduire (démonstration 
donnée dans le présent travail) que 


Himes =o, Hime, = 
Donec, il est évident, que d’aprés la solution du probléme (1), (2), (3) et les 
inégalités (1‘) nous pouvons déterminer en tout point arbitraire A du domaine D 
et pour chaque moment ¢ de I’intervalle (0, T) le vecteur de la vitesse du fluide 
visqueux, dont les composantes sont connues au moment initial et au bord du 
domaine. 


2. Résolution du probléme aux limites 


Pour résoudre le probléme proposé (1), (2), (3) on admet les suppositions 
suivantes: 
I. Les lignes fermées C,, Cz, ..., C, ont la tangente continue en tout point, 


et l’angle que fait cette tangente avec une direction fixe satisfait a la condition 
de Holder suivante: 


|So9,| S const r60,, O<a<1 


ou dg, désigne l’angle que font les tangentes aux deux points arbitraires Q et 
l 
Q, de la frontiére C. 
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II. La fonction donnée g(t) est définie et continue dans l’intervalle 0</ST 
et satisfait aux inégalités: 


le () — g(4)| SA, 
0<m, = |g (|S M,. 


e<P 31, 


III. La fonction donnée G(P, ¢) est définie et continue en tout point PEC 
et pour 0<ft<T, posséde la dérivée par rapport a la variable ¢, verifiant la 


condition 
|G; (P,t) —G,(B,O|Skcren, O0< PS. 


IV. La fonction connue {(A), donnée par la formule (3’), est continue en 
tout point 4€D+C, posséde les dérivées deuxiémes continues, vérifiant la 
condition de Hélder dans le domaine D+C. Ces dérivées sont bornées dans le 
domaine D+C. 

On admet de méme: 

af (P) 
dnp 


Gi PAO): 


V. La fonction donnée (A, ¢) définie dans la région [A €D, ¢t€(0, T)]} est 
continue par rapport a la variable ¢ et vérifie la condition de Hélder par rapport 
au point A. Elle est bornée dans la région citée. 

Introduisons deux nouvelles fonctions inconnues ¢ (P, ¢) et 7(P, t) liées avec 
la fonction donnée g(¢) par l’éqalité: 


(Pt) + M=y(P, 2) - [eg (t) + My] (4) 


ot. M,>sup|g|, M est une constante positive que l’on va préciser plus tard. 
En s’appuyant sur notre travail [2] et [3], nous allons chercher la fonction 
inconnue p(A, ¢), vérifiant l’équation (1) et les conditions (2) et (3) sous la 
forme d’une somme des potentiels suivants: 


y(4,)=— A f [fot (A, t; B,t 1) | 5” oa be a en do, dt+y,(A,t)+ 
(5) 
+f foCuiaeuston 1) dlydt 4 if hai TSO, lg de 
0 C 
ou la fonction connue wo(A, 4) est donnée par l’expression: 
t 
(A,2) = [te “s a, (B) doy — + Lf, o(A,t; B,t)9(B,t)dopzdt (6) 
D’ 0 D 


ou D’ désigne un domaine obtenu par le prolongement du domaine D A V’intérieur 
des domaines limités par les courbes Cy, .. , C,. On conserve dans ce prolonge- 
ment les propriétés admises pour la ios re 
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D’aprés les propriétés limites du potentiel de simple et de double couche, 
discutées dans le travail [2] et en tenant compte des conditions (2), nous de- 
mandons que les fonctions inconnues (Q, 7), ¢(Q, 7) vérifient les équations 
intégrales 


t 
M,-x(P,t)—M—L(P, 2) , 
(P, 1) , tai aie 


1) dly dt — [fo (P, t: 0, 1) (OQ, 1) dlodt + 


et O(P, +, Boa) 


—xfu(Ps t)dt+ [fa 4° 4 (Q,1)dlodt 


6p ocAy op oAy d 
=: GUE =H, (UL 
OyR OxB 0X%R OYE YAP) 


t 


pa tases us WO (O qydlda 


dnp dng 


Pell NED 


(supposition faite que la fonction y(P, ¢) n’est pas nulle). 


op oAy oy oAy 


d IPG 
: = dop dt — Yo ue 
OVB OxXB OxB OVR 


dnp 


ap Gilt) 


Les équations (4), (5), (7) forment un systéme de quatre équations intégro- 
différentielles, non-linéaires, 4 forte singularité, avec les fonctions inconnues 
w(A), t), C(P, t), x(P, 0), u(P, t). Pour résoudre le systeme d’équations (4), (5), 
(7) étudions le systéme auxiliaire suivant: 


~ 


WA) = OAS te Un 0, we, 2, Cy le) 
v(A, t) =0,(A;t: u,v, w, 2,6, M) 
w (A, t) ZAR Ws U,V, Wee, Cot) 


2(Ast)= 0, (AS open 2,.C, 0) (8) 


—C(P,t)+ay(P,t fee PNA cy (ee t) Q,(P (P,t; u,v, wv, 2,0, uw) —4(P,t)M,+M 


— u(P,t) 4 [fale | H(@ ,t) diy dr=+9,(P, t; u,v, w, 2,6) 


a(P,t) = [C(P, 4) + Mi/[g@ + Mi] 
14* 
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oll nous avons pose: 


Q, = = ts lil ay [vz—uwwlpdopat+ 


Ox 
: 0 D , 
a ian val ff we! ow 
+ff °° 4(Q,1) diode + { | ene 1) dlodx + 2¥e 
0c 0c 


t 
a 4 r ow 
ie St = fuz— uwlpdopdt+ 
az: ah ay | Indes 
0 D 
t 


t 
+ ff 22 wl@vdigae+ ff B(F2)t0.0) digdr + 
D eC 


oy dng. 
t 
OS -=— | [| ae |\uz—uw|pdogdtt+ 
t t P 
+ ff 22n@odigact ff £(452\0(0.2) digae + Se 
) 0 


Cc 


t 
‘= : lai ae [vz—uw]|pdogdt+ 
a, 


t 
+ f [SP u@rdigact ff F(A Fup) (0.2) dl dx + “7% 
0c 


t 
— | [ou(Q.2)dlgdr + [ [[ oloz—uw]sdogdr — yo 
C 0 D 
t 


yD é d “( dw 

C S = = a x 

OR ay Ee {/ C(Q, t) dlodt 
0 


— 


Cc 
ae ff SEs )ter—wotadenar— 2 (4) + 2 rete. a 


Le systéme d’équations (8) se compose de 7 équations intégrales aux fortes 
singularités, avec 7 fonctions inconnues u(A,t), v(A, t), w(A, 2), 2(A, 0), C(P, A, 
w(P, t), 1(P, d). 

Nous allons résoudre ce systéme par la méthode topologique, en appliquant 
le théoréme de J. SCHAUDER relatif au point invariant de la transformation 
fonctionnelle [4], de la fagon indiquée par W. PorcorzELskI dans le travail [1]. 
Considerons donc un espace fonctionnel A composé de tous les systémes de 
7 fonctions réelles U[u, v, w, z,C, M4, x] déterminées dans le domaine [AED+C, 
0<tST] ou [PEC, 0<tST], ou ces fonctions sont toutes continues dans les 
domaines cités, et de plus, les fonctions ¢ et x possédent les dérivées par rapport 
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a l’arc, continues dans le domaine [PE C,0<t<T]. Les fonctions u, v, w, z, ¢, a a 
d¢/dl, dy/dl vérifient les inégalités: 

¥414t°|u(A, t)| <9, 44 t°|v(A,t)|< 0, taigtt* | w(A, t)| <0, 
ralat?** | 2(A,t)|< 00, P|(P,t)|< 00, | y(P,t)| < 00, t+" |w(P, 2)| << 00, (40) 


fo 


dc 6 
<= ’ 


F<, OWN OE OF) 7 S05 


A, étant un point a l’intérieur d’une des courbes C,,..., C,. On définit les 
opérations linéaires dans l’espace A par les relations 


[u, v, W, 2,C, MX] + [My, ys Mr, 245015 Ma X11 
=[u+m,v+%4,w W1,% Me 5G Cis b> XZ 1 | 


MU, Vv, wW,2,6,u,x|=|mu,mv,mwv,mz,ml,mu, my]. 
On définit la norme d’un point U|w, v, w, z,f, u, x] par la formule: 
|| 7 |] = sup lars O| (A, 1)}|-+sup [7474 | (A, 4)[]+-sup [1474 "| w (A, [J+ 
st suplraal’|2(4, 9] + supe (P,. 9] + sup[|x(P. al] + (44) 
+ sup [é°*”|u(P, é)|] + sup G Fal + sup G Fa | 


dee Av esd 0 0), yi 0 (A= maxi, A,)'du travail: (2). 
On définit la distance de deux points U et U’ par la formule 
6(U, U’) =|U — U'. (11’) 


L’espace A est linéaire et normé; en outre on peut démontrer qu’il est complet, 
c’est a dire que la condition de Cauchy est a la fois nécessaire et suffisante pour 
la convergence d’une suite arbitraire de points de l’espace / au sens de la norme 
(11). L’espace A est donc un espace de Banach. 

Considérons maintenant dans l’espace A l’ensemble E de tous les points 
Ulu, v, w, 2,0, u, x] Satisfaisant aux conditions: 


[u(A,)|=R, lo(4O|SR, “| w(4,)|SR, “| 2(A,O/SR, 


dl oe , ' = 
IC(P.)|Sa<M, ar | <a t” |u| See, O0< es S% Sos, | Sb 
; , BAGE WER SE fs 
(Pe. —C(P.A)| sofia’, |S G4) SO cay, (12) 


i |W(P, ) — 4 (BR, 8)| Srarbe,, 
In(P.) —4(P.t)| Salt 4 


p! ay (P, t) dx (F, ?) 
, dl dl 


| SK, Yop, 


ou B est un exposant positif arbitrairement choisi, mais tel qu’on ait BSa (« 
étant un exposant de la supposition I), 6’ est un exposant arbitrairement choisi 
de l’intervalle (3, 1). Les constantes positives R, Q,, 01, 02, 63, Mitty How Hes %es Os 
sont arbitrairement choisies, mais doivent vérifier les inégalités que nous allons 
03+ 03 M ' M 
ple ae Fees ct ey, 
sac Coe 


préciser plus tard. On demande seulement que uM, 
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L’ensemble E est évidemment fermé, puisque les fonctions limites des suites 
uniformément convergentes des fonctions “(A,#), v(A,t), w(A, 4), 2(4, 1), C (Pio; 
u(P, t), 7(P, t) vérifiant les conditions (12) vérifient aussi ces conditions. 
L’ensemble E est en outre convexe; en effet, si [u, v, w, 2,0, u, x] et [%,%, 
Wr) 215 C1» May Z4] Sont les deux points vérifiant les inégalités (12), alors les fonctions: 


(1—m)u+mu, (1—m)v+my, (1—m)w+mw,, (1— m) z+mz, 
(1—m)C+mh,  (t—m etm, U—maztmn, 


vérifient aussi ces conditions, si le nombre réel m varie dans lintervalle (0, 1); 
cela veut dire que tous points du segment rectiligne joignant les points U et U, 
dans l’ensemble FE, appartiennent aussi a cet ensemble. 

Transformons maintenant l’ensemble & en faisant correspondre a tout point 
U[u, v, w, z,€, u, x] de cet ensemble, le point Ulu, v,w,2z,C,u, %] déterminé par 
les relations fonctionnelles: 


u(A, t) = Aste) ye, Coy 
5(A, t) =Q,(A, t; u,v, @, 2,6, m), 
w (A, ?) — 0, (A, t; U,v,W,2,C,M), 
z(A,t) =0,(A,t; u,v, w, 2,6, "), (13) 


t 
ZO(P ive ey CRs t) dt=y(P,t)Q2;(P,t; u,v, w, 2,0, 4)—x7(P,) M,+M, 
0 


t 
ee (Uwe |e) 
ae eae 
0c 
SEM 


z(P.) =F 


f(Q, 1) dlodt =~ Q,(P, t; u,v, w, 2,0), 


Tt 


Nous allons démontrer que ces relations déterminent un point Ulu, 0, we, c #4, x| 
de l’espace A. 

En effet, d’aprés les études du travail [2] et [3], les quatre premiéres équations 
du systéme (13) fournissent les fonctions #, 0, @,%, continues dans la région 
fermée [D+ C, (0, T)], en traitant les valeurs limites de ces fonctions si A +P, 
comme les valeurs correspondant aux points P de la frontiére C. 

La 5*™"° équation du systéme (13) a la forme d’une équation intégrale de 
Volterra et posséde la solution unique de la forme: 


t 
: Sx%(P,s) 


= A —71 ds x ~ 
C(P,) =4(P,) Q,—277(P,t) fe = Q,dt + M(P, t) (14) 
0 
ou la fonction Q, est bornée et intégrable dans le domaine 


[D+C, 0S¢ST], et M(P,)=-M+4Myyt+nyfe  [M—M, yar. 
0 


Kemarquons d’abord, que tous les termes du membre droit de la 5!*™° équation 
du systéme (13) possédent la dérivée par rapport a la variable ¢, ce qui résulte 
des théorémes 5, 12 du travail [2], du théoréme 7 du travail [1], et de l’expression 
(23) du travail [3]. Il en résulte que la fonction €, définie par la 5%™¢ équation 
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du systéme (13) satisfait 4 la condition de Hélder par rapport a la variable ¢, 
avec l’exposant f’, admis pour la fonction 7. Remarquons de méme, que tous 
les termes du membre droit de la 5*™° équation du systéme (13) possédent la 
dérivée par rapport a l’are vérifiant la condition de Hélder par rapport au point 
P, avec l’exposant f_arbitrairement inférieur a l’exposant « (supp. 1). Cette 
propriété résulte pour la premiére intégrale 

t 

dw 

of dng (0; T) aly dt 


UisiG, 


du théoréme 5 du travail [3], et pour la seconde intégrale 


[folA,t; Q,t) W(Q, t) dlg dt 
0c 


des théorémes 9 et 10 du travail [2]. 
La preuve que l’intégrale 


t 
f ff w(A,t; B, ct) [vz—uw)pdozdt 
0D 


et la fonction donnée yw, possédent la méme propriété, résulte des théorémes de 
W. PoGoRZELSsKI, donnés dans le traivail [/|. En tenant compte de la méme 
supposition admise pour la fonction 7 (P, t) nous pouvons constater que la solution 
€ de la 5*™° équation du systéme (13) posséde la dérivée par rapport a l’arc, 
vérifiant la condition de Hélder avec l’exposant B=«. 

Etudions maintenant les propriétés de la solution de la 6*™° équation du 
systéme (13). Cette équation a la forme d’une équation intégrale de Volterra 
a singularité faible et posséde la solution unique continue de la forme: 


t 
iD eNO a | { RB Q, 1) Q, dlp dt (15) 
0c 


ot le noyau X(P, t; Q, t) est une fonction déterminée et la fonction donnée 2, 
est bornée et intégrable dans le domaine [D+C, 0S?#<T). 

Nous allons démontrer que tous les termes du membre droit de la 6"""* équation 
du systéme (13) satisfont a la condition de Hélder par rapport aux coordonnées 
du point P, avec l’exposant f. Pour la premiére intégrale 

t 


ta) d y dw 
a | [ wae. (O;%) dl, at 
G 


dmp , 
0 


cette propriété résulte du théoréme 9 du travail [3]. 


t 
af(d 

{ff a (Gar) [02 uw|pdopdt 

0 D 


posséde la méme propriété, ce que l’on peut démontrer de la méme facon que 
dans le théoréme 7 du travail [/]. 


L’intégrale 
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apo 
dnp 
pour la fonction f, ce que l’on peut montrer sans peine par le calcul classique. 
En s’appuyant sur la continuité évidente de la fonction “7, nous pouvons constater 
que l’intégrale du membre gauche de la 6° équation du systéme (13), comme 
le potentiel de double couche pour I’équation de chaleur, satisfait a la condition 
de Hélder avec l’exposant arbitrairement inférieur a l’exposant «. Nous con- 
statons donc que la fonction /# définie par la ome équation du systéme (13) 
satisfait 4 la condition de Hélder avec l’exposant f. 

La fonction x(P, ¢) de la 7°” équation du systéme (13) satisfait a la condition 
de Hdélder par rapport a la variable ¢ avec l’exposant f’, et posséde la dérivée 
par rapport a l’arc, satisfaisant 4 la condition de Hélder par rapport au point 
P avec!’exposant f, ce qui résulte des propriétés de la fonction ¢ et g. Evidemment 
7(P, t)>0 (remarque faite que €+M>0, ob M>o,=sup |¢|). Le point trans- 
formé Ulu, 0, ,z,C, #4, 7%] appartient donc a l’espace A. 

Cherchons maintenant la condition pour que le point transformé Ul [u, 0, @, 
gaye x] fasse partie de l’ensemble E. 

Remarquons d’abord que les fonctions wu, v, w, z, d’aprés les quatre premiéres 
équations du systéme (13), d’aprés les inégalités (12) et d’aprés les inégalités 
(7), (19), (32), (36), (53), (69) du travail [3] satisfont aux conditions: 


0 , sie : 
Ouant a la fonction donnée Bt ( cela résulte des suppositions admises 


— Re r 1 aed / a / ar ~ ~ 
\u Sa ie ee + hy (0g + %) O-* + ky (P70, + 26, + t%) + Ly, 


= Re 1 mary: , , , ss a 
ge Ky tI * + hy (Qg + xq) OO* + Rg (P40, +26, + 1%) +4, (16) 


[BS yy KaP ARS (0g BF preg BOB) pg OE og tt +O) +L, 


ames rR? bes py | EE , ae = 1 — Sey \Z ra — ~ ~ 
|z Se ‘s "+ k3(0.t* ’1 +x, tt 4s *) +k; (t *yyt+%+6,+0,t7%)+L., 


L,=sup|(yo)4|;  L.=sup|(4 yo) 4]; (17) 
Cherchons maintenant la limitation de la fonction £. D’aprés la 5*™° équation 
du systéme (13), d’aprés les inégalités (5), (39) du travail [3] et l’inégalité (43) 
du travail (2), nous aurons: 

IC] S @3 (4+ 2 @,t) sup|2,| + (— M+ M, @,) +2 03t(—M + Mj ,); (18) 
ou 

os rR , 1 ; 
sup |O,| = 5 Ketek + iyi upto (19) 
Ls = sup| po| 
(—M-+ My, 9s) étant la borne supérieure positive de la différence |M —(P,t) M,| 
1 : Pri M y M 

dans le cas ae (03+ @3)/2 > ™,’ 8% = Mi’ 
La fonction ¢ d’aprés les inégalités (12) de notre travail, et d’aprés les inégalités 


(65) du travail [2], et (16), (25), (38), (56) du travail [3], satisfait 4 la condition 
de Holder: 


| (P, t) —C(P, 4,)| S for og sup || + 2x3 sup |f| ¢+ My x3 + 
x3 sup|Q,| + 03 H'(Q;)} |? = t, |" (20) 
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ou nous avons posé: 
one ; 5 ae ee ? 
H’(Q,) = hy 8-4 05 + hy gy 0-2 + KB 4 1 (21) 
4 


L, étant le coéfficient de Hélder de la fonction y). En outre la fonction £ posséde 
la dérivée par rapport a l’arc qui, d’aprés les inégalités (12), les inégalités (64) 
du travail [2] et les inégalités 7, 52 du travail [3], posséde la limitation 


d 
Ea <nG,sup|E| t+ 95sup|5>|¢ +0, up| 25] + op sup | 47" 6/4 M8, (22) 


ou l’on a posé 


Cs ~ | ec R? 
sup fi Qs (ost sea tt tte hy (Es Lea) +55 Ky pitimw tT, 


d 
iTeY Ol (23) 


Lo = sup 


La dérivée 


S satisfait 4 la condition de Hélder suivante (d’aprés (12) et (58) 


du travail [2], (9), (50) du travail [3]): 


UIE 
|" ee aA et) |< {a % sup |C|t+ 0, H (Se Ae 


ae + 030, ¢- 


dl 
H (gr She, 


+ % sup |.Q,| + 03H (Q,) + M, a+ 0, sup| 37 ai 


ou l’on a posé 
dx ~ wo R* 1 == , 
H (57 2s)\= talon tm) O74 (Pa, btm) + oq we + Ts, (25) 


kw, Paves 33 ; d 
L; étant le coéfficient de Hélder de la fonction 77 yo- 


D’aprés la 6° équation du systéme (13), d’aprés les inégalités (12) de notre 
travail et les inégalités (68) du travail [3] et (18) du travail [2], nous aurons 


ji] < Kya sup || (26) 
ou l’on a posé 


t 
Ky =1+< sup [ [|X| digdz, 
oc 


a ; § R2 Bes fess 
sup |2¢| =k, [f 4 +27 0, + m) +5, Ket POMEL Lees (27) 
0 dy 
L,=sup, |; ae M, =sup|G|. 


La fonction 7 satisfait a la condition de Hélder par rapport au point P avec 
l’exposant f et d’aprés les inégalités (12) de ce travail et les inégalités 21 du 
travail [2] et (55), (67) du travail [3], nous aurons 
1-0 _ 
|Z(P, ) —H(B,, t)| Sa hy sup |a| ¢ ? ~ + 2 H(Q,) (28) 
ot l’on a posé 


a R2 , , 
H(Q¢) =hs [PA +t7% eo, +%] + res i447 406+ k, (29) 
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oO sd 
ou L% désigne le coéfficient de Hélder de la fonction => aE Ee. 


En s’appuyant 


sur les inégalités (12) et les suppositions IIT admises pour la fonction g nous 


aurons de méme: a+M M—o, 
ieee at : wile M, +M, % 
|x (P, t) = (e t)|S = asiA “4 (M,+ My) + hg (01 eae 
qz|_ ow 


BE \— mg+M, : 
ax(P.)  d(RO| Fs 


dl dl |= m+M,~ 


(30) 
(31) 
(32) 
(33) 


D’aprés les inégalités (12) et (16)—(33) nous pouvons affirmer, que l’ensemble 
transformé EF fera partie de l’ensemble £, si les constantes du probleme vérifient 


les inégalités: 
R2 ql , , , t ~ ~ 
et TE441* + hi (Og +H) T1-* + hy (T'*0,+7T6,+74,)+L,35R, 


R2 ; 
ro 1 [es ks (Teetxo; ob Daye) -- 
Ah, (TE 44 ET? eo, ETS EP OL eS 


Ka sup|Q, a5 
i, a 
wh,sup|iT * +o H(O\= x, 
(— M+ My 03) + (1+ % @5 T) @5 sup|25| +27 e,;T(-M +My) Sa, 


M+ 
M,+m, = = Qs, 


M, x3 + 703 (1+ 2 03 T) (sup |, | + M,) + 2x3 0, T +x sup|Q,| = 
+ 05H'(,) — 05 M(1+ 21) Sx, 


Hy (Mg +My) +hg(e,+M) a! 
(m,+™M,)? et 


Mabe ete Tee a. Lo eaeu le | + esup| 49 


mi 
Se 


M,- oe 0; S63, 


MX + 7 Hs 0, T + shes T + %, sup|Q,| + 


H(Q ,) +6 03 01 T + 030, T + 


Ces = 
+ dsup| “FP a + 0H (- t Blok. 
1 bs 3 I—o 
ee rr aes SS 
M,+m, 14 = M3, M,+M, = Qs- 


Etudions le systéme d’inégalités: 
PO Opes Kym (Le + M,) < 02; mt (Le + k,) < m3 
— be 1 
M + 03(L3 + My) < a; My +m, (0. + M) S@;; 


(x 03 ie us Te a M ae 15 = ; alas #4 (mg +M,) +hg(o,+M) = 
3 3) ( 3 1) set! aM 03 < 441; (mg+M,)? Sx 


r @3 (Ls + My) 4 oa ayia oe ~ 01 SOs; 


Gots + Halls + My) + OL, <%; MA a Hy Shs; ees 
§ 1 


Mouvement d’un liquide visqueux 207 


dont les membres gauches sont obtenus par la substitution T=0 dans les 
inégalités (34). 

Les quatre premiéres inégalités (35) seront satisfaites 4 la condition que les 
constantes arbitraires R, 05, x2, soient suffisamment grandes, et la derniére a la 
condition que 93 soit suffisamment petite. Les paires des autres inégalités (35) 
seront satisfaites simultanément, si les constantes arbitraires satisfont aux con- 
ditions 

Mz ee L; (36) 
et a la méme condition que la constante g, soit suffisamment grande et la constante 
03 soit suffisamment petite. 


Il s’ensuit que les inégalités (34) seront surement satisfaites s7 l’intervalle 
(0, T) de temps est suffisamment petit et les fonctions g(t) et f(A) sont choisies de 
facon que l’inégalité (36) soit satisfaite. 

Nous en concluons que les fonctions 4%, v, W, z, jo i, 7, Ssatisferont aux condi- 
tions (12) et l’ensemble transformé F fera partie de l’ensemble E. 


Nous démontrerons maintenant que la transformation (13) est continue dans 
Yespace A. Supposons que 4%, U,, Wy, 2n>Cn» Mn» Xn St la suite des éléments 
de l'ensemble £ qui tend vers le point U[u, v, w, z,¢, u, x]. Il en résulte que les 
suites {u,‘, {v,}, {w,}, {z,}, {C,}, {un} {y,} tendent uniformément vers les fonc- 
tions u, v, w, z,C, uw, x. D’aprés les propriétés connues des intégrales a singularité 
faible nous pouvons constater que ces intégrales tendent uniformément vers les 
intégrales analogues des fonctions limites. Les intégrales a singularité forte 


t 


a : d (Morte 
pen a daca talangle Oe 
0c 37) 
avi, ole=| ff A<2. ¢(Q,2) dig 
dé 


possédent la méme propriété, ce qui résulte des études des travaux [2], [3]. 

Nous en concluons, d’aprés les propriétés connues des solutions des équations 
de Volterra, que la suite des points U,,[,,, 0,5 Wn» 2n1Sn» Mn» Yn] transformée par 
le systéme (13), tend uniformément vers le point U. La transformation est donc 
continue. Il reste 4 démontrer que l’ensemble transformé est compact. Con- 
formément a la définition, l’ensemble E est compact, si de toute suite {U,,' de 
ses points on peut extraire une suite partielle orm: convergente au sens de la 
norme (11). 

Soit donc une suite arbitraire de points U,, ns Vers Orns Sons rine, Mane Nes de 


= : lS ee eS eee eT 
ensemble EF. Les fonctions 4,,,Um>@ms2m> Sm» Lms Xm» a a vérifient les 


inégalités (16), (18), (22), (26), (30), et (33). 

Considérons un cercle C, de centre A* et de rayon k,>L, A* étant un point 
4 l’intérieur d’une des courbes C; et L le diamétre du domaine embrassé par les 
courbes C, et désignons par D, le domaine embrassé par la frontiére C et le 
cercle C,. D’aprés la définition de la norme (11) et les inégalités récemment 
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citées, nous pouvons choisir le rayon R, tellement grand pour quon ait 


|Gnl< ze (38) 


si AC D—D,, 0<tST, m quelconque. 
Ensuite, au nombre « faisons correspondre un nombre positif ¢, pour qu’on ait 


|G.l<te (39) 
SACD, 0<tSi.. 


Le cercle C, et l'intervalle (0, ¢,) étant fixés, il est facile de démontrer que 


ee ahs Keenan dx 
toutes les fonctions %,,, Un, Wn» 2m Rigstyrac Ts, se ae, rs sont également continues 


et également bornées dans la région [A € D,, ¢, St T]. Donc, d’aprés le théoréme 
connu d’Arzela, on peut de la suite {U,,} extraire une suite partielle {U,,} con- 
vergente uniformément au sens habituel dans la région: [A€D,, 4,.StS<T]. 

Nous pouvons donc au nombre «¢ faire correspondre un nombre WN, tel, qu’on 
ait dans la région AED,, ¢,5t<T: 


5 (Oy — Un.) = Tin — Uiall < € (40) 
si m, s>N. En somme nous aurons d’aprés (38) et (39) 
6 (Op, Ve. U,,) = | a U,, | = é 


si AED, 0<tST, m, s>WN, ce qui confirme la convergence de la suite partielle 
{U,,,+ au sens de la norme (11) dans l’espace A. 


Toutes les conditions du théoreme de J. SCHAUDER étant satisfaites, nous 
pouvons constater que dans l’ensemble E il existe au moins un point U*[w*, v*, 
w*, 2*,¢*, u*, y*| invariant relativement a la transformation (13), c’est-a-dire 
vérifiant le systéme d’équations intégrales (8). 

Considérons maintenant la fonction yw*(A, ¢), déterminée dans la région 
[AE€D+C, 0StST] par la formule 


y*(4,)= of {feta B, 1) [v*(B, t) 2*(B, t) —u*(B, 1) w*(B, 1)] doy dt + 
mellinge Q, t) w*(Q, t) dlydt + (41) 


COAL, ES : 
off - ae C*(Q, t) dig dt + (A, i). 


On voit, d’aprés les équations (9), que les fonctions u*(A,t), v*(A,t), w*(A;a), 
z*(A, t) vérifient les égalités 


u*(A, t) sadn v*(A, t) = Halbtedh 
i y 


wh{A, i) = SEP 
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Les fonctions obtenues ¢*(P, ) et w*(P, t) vérifient donc le systéme d’équations 
t 
= C*(P,t)+4*(P,t) [ C#(P, ct) dt 
0 
= 4*(P, t) Q5(P, t; ut, v*, w*, 2*,C*, w*) + M— M, y* (P, t) 


SAS 5 _ Oy dy* dAy*t eéAy* 
= 7*(P, Sune ae ee 


— *(P, t) 4 < ff (eer 1) dlp dt 


=-9 a (P, t; u*, v*, w*, 2*,C*) — 10 6 Pats oy* as eAy* @aAy* cr) 
n Ox ov OX oy 


et de méme le systéme d’équations (7). En effet, la seconde équation du syst¢me 
(42) est obtenue par la différentiation de la seconde d’équations (7) par rapport 


a la variable ¢ et la relation est reversible puisque Lu =0, grace ala sup- 
Wp 


1 
position admise sus Let (20): 
P 
Enfin nous pouvons remarquer que la fonction y*(P, t) vérifie l’égalité 
C*(P, t)+M | 
+«(P, j= VU. 


Il en résulte que les fonctions w*(A, t), ¢*(P, t), u*(P, t), ~*(P, t) forment une 
solution du systéme d’équations (4), (5), (7). 

Nous allons montrer que la fonction obtenue w*(A, ¢) est une solution de 
Péquation (1) et qu'elle vérifie les conditions limites (2) et (3). En effet, d’aprés 
l’équation (5), nous pouvons constater que la fonction p*(A, ¢) vérifie la condition 
initiale (3). Kemarquons maintenant que toutes les intégrales des membres 
droits des quatre premiéres équations du systéme (9) satisfont a la condition 
de Hélder par rapport aux coordonnées du point A (voir [2] et [3]). Nous en 
tirons que la fonction (v*z* —u*w*) vérifie la condition de Hélder dans tout 
domaine borné D* cD, et par conséquent la fonction w*(A, #) admet les secondes 
dérivées en tout point intérieur A €D, alors elle vérifie l’équation donnée (1). 

Il est aussi évident, d’aprés les deux premieres inégalités (16), que les fonctions 
u* et v* possédent la propriété limite 


lim w*(A, t) = f(A), lim v*(A, t) = fy (A) 
t— t—> 


puisque les derivées premiéres de l’intégrale de Poisson-Weierstrass possédent 


cette propriété bien connue et les autres intégrales dans 2, et Q, tendent vers 
zéro sit—>0O. 


Théoréme. Si les fonctions qui figurent dans ]’équation (1) et les conditions 
(2), (3), vérifient les suppositions II, III, IV, et sont choisies de fagon que l’in- 
égalité (36) soit satisfaite, si le contour C, limitant le domaine non-borné D vérifie 
la condition I et si l’intervalle de temps est suffisamment petit, alors il existe 
au moins une fonction y(A, ¢), qui satisfait a l’équation (1) en tout point A €D 
pour 0<?<T, qui vérifie la condition limite (2) en tout point PeC, pour 


: i cake alive ie. ys any 7) 
0<t<T, qui vérifie la condition initiale (3) et dont les dérivées = ' 2 tendent 
vers les valeurs limites /,(A) et f(A) conformément, si¢—>0.  ~ 
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Au point de vue de I’hydrodynamique nous avons donc démontré qu'il existe 
dans le domaine non-borné D et dans I’intervalle de temps suffisamment petit 
un mouvement d’un liquide visqueux, incompressible, dont les composantes du 
vecteur de vitesse sont données a priori au moment initial et au bord du domaine. 
Il faut remarquer, que la condition exprimée par Vinégalite (30), epnceman i 
fonctions g(t) et f(A) arbitrairement choisies, ne diminue pas le caractere general 
du probléme de l’hydrodynamique. Il faut remarquer aussi que la solution 
précédente concerne de méme le cas du domaine borné, multiplement connexe. 


Corollaire. En s’appuyant sur le théoréme précédent et l’expression (41) 
nous pouvons déduire que 


9,,0 a 0 
lim ee os | = ace Pe (44) 
t+>0\ 0x oy Ox oy 
c’est a dire 
lim rot, = rot,v° (44’) 
t—>0 


En effet, d’aprés les expressions (1’) et (41) nous aurons: 


rot creeps! 


al Aw(A,t; B, 1) [v*(B, t) *(B, t) — u*(B, 1) w*(B, 1) ]dogdt+ (45) 


t 


_ [ fidel4.n0 Np4(Q,1) dlgdt— [ [(4,650,7) w*(Q, 2) dlodt—Ay A, 


0c Cc 


tt) 
En s’appuyant sur les inégalités (17), (25), (56) du travail [3] (remarque faite 
que V/ = AV, U/ = AU) nous pouvons constater que les intégrales dans l’expression 
(45) tendent vers zéro si t—>0. 
La fonction 
TAB 


Ay(A,t) =A as ff re * FB) dog 


; (46) 
“yA eff (A, t; B, t) 9(B, 0) doy dr 


se compose de deux intégrales dont la seconde posséde la propriété précédemment 
citée, la premiére a la forme du laplacien de l’intégrale de Poisson-Weierstrass. 
U1 est facile de démontrer d’aprés les suppositions admises pour la fonction f(A) 
que: 


VAB 


im | alll 7 1B B) da] = ANA), (47) 


Nous en concluons enfin que 


lim rot,¥ = lim(— Ay) = — AT (A) == Tot, 0" (48) 


t—>0 t—-0 


d’ou l'on obtient l’égalité (44’) ou (44). 
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The Influence of Edges and Corners 
on Potential Functions of Surface Layers 


RouF LEIS 


Communicated by C. MULLER 


Let S be an analytic, orientable surface in three-dimensional space; let p 
and g be two points, 7(f, q)=|p— | the distance between p and gq, and e(g) 
a function defined for all géS. By 4/én, we denote the derivative in the normal 
direction relative to g€S. The potential function of the N-fold layer @(g) is 
then given by 


Uy () = f o(a) oe. apy ts 
S 


The functions Uy(p) are analytic in the neighborhood of all points # not 
lying on S. Uy(f) jumps when # crosses S at an interior point. Papers by 
LiapounorFF, Poincare and E. Scumipt (cf. [7]) deal with these jump relations. 
They are explicitly stated for potential functions and their derivatives of arbitrary 
order in a paper by C. MULLER [2]. While these papers deal only with the 
behavior of potential functions near an interior point of the surface, for many 
problems it is of interest to know the behavior of potential functions and their 
derivatives at points # on the boundary C of the surface S. This behavior will 
be investigated here. We shall see that boundary points are singularities of these 
potential functions. The degree of the singularity will be given explicitly. 

So as to avoid burdening the exposition with differentiability assumptions, 
let us assume for the present that the surface S and the layer function S are 
analytic and that the boundary curve C is piecewise analytic. In the last section, 
this assumption will be weakened by requiring only differentiability of finite order. 

The formalism of differential geometry needed will be given in the first 
section, where we follow the notations of C. MULLER [2]. The second section 
deals with the functions U), U, and their gradients. In the third section, the 
singular behavior of the functions Uy and V Uy will be reduced to line integrals 
over the boundary C. These integrals will be discussed in Section 4. Finally, 
Section 5 weakens the requirements of differentiability. 


1. Introduction 


Let the orientable surface S$ be bounded by the piecewise analytic Jordan 


curve C (consisting of corners and regular points). Let S have the following 
properties: 


Potentials of surface layers 243 


1. There exists a coordinate system for every point P of S (a closed domain). 
2. P is the origin of the coordinate system. 


3. A constant c>0 exists such that the part of S lying in the interior of the 
sphere (x7)?+ (x*)?-+ (x%)?<c can be represented in the form x8=@(x1, x2). 
(x1, x?) is an analytic function for all (x1)2-++ (x?)2<c. 


4. The derivatives aoe @), and ®,). vanish at the origin. 


Thus the «°-axis is normal to S at P; the x1- and x2-axes are the directions 
of principal curvature. For brevity we call this coordinate system a tangent- 
normal system. 


Let e; be the unit vectors in the direction of the coordinate axes; 
(1.1) ey arena Dia 47) 6,7 (4) 2x") 


then represents a point of the surface. { denotes a vector with the components 
f, 72, #8. At the origin we get 


(1.2) ihe sens fiz = Pre es 0. 


In (4.2) and in the following equations, = denotes “equal at the origin’ (at P). 
The range of Greek subscripts is 1, 2; of Roman, 1, 2, 3. 


Now we introduce a new coordinate system, which we shall call a u-system, 
by the formula 


(1.3) t= {(4, vw) + n (4, u?) 


where 1 denotes the unit normal to the surface. Then we have 


(1.4) = bi, 

and, with 

(1.5) sea a Ly ee Vey — (iGatjeys 

also 
‘Sig (Uh. Lp) oe (Fu — 48 ibs i) (F1, eh Les, To) 

(1.6) 


Sig 206 Ly (U8) oes 


Let 2H=hk,+k, be the mean curvature and K=k,k, the Gaussian curvature. 
Then, with 


(1.7) Vig Ry Liga, Ly. == 05 We Og 
we have 
(1.8) C(t (@ PR) ype 240 #) Le. 


Equation (1.8) is a tensor relation valid at the point P. However, the point 

P is in no way distinguished among surface points; rather, at an arbitrary point 

of the surface there exists a tangent-normal system and, consequently, a #-system 
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in which equation (1.8) is valid. Thus (1.8) is valid at every point, and we get 
(1.9) Sur — (14— (u®)? K) Vyy— 243 (1 — uv H) ies 


and 


(1.10) £33 = (js 2s) =? = 1; lsu = (Lista) = (tlh. — # Ly fe) =0- 


Let g stand for det g;, and y for det y,,,; then g/y is an invariant of the surface 
(w3 fixed), and we get 


(1.11) g=y(1— 20H + (uw)? K)? 
or, with 
(1.12) G=1—2H + (v’)*K, 
simply 
(1.13) so Ge 
The Delta operator 

(a, a \2 a \2 
ely - (sx) | (ea) . (33) 
is given in the u-system by the formula 

bead oa 3R oo 
(13) Aix qe gaelSaiiagg 
or 
me oe pa ae, TOE RL: é 

se) Bigs G Vy eur ole G 6u3~ éus* 


and we write = 

(48) A= WS ae (40-75 soe Wr gee) 
y Out ow” Vy our ou” 

Thus we get 


Lemma 1. For sufficiently small v3, 


aa : 


| . 
y ous eur 


From Lemma 1, for u?=0 and for harmonic functions with 


14 de Oe 

(1.19) oor sanea on S 

it follows that 

1.20 ae e) a 
heed) AU +(5,-2H*\u=0 
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In analogy to Lemma 1, by differention with respect to u® we get 


Lemma 2. Let the function U be harmonic in the neighborhood of a point P. 
Then the derivatives of U satisfy the recurrence relations 


0={3 (9) aya +(e) — amor o (Sy an (Myla e. 


j=0 


The Nabla operator 


0 a) 7] 

(1.21) y= Pia rea 
is given in the u-system by the formula 
(1.22) VU =¢"'4,U, =" fj, +e n,) Uj, +n Up. 
For small values of u? we have the development: 

v A — 2 7 ro) 
(1.23) eye + 0 m,) = eM (88 — W LE) f= fe > Off. 

j=0 
With the notation 
(1.24) Vi= OS frgar (P= fe gor 
: i (7) 1” Baye 0 eli aul |? 

follows 


Lemma 3. For sufficiently small u?, we have 


Finally, we wish to derive from Gauss’ theorem certain formulas which we 
shall need later on. Let Y be a continuously differentiable function defined 
on S. Then we get for a sufficiently small surface S’ with the boundary C’ 


[UG A4S=f Wi Fu ly du ae 
? 


(1:25) Ore ey a. Paleeg ae (Ole ea at 
7a =| Paz (Qf) Ty Vy) der du + [Vr F fw(Qhe — Qi #) de. 
Sf Cc’ 
With 
f a oe) = 
(1.26) 1 = fp (Qs 7 — OG) #) Vy, Uae ae (OG) Tn V7): 


we find, after summation over S’, 


Lemma 4. Let Y be a continuously differentiable function defined on S. Then 
we have 


SV, fZ4S=—Jf AvaSt+f Fnu,de. 
s s C 
For ;=0 we get from Lemma 4 


Lemma 4*. Let Y% be a continuously differentiable function defined in S. 
Then we have 


(Vj fdS=—2fnuH JdS+fu fade, 
§ s C 


where Ny 1s a vector normal to C and orthogonal to n. 
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Let % and ¥ be twice continuously differentiable functions defined on S, 


Then we have 


ee AN Bh: 2 (Vy Sti L,) due du? 


- elke aly i tS L,, (St) w — St} w) de 
(e277) ; 
> =f 24 Gleam) Vy Z(Sij @ — Si} w) de — 
— [2 poy lStt it — Sit) ae. 
ie 

Since 
(1.28) S(i) = Vy (sje Siw), Sy a 3(j) Te» VL=Lf, ia 
we obtain 


Lemma 5. Let the functions H and LF be twice continuously differentiable 
on S. Then we have 


EER Tee > aut (4V,£ —LV.H) 3,dc, 
or, in particular, 


Lemma 5*. Let the functions H and L be twice continuously differentiable 
on S. Then we have 


[HA LaAS=S LA HAS +f (HmyVYL—LuyVH) dc 
S S Cc 


2. The behavior of potential functions arising from single and double layers 


The potential function of the single layer o(g), namely, 
2.4 Uy (p) = ae Se: 
(2.1) 0)= fe 76.9 


is continuous everywhere [1]. The behavior of the potential function of the 
double layer j (q), 


(2.2) TiO) = ful) en, vipa eS 


when # approaches an interior point of the surface S is characterized by the 
jump relations [7], [2] 
UO, lext. = a 
et ge lili: ane he 
(2.3) 


Chiat = — 
lint, ant fit an a Ix 


The integral 

2.4 ee ae 

(2.4) fr@ pat Se 
s 


PES 
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represents a function continuous for all #€ S. Thus U,(f) is not continuous as 
p crosses S, but it is bounded for all #. 


The gradients of these functions become singular, however, when approaches 
a point of the pounaty C. To show this, we form, according to Lemma 3, the 
quantity 


a a1 Sart 1 
Ds) Vu,= ~ fea ise fen tas fe%-4s. 
s s 
From this formula, according to Lemma 4*, we get 


Ones Al 1 
(2.6) VU,=— fons dS+ f Wyo) taS+2 [no tas— f motac. 
iS § § c 


The line integral in (2.6) becomes logarithmically singular [3], while the other 
surface integrals are regular. The singularity can be described by the formula [3] 


1 
(2.7) [mee t-de=—2n9(b) 0(v)g|z— | +01) as |z—y] +0, 
C 
where x approaches the regular point ) of C. Thus we get 


Theorem 1. The potential function Uy(p) is regular when p approaches the 
boundary C, while VU, (p) becomes logarithmically singular as stated in (2.6) and (2.7). 


We shall now discuss V U, (#) and, in a similar way, get 
= 0 1 _ a 4 
—VU, fu ae Veen funyatas+ 
s s 


C= al 1 
+ fulz cast funtas. 
Ss Ss 


(2.8) 


We discuss the resulting surface integrals individually. According to (1.20) and 
Lemma 5*, we < 


Be ttomer sss Sst 


(2.9) 2 fans 1a5— fdynn tas 
S) 


— pin, Vv) == — (119 Vo) wtp dc. 
i Y if 


The line integral has a singularity like 1/r as p approaches C [3]; explicitly, 


[fem Pa) } = 5 (to Ma) wen} de 
C 

2a(r—y) 
le—»]? 
when x approaches the regular point y of C. f(s) and 6(s) are the normal vectors 
of C, x is the curvature; yj Vo=aV. 


(2.10) 
— (wnx(§a) — 2(9%) wn) lg|y—y| +01 


=n 
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According to Lemma 4*, for the second integral in (2.8) we get the formula 


oA on 4 Bp A 
=— fom 2 tas—2 fntys Fast [mua 7d. 
s s é: 
The line integral in (2.11) also has a singularity like 1/7; explicitly [3], 


furan 7 8H) Tee 70) (5) Ig |x — |} + 


+0(4) as |r—y|—O. 
Finally, for the third integral in (2.8), according to Lemma 4, we have 


(2.13) fuvtas=— [Cin Pas— fupradst | uprede. 
S Ss S: C 


(2.12) 


Again the line integral is logarithmically singular: 


(2.14) f ny 7 de = — 2p mglg|x — | + O(1) as |xr—y|—O. 
C 
Thus we get 
Theorem 2. The potential function U, is bounded as p approaches the boundary 
C, while VU, has a singularity like 1/R (R being the shortest distance between p 
and C) as stated in (2.8—14). 


3. The behavior of potential functions of multiple layers 


In the last section we discussed the behavior of potential functions of single 
and double layers and their derivatives. Now we wish to extend these discussions 
to the potential functions of N-fold layers, 


py) A= fel 2) (re) Way Se 


To do this, first of all we wish to express the operator (@/@n)* through lower 
derivatives. Making the assumption 


OI\N St ase: 
(3.2) (sn) =OVt QS, 
by Lemma 2 we get the recursion formula 
u =f ~>(N+1 = N 
(3.3) Nya 2H(N +1) Ofrys — 2K ( ; } Qu 2G ) nv 


with the initial values 
(3.4) Qk=1, M=0, 10, a4, 


It follows from (3.3) and (3.4) that the operators Q!y and Q,} n41 contain deri- 
vatives of maximum order 2N and that the operators Q3y_, and Q3, contain 
derivatives of maximum order 2N — 2. 
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Now we form the adjoint operators OX to Qk according to 
(3.5) NAD LAS = | LOCALS (XL, S) de 
S s G 


where # and # denote sufficiently differentiable functions. Since the operators 
Ay _,; are self-adjoint, for ®{ we get the recursion formula 


N 
6 DX = ree = N-+1 HK == Nh L 

(3.6) N+e= 2(N +1) Ohi 2( 3 )OKK d(5) Of 49- 

with the initial values 


(3.7) C1, 8 (DPS OE Ot S01 OFS 4. 


According to Lemma 5, for the X we get 


N+2=2(N +1) Xt (HH, L) = 2 (NF) XaKoe, L)— 


(3.8) Nv 
— D1 Le MoO L — (Of L) MH] 8,_7+ X$ (Aw, D)} 


with the initial values 
(3.9) Xo = Xi = X2 = X?=0. 


From (3.8) and (3.9) it follows that Xiy(%, VY) and Xiyii(%, Y) contain 
derivatives of maximum order 2 N —1 with reference to Y, and that X}y_,(%, ¥) 
and X}y(#, Y) contain derivatives of maximum order 2N —3 with reference 
to LY. Xin(%, #) and Xoyii(H#, L) contain derivatives of maximum order 
2N —1 with reference to #. 


Thus we get 


Lemma 6. Let the function H be analytic and the function # be harmonic. 
Then there exist operators ee NK, Xx, defined by (3.3, 6, 8), such that 


. [ele (avrot f)@ 
an 


fHQRLAS= if LD Ha Safi Xi (L,H) dc. 
§ § C 
It follows from Lemma 6 that 


EM eee Oe et andl “pela 
we fels) ‘as feQhtas+ fo as 
& Ss S 


1 a 1 
(3.10) ae ve) aS+ [ (Phe) sas + 
S) S 


+ [ Xk(e, t)ac+ { Xie, sr 7) ae. 
Cc 5 


Thus we have expressed the function Uy through functions of the type U, and U, 
and boundary integrals. The line integrals in (3.10) therefore describe the singular 
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behavior of the functions Uy(p) as p approaches the boundary C. We shall 
discuss the singular behavior of these integrals in the next section. 


For the gradient VUy, according to Lemmas 3 and 4 we get the formula 


yas — 


é 
on 


(3.11) = forte)” : dS — S f Cx-s0)( 


7=0 § 
at >| frelan) $45— fvelgu) +44. 
HM com 


The surface integrals in (3.11) are again functions of the form Ux. By means 
of (3.10) we can express their singular behavior in terms of line integrals which 
we shall discuss in the next paragraph. 


4. The singular behavior of potential functions arising from one-dimensional 
layers 

In the last section we expressed the singular behavior of potential functions 

of surface layers through line integrals over the boundary C of the surface S. 

We shall now discuss these integrals. As the derivatives with respect to all 

coordinates u‘ occur in the integrands, we shall discuss in particular the integrals 


(4.1) Vw= f o(z) tae 

and ‘i 

(4.2) Wy =f ofa) tac, 
Cc 


where c denotes the arc length and n a vector normal to C. 
k 
We assumed C to be piecewise analytic, so let us assume that C= )'C,, 
i=1 
where the C; are analytic curves. Through partial integration, it follows from 


(4.1) that 

(4.3) Var 3 | olay eee afte [Gey 4 ac} 
i Ci 

or 

Sa a ne p[SO NED Be NN =F 43 w(far 4 

(4.4) Vy >| 2, 1) Ise) alee) v |0 eo) {Ge eae 


Thus we get [3] 


Lemma 7. The function Vy (p) is logarithmically singular when p approaches C. 
For N21 there are additional singularities at the corners of C of the order R~N 
(R being the distance trom p to a corner) 


Potentials of surface layers 221 


To discuss the behavior of the functions Wy, we represent the curve C in 
the neighborhood of an interior point P in a tangent-normal coordinate system 
in the following way: 

1. P is the origin of the coordinate system. 


2. A constant t>0 exists such that the part of C which lies in the interior 
of the sphere x?+ y2+ z2<1 can be represented in the form y=f(x), z=g(x). 
The functions f and g are analytic when x?<T. 

3. The derivatives /’ and g’ vanish at the origin. 

4. In P the direction of the normal vector n is the direction of the z-axis. 


If P is a corner point of C, we can con- 
tinue the analytic pieces of C over P, so that 
P again may be considered as an interior 
point (for each segment separately). 

To discuss the function W, as Q (x9, Vo, Zo 
approaches the point P(0, 0, 0), we proceed 
in the following way (see Fig. 1): Let Q be 
sufficiently near to P. Then choose a point 
P,(c=¢) so that Q has the coordinates (0, y,, z,) in the tangent-normal system to 
P,. Let P lie in a t,-neighborhood of P,(t,< 1). Then we can write W, in the form 


Fig. 1 


Ty 


a : Dy 
(4.5) W,=We+W* with W; =a} 05 —-de 
The function W,* is bounded; for W? we get ; 

wejsal f Bax =o [ eth aly—m) Hi=al aa| 

ee ot  VAF =) +e—4) 
(4.6) ‘ Pek 
ax 
“ade 


|= [R= ly] y= | +1) +e (eal Hla) 
y R| yvR(v+R) yR(y+hk) 


== O(4) jor 0-0, 


(4.8) 


and, by induction, it follows that 


(4.9) a ae |=0 (art) 
Thus we get 
(410) Wwe =0( f Fx)=(3): 


o being the distance of Q from the curve C (in the normal plane). The example 
of the straight line shows that in general the estimate in (4.10) cannot be improved. 
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In a similar way, in the general case (V=1) we get the result 


(4.11) wie = O( [ purer) = (Gy) 
Thus we have 


Lemma 8. The function Wy(p) has a singularity like o ™ as p approaches 
the curve C. @ is the distance from p to the curve C taken in the normal plane of 
the curve. 

We can now describe the singular behavior of the functions 


(4.12) Uy = f oz) as 
S 


Since the line integrals in (3.10) contain derivatives of 1/r of maximum order 
N —1, we get 


Theorem 3. The potential functions 


N 
Ux) = f o(s) Ld Sayin eA 
s 
have a singularity like Rt~N as p approaches the curve C (R, being the distance 
from p to the curve C taken in the normal plane of the curve). Besides this, at the 
corners there ave point singularities of order R$“ (R, being the distance from to 
the corner). 


Similarly, from (3.11), we get 


Theorem 4. The potential functions 


ae, a V f oz) +as N=1 
S 
have a singularity like Ry“ as p approaches the curve C. At the corners there are 
point singularities of order Ry. 


5. The reduction of the differentiability assumptions 


In order to simplify the presentation, so far we have assumed that the surface 
S, the boundary curve C and the layer function @ are analytic. However, 
analyticity was not used. Actually, it is sufficient to assume differentiability of 
finite order. In discussing the behavior of the functions Uy at interior points 
of S these assumptions are stated in [2]. Thus we need consider only boundary 
points. In discussing the functions Uy, only derivatives of order N —1 occurred. 
Thus we assume C, S and @ to be (N—1)-times differentiable. To be able to 
introduce the tangent-normal coordinate system, however, C has to be at least 
twice (or once Hélder-continuously) differentiable. The layer function @ has at 


least to be bounded. For the discussion of VU, the derivatives of order N should 
exist. 
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Singular Perturbations of Eigenvalue Problems 


W. A. HARRIS, JR. 


Communicated by A. ERDELYI 


1. Introduction 
This paper deals with the eigenvalue problem, 


(1.1) Lly, a] = Sy — AQ +B a} y=AROy, 
(1.2) s[y, a] = M(x) y(a, x) + Nla) y(b,a) =0, 


where A, B, R, M, and N are square matrices of order n-++-m, y is a vector of 
order n--m, « is a large parameter, and —coo<aXt<b<-+oo. The matrices 
B(t,«), M(«), and N(«) have convergent representations of the form 


B(t, 0) => B,(t) a7 
j=0 
(1.3) MEG Vw pt ie 
i= 
Nw) = Na, 
j=0 


for sufficiently large «, and all functions are infinitely differentiable with respect 
WO) Uy (SUS: 

The character of this problem as «-—>co depends on the matrix A, and we 
shall assume that A has the partitioning 


(1.4) A=(? a 
Ag, Age 


where Ay, is a square matrix of order m which is nonsingular for a<t<b. The 
problem is to investigate the spectrum of L[y, «], s[y, «] in its dependence on 
the parameter a. 


If we formally let «—>oo, that is, multiply (1.1) by ss 1 
Oy tears Ih 

are identity matrices of order 2 and m respectively, and let «—>oo, we get the 

degenerate differential system which can be written, after suitable partitioning, 


where J,, and J,, 


(1.5) Sapna (By Ps ie Bis Yo) =A(Riy i + Rye Vp) — (Aaa Vi a Ano Walia Oe 
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If Ay, is nonsingular, the last equation may be solved for y, in terms of ¥,, 
and y, may be inserted in the first equation to obtain a differential system of 
lower order. 

If the eigenvalues of (1.1), (1.2) have finite limits as «oo, one would expect 
that the limits are eigenvalues for the degenerate differential system (1.5) with 
appropriate boundary conditions related to s[y, oo]. Since the degenerate 
differential system together with appropriate boundary conditions is essentially 
an eigenvalue problem of lower order, one speaks of singular perturbations. 

These problems have a close connection with singular perturbations of 
boundary problems. There, the corresponding problem is the relationship of 
the solution of a boundary problem of the type 


Ly, o] = 2 — (2A +B) y=0, 
s[y,4] =My (a,x) +N y(b, a) =c(a) 


as a->oo to the solutions of a related degenerate problem obtained by setting 
a=oo formally, where as before A has the same structure. The boundary problem 
has been treated by many authors, especially W. WAsow [13], G. CHAPIN, Jr. [2], 
and the author [4]. The two problems have similar features, and the methods 
for the proofs are closely related. 

A typical feature of both problems is that for «<oo we have a differential 
system of order n+ m and a boundary form of rank +m, but in the degenerate 
case the differential system is essentially of order m while the boundary form 
is of rank »+m. However, if the limit as «—>co of a solution of the boundary 
problem exists, it is in general a solution of the degenerate differential system 
and ” appropriate degenerate boundary conditions. The determination of the 
appropriate boundary conditions is quite similar for both problems. 

It will be shown under suitable restrictions that corresponding to each simple 
eigenvalue A, of the degenerate (or unperturbed) eigenvalue problem there exists 
a uniquely determined eigenvalue of the complete eigenvalue problem, A(q), 
which has an asymptotic expansion of the form 


Aa) = Ay + aay 
Further, an eigenvector y (¢, «) corresponding to A(«) may be determined so that 
1 
y(t, &) = Yo (t) + i ian 


on the open interval a<t<b and uniformly on any closed sub-interval, where 
jo (¢) is an eigenvector corresponding to A) for the degenerate problem. 

It is not apparent that this method will yield all the eigenvalues of the 
complete problem. However, it will be shown that for a class of self-adjoint 
problems this procedure characterizes the asymptotic behavior of all the eigen- 
values of the complete problem. 

Similar perturbation problems in which the eigenvalues can be obtained as 
asymptotic series have been treated by T. Kato [5], J. Moser [7], and V. 
KRAMER [6]. 
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2. Preliminary Transformations 
We shall use the symbol 0 indiscriminately for the zero matrix of any di- 
mension, J for an identity matrix of any order, and J, for an identity matrix 
of order k when we wish to emphasize the order. We shall use the symbol diag 


(A: B) to represent the matrix ie a): All partitionings that occur are such 


that the diagonal blocks are square. The norm of a matrix A is denoted by || 
and is the maximum of the magnitudes of the elements of A. 

Hypothesis 1. 

(i) The matrices M(«) and N(a) have convergent representations of the form 
(1.3) when a ts sufficiently large. 

(ii) The matrix B(t,«) has a uniformly convergent representation of the form 
(1.3) when « is sufficiently large, where the elements of B;(t) are infinitely differen- 
tiable on the finite closed interval [a, b}. 

(i11) The elements of A(t) and R(t) are infinitely differentiable on [a, b]. 

(iv) A(t) has the partitioning (1.4). 

(v) There exists a square matrix Wy ) of order m whose elements are infinitely 
differentiable on [a,b] such that T33' (t) Ay» (t) Tr. (t)=(6;; B®), 


Pit) = Bo) =--- =8,, (4) $B, = one Hehehe Bao) i) sep 
Re{f,} < Re{f,} <-- <Re(f,) <0<Re{Bys}S---<RefB,}, a StS. 


Since Re {8,;'=-0, A,. is nonsingular, and the following transformation may 
be made: 


ZF 10 
(2.1) y (é,a) = | z(t, a) = 7 (2) 2(¢, a), 
Ty To. 
where 7, is the matrix given in H 1—(v) and T,,= — Agq' A,,. Under this trans- 


formation (1.1) becomes 
az 0) O = D 
——ly +B +aROhe 
at | é 6,89) : 


(2.2) . = 
B(t, «) = B,(t) a7 


The degenerate differential system corresponding to (2.2) may be written 
(after suitable partitioning) 


1 Ex 
a cat i oda Oy ee, AO eRe (t)} 2 


= {(0,;5; (0 eee 


or ea 


az EN 
(2.3) Fi = {Ba xdl) +AR, (t)} % 
O= 2, 
where 8 7 shi hg 
By = (om ssa Ry =( Ry 4 
Bo, a1 Bo, 22 Faye 


na 


Bou Ra are of order n, and By Foy ae are of order m. 
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Lemma 1. Let Hypothesis 1 hold. To every positive integer o and A>0O, for 
sufficiently large x, “>, there exists a fundamental matrix of the differential 
system (2.2) of the form 


o 


(2.4) L txcvA), = Px (t, A) See = Zor (8, A, 2)) E (00), 


where 


B(,0) =(8, jexp{a Sol at}), ae Sesh si aye a 
20414, a)I|< co, 


c, independent of t,x, A in aXt<b, a>a, |A|<A. Further Z,(t, A) =diag (Z}): 
Bes Dy ntyaa)y 24; $5. Of OVE? Ty iy — Ti o(Tg—=0), 1222,. det Z;;--0, 724, 
astsd, |A|<A, and in particular Z{, is of order n and a fundamental matrix 
for the degenerate differential system, 


d 
dt 


Tee {Bs (0) a5 AR (t ba 


Note. This method of describing the asymptotic behavior of solutions of 
differential systems by divergent power series is now standard. The investi- 
gations of TuRRITTIN [1/2] are especially appropriate, but since we have an 
additional parameter /, and since we need the dependence on A of the constructed 
solutions, we give the following proof of this lemma. 


Proof of Lemma 1. We first construct a formal fundamental matrix and then 
prove that this formal matrix is asymptotic to a ¢vue fundamental matrix. 


The equation under consideration is of the form 0 gee Leis) Cir?) 
where Cy=diag (0:(6;;B))), Ci, a)=Bo+4R, C,Q=B;,(, or Sub- 
stitute a founpl series Z(t;a,A)= 2, (t, A) + —Z, (t, A)+ LEG, «), where 
x a) ) = exp {a i Cet 7) an}. (Since Cy(¢) is a diagonal matrix, we have 
£ Et, See nea ; We obtain 


Zo +4 Z,++--)aCyE+(Zo+ EZ +---)E=(aCot Gt Cot) Zot) E 


or 
ZyCyo=CoZo; 
(2.5) Z,Cy+ZM=QO4+GQZ, 
Za 2g Co Oy 2, Calg ele: 


CG, has the form diag(6,,J,:6.,J,,: Bal Me EE Bel hn) Be) ience, 
partitioning all matrices in (2.5) in this manner, we see that the off-diagonal 
blocks of Z) must be identically zero. Now consider the diagonal blocks in the 
second equation. The contribution from Z,C, will cancel the contribution from 
C)Z,, and hence we have in this partitioning, if (Z);; represents the 7 diagonal 
block, (Z);;= (Bot AR);; (Z,);;- We can choose the initial conditions to be 
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independent of 2 and determine (Z,);; as a fundamental matrix whose elements 
are entire functions of A. Thus Z, is completely determined, det Z)+-0, aStS0. 
Consider now a off-diagonal blocks of the second equation. (2,);;B,= 
By, (Z1)i; + (Cr); (Zo);;- Since B,,+-B,,, this determines (Z,);; as an entire function of 
A. The diagonal blocks of the third equation give (Z1),;= (Cy),, (Z1);;+ (Co);; (Zo);;> 
and hence (Z,);; may be determined by solving non-homogeneous differential 
systems, and if the initial conditions are independent of /, the elements of (Z,),; 
will be entire functions of A. 

This process may be continued to obtain a formal fundamental matrix with 
the desired properties and of the form Z=UE. 

Since U is a formal power series with nonsingular leading term, Z), Ut is 
a formal power series with nonsingular leading term, Zo. As Z’=aCZ formally, 


' 7’ ZA=C formally, and t U'UA+UC,U4=C formally. Hence, if we 
a 


replace U by a finite number of terms 


1 1 
U=Z,5 741 + 5 2p, 
and define C by 
C=—_UT++ kG 


then C is a convergent power series in - for « sufficiently large and the coefficients 


Gi ay soapy! 
constant ¢,. 


Z=UE solves the differential equation Z’=a«C Z by construction. Writing 
the original differential equation in the form Z’ —«CZ=«(C —C) Z, we obtain 


— tae ‘<a = a. 
the integral equation Z=Z+« {Z()Z14(€) (C—C)Z(é) dé, or U=U+J(U), 
where 


eae) (a, exp {oJ 0, (t) az} T-(2)\ (C—O). 0) [s:sexp{—afo, (x) az})aé. 


J(U) is a matrix which is linear in the elements of U=(U,,), and the coefficient 
of U,, in J/(U) has the form 


fab (t, &, a) exp {a Sot) — o.(r) dr} dé 


k 
where a} are convergent power series in ".. We choose a as the lower limit of 
a 


integration if Re{oe,(t)—@,(z)}<0 and bd as the lower limit of ues if 
Re {o,(t) —@,(t)} 20. Then, since |C —C}<efo? +, ie |<c||U|/«? with an 
appropriate constant cy, 3|U I. 

We can solve U= Ue Fii U) by iteration, wee 0, ieee aor This 
sequence converges uniformly for |A|<A because 


Gu -G,|S2]G—-Guils2~. 
Further, since 
ul 7 ee: US Dea, U;) + (U3 — U,) +: 


. 
> 
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and = ae 
lUlslUl+e+e+---=Ul+1, 


hence s 
7(U)| < c,(|U| = A Vink Cala? 


This inequality implies that the coefficients of 1,...,«~?*? are the same in U 
and U. 

The iteration shows that the elements of U,, U,,... are entire functions of /, 
and the uniform convergence, U,—>U in |A|<A, shows that U is analytic in 
|A|<A. Hence we obtain a solution Z(t, «, 2) = U(t, «, 4) E(¢, «) with the desired 
properties. This completes the proof of Lemma 1. 

The solution constructed in Lemmai1 depends on o. One can construct 
solutions which are independent of o, but we shall not need this result. 


3. Permissible Boundary Forms 
The boundary form s[y,«] is rather arbitrary since multiplication on the 
left by a nonsingular matrix gives rise to an equivalent boundary form which 
gives the same solution to the eigenvalue problem. 
We have made a transformation, y=T7z, of the differential system. This 
transformation will also affect the boundary form s[y,«]. In particular, the 
boundary form for z will be 


(3.1) t[z,a] =s[T 2,0] = P(a) z(a,a«) + O(a) z(0, a), 
where 
Pla) =M(@) Tl) =D Ba, Qe) = Ne) TO) =D yor”. 
i= j= 
Let A and Q, have the partitioning 
BA Pa Pt Qtr Ole Os 
Py=| Pn Pr Pls}, Qo=| Qtr O22 O2s | 
P31 Py, P3s Q31 O32 Q3s 


where the diagonal matrices R,, P,, Ps have orders n, k, m—k respectively, 
and similarly for Qy. 

One would expect that the degenerate boundary form would have associated 
matrices of order ~ whose rows would be linear combinations of the rows of 


0 0 
Pi 11 
By and B1 

0 0 
P31 31 


It turns out that the matrix I’, 
Phe QYs 
(3.2) =| Pe Qis ], 
Br Q3s 
plays an important role in selecting the proper linear combinations of rows for 


the degenerate boundary form. 
Arch. Rational Mech. Anal., Vol. 7 16 
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If / has maximum rank, m, then there exists a nonsingular constant matrix F 
of order n-++-m such that 


pear 
(3.3) hk Pee, ao ; det P9,+0, det 02,+0. 
0 Q3s 


The boundary form tz, a|—Ft[z,«] is an equivalent boundary form for which 


(3.4) P= Pe P PS, Oh Qs O 
1 etal ie 2 ‘OR 2, ae 


With this adjustment, the degenerate boundary form will have Py and Q?y 
as associated matrices, and the degenerate eigenvalue problem is defined by 


ee es Ors Qi 0) 
Qo= 


fae = {B, 11(t) ARGO} a 


6-5) P?, z(a) + 0%, (6) =0. 


The eigenvalue problem (3.5) is reasonable only if the n by 2” matrix (PR, Q,3) 
has rank n. 

Definition. The boundary form S[y, a] of rank n+ m 1s said to be permissible 
if the boundary form t(z, al|—=s[Tz,a] 1s such that (3.4) holds and (P°, Ons has 
rank n. 

4. Eigenvalues 

A(«) is an eigenvalue of the problem (1.1), (1.2) if there exists a non-identically 
vanishing vector y(t, «) which satisfies (1.1), (1.2). If Y(é,«, 4) is a fundamental 
matrix for the differential system (1.1), then any particular solution y (¢, «) must 
be of the form y(t, «)=Y(é, a, A) /(«). Thus, if y(¢, «) is to satisfy the boundary 
condition, we must have 


{M (a) Y(a, a, A) + N(«) Y(b,«, A)}l(@%)=0, andif I(x) =0, 
we must have 


(4.1) A(a, 4) = det {M(«) Y(a, «, 4) + N(«) ¥(b, , a} = 


that is, A=A(«) is determined by (4.1). 


Theorem 1. Let Hypothesis 1 hold and let the boundary form (4.2), s[y, «], 
be equivalent to the boundary form S[y, «] which is permissible. Let the degenerate 
eigenvalue problem be defined by (1.5) and the first n rows of the boundary form 
S[syncol|s 

If Ay ws a simple eigenvalue for the degenerate problem, then, for sufficiently 
large «, there exists in a neighborhood of Ay a uniquely determined simple eigenvalue 
A=i(«) belonging to the complete eigenvalue problem (A.A), (1.2), and A(x) can be 
written in the form 

A (a) = Ay 4 s 


ar 


A(a 
+2), |A@l<e. 


ppt ie 


1 
a 


Singular perturbations of eigenvalue problems PUBys| 


An eigenvector y(t, «) corresponding to A(x) is determined uniquely up to a 
factor depending on « with a suitable choice of this factor y(t, a) can be written 
in the form 


y(t, @) = yo) +21) +--+ 25 vo) + a oa), 


a@< 0,S156,<), (t, x)||<¢,(0,, 62), for every positive integer a, where y(t) 
ws an eigenvector for the degenerate problem corresponding to hy. 


Proof of Theorem 1. Under Hypothesis 1 there exists a transformation T(t) 
which changes the differential system (1.1) into the canonical form (2.2) and the 


boundary form $ [, «] into fi [z, x] for which (3.4) holds. Under this transformation 
the degenerate differential system (1.5) becomes (2.3), and the first 7 rows of 


tz, a] become Pia t) 2 (a, & )+ Ax («) % (0, x). 


We wish to determine the asymptotic expansion of 
A(x, A) = det {M (x) Y(a, «, 2) + N(x) ¥(b, «, a)}. 


If Z(é, «, A) is a fundamental matrix for the differential system (2.2), Y(é, «, A)= 
T(t)Z (t, «, A) is a fundamental matrix for (1.1). Hence A(«, 4) can also be written 


(4.2) A(o,, A) = det {P(«) Z (a, «, 2) + O(a) Z(b, a, A). 


Choose for Z (¢, «, A) the fundamental matrix indicated in (2.4) and let Q be the 
matrix 


(4.3) Q(a, 2) = Pla) Z (a, «, A) + O(a) Z (B, a, a). 


Let the matrix Q have the partitioning (Q;,), 7,7=1, 2, 3, so that Q;; has order 
n, k, m—k for 7=1, 2,3 respectively, and let the other matrices nae similar 
partitioning. 

Since E(a, «)=I, we may write, with this choice of fundamental matrix, 


ie ual 12(@) Z13(a) Es * (0, a) 
Q= 5 P(a) | Ze1(4) Z22(@) Z23(a) Es *(b, 0) | + 
| Bais 32 (4) Zyal(a) Es (b, x) 
Z11(0) Z12(b) Eo (0, a) Z3() Ly 0 0 
+O(@) | Z21(0) Zo2(0) Ea (0,2) Zx3(0) |p {0 I 0 
Z31(0) Z32(b) E (0, a) Z33(0) 0 0 E5(b, x) 


From the nature of the elements in F,(b, «) and E3(b, «) we infer that 
|Z;3 (a, «) E;*(b, «)|>0, 
|Z:2 (0, «) £2 (0, «)| > 0 


exponentially fast as «—>oo. Hence the asymptotic expansion of these matrices 
will be zero, and the asymptotic expansion of A(«, 4)=det Q(«, A) will be the 
16* 


232 W. A. Harris, JR.: 


same as that of 


Z11(4) Z12(4) 0 of Z,1(b) 0 243(0) 
Aa, a) =det !P(a) { Zs1() Zeo(a) 0 )+O(#) | Zr2) 0 Ze0(0) |p det E5(0, a). 
Z31(4) Z32(4) O Z13(b) 0 Z33(0) 


If we let A=D det E;(b, a), D will have an asymptotic expansion 


D(A, «) =Do(A) + Dy(A) +>, 


where oe a 
PS, 24; (a) = Q21 231 (8) 0 0 


DA) = det PY 23; (a) + 21241 (6) Pee Z22(a) t 0) 
Py, Z$1 (a) + Q31 241 (2) 0 233 Z33 (0) 
Since the boundary form is permissible, the determinants of P2, and Q$3 are 
not zero, and since Z$, and Z$, are fundamental matrices for linear differential 


systems, Z$.(a) and Z3;(b) have a non-zero determinant. Thus, Dy (A) =C, 4 (A), 
and A(«, A) can be written 


(4.4) A(u, 4) = C, F(A, «) det Ez (b, «), 
where 
(45) F(a) =A) +R) ++ +E) + er Golo), 


F, (A) =det {P°, 29; (a) + 09, Z21(B)}, ||pol<cy, ¢ independent of 4 and a, a>ap, 
|A|<A. Further, F(A, «), F(A), ..., (A, «) are analytic in A, |A|<A. 

We note that the zeroes of f(A) are the eigenvalues of the degenerate eigen- 
value problem, and the zeroes of F(A, «) are the eigenvalues of the complete 
problem. If A, is a simple eigenvalue for the degenerate problem, then Fo (A) =0, 
OF 
is 
the equation F(A, ~)=0O possesses a uniquely determined solution A(«) which 
satisfies 1 (co) =A, and can be written in the form 


(Aj) ==0, and one can show (e.g. by iteration) that for sufficiently large « 


1a) = Ay el ap ah Malis eS 


ae | ord d 
Cc, independent of «, «> %. 


An eigenvector y(t, a) corresponding to the eigenvalue A(«) is y(t, «) = 
T(t) Z(t, «, A(«)) L(x), where /(«) is a solution of the matrix equation 


(4.6) Q (a, A(a)) L(x) =0. 


We know that A (a, A(«))=det Q(«, A(x))=0, hence one solution is given by 
the cofactors of the elements of any row for which all the cofactors do not vanish. 
The asymptotic expansion of the cofactors of any element in the first rows 
and columns of will have for leading term the corresponding cofactor of 
the same element in (BN Zo ai Olea) multiplied by det (PP, 282 (a), 
det (233 233(b)) and det £3(b, «). All such cofactors could not vanish, for this 


: PF ‘ : : 
would imply that if ‘ (49) =0, “hich contradicts the assumption that A, is a 
simple eigenvalue. ~~ 
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. A calculation similar to that used in evaluating A(«, 4) will show that, after 
division by the factor det E,(b, «), /(«) may be written in the form 


I © 0) 
(4.7) Woes | Of 0 fot tat, 
0 0 E3*(b,«) 
ly 
where Jy>=| 02], Zi (¢,a,A(«)) 2? is an eigenvector for the degenerate eigen- 
13 


value problem (3.5). Thus y(é, «) = T(t) Z (t, a, A( (x)) 2(«) = Pi) Zot LZ} 


on 10 0 Te dOtarty 6 
OnE Emr (a0 OF dijiashd {n+ tat-ot. 
0 0 E3(t, a) O © 1 Re (i a) 


Since | g(t, «) E3(b, «)4|+0, a<t<b, exponentially fast 
as a—>oo, for a<t¢<b6 and uniformly in any closed subinterval a< 6,Sit<6,<b, 


¥) (t, x) = T (2) {% (t) + < 2, (t) + oh, where 
Za (t, 0, A (cx)) ih 
Za\t) = 0) 
O 
Hence 
y (é, x) = Yo(t) + on) ++ ire ai Vo+1 (t,o), a<o,StSd,<), 
llyo+1 ( t, m )| — Cg (01, d») , 
and y,(¢) is an eigenvector for the degenerate problem. This completes the proof 
of Theorem 1. 


Corollary. Jf we assume in addition to the hypothesis of Theorem 1 that the 
elements of A(t), B(t,«), R(t), M, and N are real and that Y(t, «, A) can be chosen 
so that Y(t, aAienY (tc: A) C, where C is a constant matrix with determinant 1, 
then if Ag ts veal, A(x) ts real and y(t, x) can be chosen real. 

Proof of Corollary. Clearly by construction Y(f, «, 4) is a fundamental 
matrix for system (1.1) with A replaced by A. Also Y(t, «, A) satisfies (1.1), and, 
taking complex conjugates of both sides, we see that Y (t,«, A) is also a fundamental 
matrix for (1.1) with A replaced by 2. Hence 


Soir a sui Gioia ,A) 
= det {M (a) Y(a, «, a) + N(«) Y(b, 
= A(a, a). 


Y 
j 
d)}- det C 


Thus, if A, is real, A(«) is determined uniquely by Ay=4y and A(«, 4)=0. Further, 
an examination of the way in which y(¢, «) was constructed shows that y (é, «) 
may be chosen real. This completes the proof of the corollary. 
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5. Self-adjoint Problems 


For a two-point boundary problem of the form 


y=Ahyt+ABy, 
M y(a) +N y(b) =0 
G. A. Biss [/] introduced the concept of self-adjointness under a nonsingular 


transformation z= T() y and considered in detail a special class termed defimitely 
self-adjoint. This concept was modified and extended by W. T, Rew [9-H] 


Consider the differential operator L, 


(5.1) L[y, €] = Az (t, €) y’ + Ao& «) y, 


where the elements of A,(¢, ), A(t, €) are infinitely differentiable with respect 
to t, a<t<b, and analytic with respect to ¢, OXeXeé,). Corresponding to this 
operator is the classical Lagrange adjoint differential operator L*, 


(5.2) L*[z, 2] =(— A¥@ 8) 2)’ + 488) z, 


where A* represents the conjugate transpose of A. Hence, corresponding to the 
boundary problem 


(5.3) L[y,eJ=AB)y, sly]=My,e) +N y(,2)=0 
is the adjoint boundary problem 
(5.4) E*|z,el=A Bt) 2, — tl2) = MM *2 (G2) Lana (bce) OG: 


Following the terminology of REIp [//], the boundary problem (5.3) is said 
to be symmetrizable under a transformation 


(5.5) 22) = Tere) 2) 


if: (i) T(é, e) is infinitely differentiable with respect to t, a<t<d, analytic with 
respect to ¢, OS eSéq, and nonsingular, a<t<b, 0<eXe@Q; 

(u) for arbitrary values of A and ¢, O<e<é , a vector y satisfies the diffe- 
rential system or boundary condition of (5.3) if and only if the corresponding 
vector z of (5.5) satisfies respectively the differential system or boundary con- 
dition of (5.4), in which case we say that (5.3) and (5.4) are equivalent; 

(iil) the associated matrix S(, e)= T*(t, e)B(t) is hermitian, aX<t<b, 
OS eS. 


An important instance of symmetrizable problems is the case of problems 
that are self-adjoint in the classical Lagrange sense; that is, 


Lly,e]—ABy =L*[y, e] —AB*y, 


while s[y]=0 if and only if ¢/y]=0. In terms of the coefficients, the conditions 
are d;= — Af#, B= B*, Aj =A,—A¥, M* Aj *(a, 2) M=N* AP? (0, e) N. 
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Rei [//] has shown that if (5.3) is equivalent to its adjoint under the trans- 
formation (5.5), then there exists a suitable transformation T, (¢, €) such that the 
boundary problem 


(5.6) T* (t, 2) L[y,e)=AT* By,  s[y]=0 


is self-adjoint (in the classical sense). However, if the coefficient matrices Ay, 
A,, B are real-valued and T is real-valued, JT, may not always be chosen real- 
valued 


Let 
<y (tj, 2()> = f *@ yO dt. 


Definition. Let & denote the linear space of vector-valued functions satisfying 
s[y]=0, L[y, e]=Bg for all g(t, e) infinitely differentiable with respect to t and 
analytic with respect to ¢, O< eX. 

For arbitrary real constants c,, c, (or real-valued functions of the parameter ¢) 
the formal operator 


Ly, 23 05 65 | = Pe lerl ly, 2-6 c By) 
is hermitian on Y in the sense that 


(5.7) Ll, €3 G, Co; T], Yoo = (1, L[ ve, €3 1, Ce; T]>, 


V1, ¥2 members of #. We note that if 7* A, is skew-hermitian, then (5.7) will 
also hold for y, and y, such, that s[y,]=0O=s[y.|. In either case, <L/y, ¢; 
Gi» Cosel |,.9)> 1s neal \valued. 

The problem (5.3) is called definite [c,, cy; T] whenever (5.3) is symmetrizable 
under (5.5), and for this J and suitable constants c,, c, the functional 


(5.8) ILy, €; G, €,; T] = <L[y, €; G, Cy; t1,y> 


is positive for all y in # with By =0, aXi<b. 

The problem (5.3) is called normal if there are no solutions of L[y, ¢]=0, 
s[y]=0, y =0, such that By=0. 

A problem that is normal and definite [0, 1; 7] may be treated by methods 
similar to those of Briss [7]. A problem definite [c,, c,; 7] with c,=-0 is equi- 
valent to a problem definite [1,0; ZT] through replacing A by A—c,/c, in (5.3) 
and a problem which is normal and definite [1,0; 7] may be handled by the 
method of Re1p [9—J1]. In all cases the eigenvalues are real. 

It has been remarked by REID that if the eigenvalues of (5.3) are bounded 
below then by replacing 4 by A+ A, such a problem can be reduced to one for 
which 


DVeN Os a Level, vs > 0, 


(5.9) T9004 a Nan S. yyy > 2.0 


for yin Y, By +0, S=T*B. Conversely, if we assume (5.9), then it can be 
proved that the eigenvalues are bounded below. 
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Consider a simple self-adjoint case, T=J, 
{AS (é) y' + A(t) y} + {A 9’ + Aoy} =A BY y 
L{y, e] =Lol[y] +eLily]=A4By. 


To insure semi-boundedness, we could ask that both Ly and L, be semi- 
bounded, that is, 


(5.10) Ly], 24 <By. >, Oly, w2a<By, y, 
for yin Y, do, d, finite. However, if we set A=A+ d,+ ed, Lo [y] +d, By=Lyy, 
and fe [y]+d,By=L,[y], we get an equivalent problem L,[y]+eL,[y]=ABy 
for which <Ly[y], y>20 ue <L{y, y>20. 

Similarly, if L[y, e]= Ds L,[y] &, <L;[y], y>=4;<By, y>, where d(e)= Pa d; é 
is convergent forOSeS ae then replacing A by A+d(e), L,{y] by B ily]- ees "B y, 
L[y, €] by Lly, é| +4 (e ) By, we get an equivalent problem Lily, é| =AB Ms for 
which formally <L[y, e], y>20, yin &. Further, if L’ is the operator a 


or 


aH aA , 
L Ly, é| = au y 
then formally Om [y, €], Yy>=0, y in & 


We shall, however, only make assumptions on the operators L and L!, where 
L{y, €] = Ay(t, €) v + Ag’, €) y, 


oA oA 
(5.11) L'[y, 2] =>4 be) se he) y- 


Theorem 2. Let the conditions of Theorem 1 apply to (1.1), (1.2). Let (1.1), 
(1.2) be such that 1t may be written in the form (5.3) with A=1/e for which the equi- 
valent permissible boundary form is independent of e. Let the eigenvalues of (1.1), 
(1.2) be analytic functions of &¢, O< eo, and let the degenerate eigenvalue problem 
have only simple ergenvalues. Let (5.3) be normal and definite [1, 0; I] and normal 
and definite [0,1; I]. Let the operator L' as given in (5.11) be such that <L[y, €], y»=0 
for yin £. 

Then lim A, (€) =A, (0)<00 for every n=1,2,.... The values A,,(0) are eigen- 


values of the degenerate problem. For small values of ¢ the behavior of all eigen- 
values j,,(€) 1s given by an asymptotic series with respect to powers of e. 

Remark. A problem that is normal and definite [0,1; 7], and normal and 
definite [1,0; 7], can always be adjusted to be self-adjoint, that is, normal 
and definite [0, 1; 7] as well as normal and definite fay OF ie 


Y mln — 1), 2. (0 EAL ee _ 


Proof of Theorem 2. If A(«) is any eigenvalue that is analytic for 0<e< Eq; 
and lim A (2) =p, then A(e) is a solution of A(e,2)=0 and 4(0)=A,; hence 2, 
is a solution of the degenerate eigenvalue problem, and A(e) is one of the eigen- 
values that we have constructed. Thus the proof of Theorem 2 is reduced to 
the proof of 


(5.12) lim A, (e) = Ao 
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Let y(t, e) be an eigenvector corresponding to the eigenvalue A(z) (which 
is real) and analytic in e, 0<eSéy, and normalized so that <By, y»=1. Thus, 


(5.13) A(e) = <L[y, €], y> =I[y; 1,0; I]>0. 
Further, 
d y 
= (12, e], y+ (Ly, 2], 2) + Ly, el. 9». 
Since 
nat, (Bs)+(BxB)=0, sp1=o=[2) 
and 
Lge 4-9) = Goer El el) = Ge Mle) By) = M0) GE BY). 
Thus, 
e1 Ben) + (eo.e Bt) —200 (2.89) +90 Be} 0 
and dA/de becomes 
(5.14) oh = <L'[y, €],. y= 0. 


Hence lim A(e) exists since A(e) is non-negative and a non-increasing function 
e—>0 


of ¢ as « decreases. This completes the proof of Theorem 2. 


6. Examples 


To illustrate the results of this paper consider the following eigenvalue problem 
considered by Lord RAyLeicH [8], which is a special case of a problem con- 
sidered by J. Moser [7]: 


sae ea CSta05 


x (0) = x'(0) = x(1) = x’(1) =0, ~-e-small and positive. 


Let y be the vector with components y;, where y,=¥%, Yg=—%’+ 67x", 
Vg=Ex"', Yy=erx'". y will satisfy the differential system 


== jl 


Ci = + ay 


Cm. 0) 04 0-0, 0.0 
000 0 0 Wes ertey x0) 
€ 001 00 0000} 
01.0 00 1 @, 0-0 


Qa @ Se 
qr OF. SO 
= eS 


and the boundary conditions for y corresponding to those for x will be 


(6.2) s[y,e]=My(0, 0) + Ny (1, 6) =0, 
where 
45 0) 0) (0) ORO ROCs () 
Rp 4 Oe yt I0.,0 0 0 
OR Om One i @ @ 
Om OMORO) 010-1 
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We see that Ayo= (2 al which has +1 as characteristic roots, and the 


transformation y= Tz, where 


AO Oy 
Oat OF x0 
hax Oe S183 
00-11 
will change (6.1) and (6.2) into 
OQ © O @ Opt OPO OsOMG 
et) Ag Om OL LOY 80 Oo O OmO ] AO FORO - 
CP anion emesis Arr oe Wma —1)0 0 OS 
OO... Ow at ONG 00) + 000 
 OY.O OR OR Oras Ome 
OF OF. O76 
= O,ée) + EGU) et) 
Cad T= | Ging ovo 1141 Sra hen en knee 
OOOO One tee Ty, 
Thus 
0.40 
pau 1 2m) 
OF70 
QO —1, 
which has maximum possible rank, and we may choose 
DA 2 ORO 
Pes OM Om ten ; 
Or Oy ae 
OOO Rat) 
so that 
Or OQ C Ona Om.0 
i OROMOREO ee) O O 
OS) “ls el" ae el = 2(0, 1,€ 
(65) iz 6]=Felze]=(P 9° ° lenet( fh 0° P lette), 
Om a © Ooh, Saal 
O O 
he 0 O 
1t 10) 
Old 


. vy ‘ ; 
Since (4 i ; 3 has rank two, Theorem 1 applies to the eigenvalue problem 


(6.4), (6.2) witha= 1. 
coy 
The problem (6.1), (6.2) is symmetrizable under the transformation U, 
a1 0 0) 


Oo —eée O 


@) 
6.5 pee 
oe Ome QO —e]?’ 
0) 
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where U*=—U, MU141M*=NU14N*=0. Hence (6.1), (6.2) can be written 
as the self-adjoint problem 


Ome O40 (OP OO) T0020 
Seely See lees Oey ea leae) OOo 6 

6.6 : —=J 

9) SERCO Keto Vi cams Od ‘ia net Od ed 
OME P7210 Oe OrniOren 4 0000 

(6.2) M y(0, 2) +Ny(1,0) =0. 


/ 
In this case, 
EAS OS 13 =fx{e M9 — AM dt af aae, 


L 
and if f x?=1, then an integration by parts shows that 
0 


1 1 
A(e) =e? f (x22 dts f (1)? a7; 
0 0 
hence dh 


The more general eigenvalue problem of the form (P+ ¢?"-?"Q) x=Ax, 
where P and Q are linear differential operators of order 2m and 2m respectively, 
n>m, with n boundary conditions prescribed at each end point, t=a, t=), 
has been treated by J. MosER [7]. 


2m 2n 


k 
Px= Yan) x), Ox=Yin 2, w= FR, astsd, 
k=1 


p(t), 9, (t) are infinitely differentiable real functions, with boundary conditions 
A, eR gon (a) a >; ays a) (a) — 0 
Sor =a Wee 


B, % = x‘) (b) + >) b,, x!) (b) =0 


$< By 
where a,, and 0,, are constants. 
If we let y be the vector with components y,;, where 


y, = xt), VA wee 2H 
Vom ane Ns toay oe GanlPam Veons 
Yam+j = pik a 7= 15 2n — 2m, 


then y will satisfy a differential system of the type ae with «=1/e where B (i, «) 
is a polynomial in e of degree at most 2m —2m and A(t) has the form 


[a0 
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The characteristic roots of Ay, are the roots of the equation ¢,,, B?”-°" + pom=O0, 
and if (—1)"¢.,>>0, (—1)"Pem-> 9, the characteristic roots will have non-zero 
real parts. The boundary conditions A, and B, may be written in the form (1.2). 
In particular, if %,,<2m, By<2m, A, cae), Boxma), r=m+1,.., , then 
the boundary form may be chosen independent of «. 


m n 


If P and Q are formally self-adjoint, Pa=> (p, x), Ox=> (G, 2"), 


k=0 n k=0 
then P+e2"-2"Q is formally self-adjoint, (P+e?"-?"Q) x=) (r,x). Let 
k=0 


g be the vector with components g;, where Tae eae 1 ue, } Onti= 
(1) (7,09 4 (ty 84 EO, 741, -, n. The vector g will 
satisfy a first order differential system that is symmetrizable under the trans- 


formation ] = 0 ‘| (see [11], p. 419; [3], p. 206). Further, p=V(é, «) y; 


—I0 
hence the differential system for y is symmetrizable under the transformation 
U(t, e)s=V*it, e) J Ve, e)=— UF. 

The differential system for y may be written in self-adjoint form by multi- 
plication on the left by U*, and in the resulting self-adjoint form the matrix 
coefficient of A has 1 in the first row and column and 0 elsewhere. Thus 


b b b 
TlyyoryT) ACY edt = fe Pea ee Oe PO eae: 


Further, if (—1)",,>0, (—1)"¢2,>>0, the problem could have been adjusted 
b b 


so that fu Pudt=0, fuQudt=0 for all wu satisfying the boundary conditions 


ALu==0, Bye=0; Tel, a, fs Elemce, Ale) 22%, —_ E210) 
é 
This paper was sponsored by Mathematics Research Center, U.S.Army, Madison, 
Wisconsin, under Contract No. DA-11-022-ORD-2059. 
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Uber Instabilitatshereiche eines numerischen Verfahrens 
zur Losung des Cauchy-Problems 
fir hyperbolische Differentialgleichungen 


R. ANSORGE und W. TORNIG 


Vorgelegt von L. COLLATZ 


1. Einleitung 
Die folgenden Betrachtungen beziehen sich auf die von einem der Verfasser 
mitgeteilten Verfahren zur numerischen Lésung des Cauchy-Problems bei par- 
tiellen hyperbolischen Differentialgleichungen 2. Ordnung [2]. 
Vorgelegt sei das Cauchy-Problem 


Bax y = f(x, y, 4; Ses ahs 


2(%,.4) = O(X) 2. 2%) OA eal Bae Oh 


(1.1) 


(In [2] wird gezeigt, wie man das allgemeine Problem auf (1.1) zuriickfiihren 
kann.) Zur Lésung werden die in [2] angegebenen Verfahren benutzt. Bei 
exakter Rechnung (die infolge unvermeidbarer Rundungen im allgemeinen nicht 
durchfiihrbar ist) liefere das Verfahren bei Berechnung der Werte oberhalb y= x 
auf den zu der Anfangsgeraden g): y= parallelen Geraden g,: y=x+nh die 
Lésungen 


Zz LN aN (AO Nien oc (es Se at a Re ener, 28 (1.2) 


1,4+n? 


Diese unterscheiden sich von den exakten Lésungen des Problems (1.1) noch 
durch den hier nicht betrachteten Verfahrensfehler. Benutzt man nun auf den 
Geraden g,, fiir die weitere Rechnung statt (1.2) die durch Rundungsfehler etwas 
verfalschten Werte 


Sy y+n =I Ny, r+n? Bx», r+n at Ugcencn ’ Ayn rn =f Nyy pen? (1 3) 


so liefert das Verfahren auf den Geraden g, (@«=N+1, N+2, ...) unter der 
Voraussetzung, daB exakt weitergerechnet wird, veranderte Werte 


An7 +6 aie "rrt+o> 2 xp, r+o oh "er, ro? Buy, r+o ae Ny, rro° (4 4) 


Wir nennen das Verfahren stabil, wenn die Fehlervektoren 


/ 
br,r-+o ar (jp s-0 Necr, r+o Nyy, ae) (1.5) 


fiir alle y mit wachsendem o bei beliebiger Anfangsfehlerverteilung beschrankt 
bleiben. Andernfalls hei&e das Verfahren instabil. 
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2. Voraussetzungen 


Es st6Bt auf groBe Schwierigkeiten, die Frage nach dem Instabilitatsbereich 
in voller Allgemeinheit beantworten zu wollen. Aus diesem Grunde wird in der 
Literatur bei ahnlichen Problemen im allgemeinen davon ausgegangen, daB das 
Differentialgleichungsproblem linear ist und die Koeffizienten im betrachteten 
Bereich als konstant angesehen werden kénnen. Weiterhin wird in der Regel 
von einer bestimmten Anfangsfehlerverteilung ausgegangen. Auch in dieser 
Arbeit wird von diesen einschrankenden Voraussetzungen nicht abgegangen. 
Es gelingt dann, Bereiche anzugeben, in denen das Verfahren fiir die gewahlte 
Anfangsfehlerverteilung instabil ist und daher auch allgemein als instabil be- 
zeichnet werden muB. 


Das Problem (1.1) laute also 
Z,4y=2@(%, vy) +Az+ Bz,4+Cz,, 
a Me A) — Gta), 2A Xa ed (5), By (%5 a) =o (x) 


(2.1) 


mit konstanten A, b, C. Auf einer ersten Geraden g,, seien alle Fehlervektoren 
gleich groB, d.h. 3,,,, ist unabhangig von 7. Bis zur Geraden g,_ sei exakt 
gerechnet worden. Diese Voraussetzung ist nicht einschneidender als die, dab 
exakt weitergerechnet wird. 


3. Instabilitat des Verfahrens 


Es geniigt, die Instabilitat des Interpolationsverfahrens fiir die fortlaufende 
Rechnung zu untersuchen, sofern nicht ausschlieBlich mit dem Extrapolations- 
verfahren gerechnet wird. Wir beschranken uns dabei auf die Formeln zur 
Berechnung von Werten oberhalb y=x, weil man bei der Untersuchung der 
Werte unterhalb y= ~ letztlich nur # durch —/h zu ersetzen hat. 


Das Interpolationsverfahren m-ter Ordnung zur Berechnung der Werte ober- 
halb y=~ lautet: 


m mm m— 
totaal 2|,*m™ s[*+1] vy BA 
ey ( Y, Rais h 2, Cost nese sty h oa r,st1 sre = On v alte st+li—p =0 


u=0 size v=0 
; L, v0, 0 
mm 
Pret — {Pr thot oft + > nee ll =i (3.1) 
y=1 


m 
1 rik = 
ee — aes —h mote cia 7 2. en ee eal 0), 


u=1 


Dabei sind die iiblichen Abkiirzungen /=2,, g=z, benutzt worden. Die Werte 
der auftretenden Gewichte entnehme man [2]. Nach (2.1) ist hier zu setzen 


jf (Gey cs Padre (*.V)- Aee BD 4 Cq. 


Die Schrittweite h sei so gewahlt, daB das Iterationsverfahren (3.1) konvergiert 
(s. [2]). Die Iterationen seien jeweils bis zum Stillstand durchgefiihrt, so dak 
die Indizes [x] fortgelassen werden kénnen. Dann kann (3.1) als implizites Glei- 
chungssystem fiir die Werte z,.:1, Py si1» Gs¢1 aufgefaBt werden. Es liefert 
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ay s+1 2p s41 1 Yr, st 
Dy sn verfalschte Vektoren (apa ae ak) abet ney 
Yr,s+1 r,s+1 5 § "yn, s+1 


sofern die in (3.1) auftretenden Werte langs der verschiedenen vorangehenden 
Geraden g, (o<s+1—7) mit Fehlern behaftet sind. Wie man leicht einsieht, 
errechnet sich dann der Vektor 


, 
ons = (Geass Nar, +1 Tica) 


aus einer Matrizengleichung der Form 


m m—[h 


ie br,s+1 = 2 > ie br+mu,st1—» fiir m = 1 , 
pu=0 v=0 
0,0 
4, v9, (3.2) 


tol Maya ON by, + eb 4a, 541 fiir m= 0. 


Ist g, die erste Gerade, auf der die Werte mit Rundungsfehlern behaftet sind, 
und sind die Fehlervektoren 3, ,,,, langs g, alle gleich, also unabhangig von 7, 
so sind nach (3.2) auch alle Fehlervektoren 3, ,,,,,, langs g,+, (k=1, 2, ..-) gleich. 
Wir schreiben deshalb 3, statt 3,,,, und erhalten mit (3.2) 
m—1 
va eae pS i ba» fiir Mm = a 
v=0 
110] 4a NBS hs fir m=0. pe 
Fiir die weiteren Betrachtungen ist es im Falle m= 2 zweckmabig, die folgenden 
Fehlervektoren u!”! einzufiihren: 
Ns 
Nas 
"yy, 
Ns—1 


m Nee 
ut ee a e (3.4) 


Ns—m-+1 
Uke d 


Ny .—m+s 


Dann ist (3.3) aquivalent einer Gleichung 
ato) ll, — tml yl 6.5) 


Fiir die Determinante der Matrix @!”) gilt: |@”™| =|]. Sie sei von Null 
verschieden. Dann ist 

al, gylm™® Lom ym) gato) gym 
und daher 


rs a (3.6) 
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GemaB8 unserer Definition ist das Verfahren somit héchstens dann stabil, wenn 
die Folge der Matrizenpotenzen {8"""! (k=1, 2, 3,...) beschrankt bleibt. Im 
anderen Fall ist das Verfahren sicher instabil. Fiir die Beschranktheit der Folge 
{pen} ist es notwendig, daB kein Eigenwert der Matrix 8"! dem Betrage nach 
groBer als 4 ist. Sind alle Eigenwerte dem Betrage nach kleiner als 1, so konver- 
giert die Folge {8'""} gegen die Nullmatrix. Gilt |A;|=1 fiir nur einen Eigen- 
wert, wahrend alle anderen Eigenwerte dem Betrage nach kleiner als 1 sind, 
so bleibt die Matrizenfolge beschrankt. Sind mehrere Eigenwerte dem Betrage 
nach gleich 1, so kann die Folge sowohl beschrankt als auch unbeschrankt sein. 


4. Der Fall m=0 
Die in (3.3) auftretenden Matrizen ! und $f"! ergeben sich zu 


h? h? h? 
Te ee {toe 
Chir Si awe PB — pC te oe 0st 0 
ha (pBae. LRG OeOhsd 


Zu untersuchen sind die Eigenwerte von "BY. Diese Eigenwerte sind 


1—h? A 1 1—-hA\? WA 
pay Wee 2] at°)] = Diois (4.1) 


Dabei ist 


[WO] =144C—2B +" A. 


Abkiirzend setzen wir h?A=a, hC—hB=f und ermitteln den Instabilitats- 
bereich in der «-f-Ebene, da nach (4.1) die Stabilitat im wesentlichen von 
diesen beiden GréBen abhangt (im Falle «=0 mu die GrdBe von B noch in 
Betracht gezogen werden). 

In genau folgenden Bereichen sind A, und A, dem Betrage nach kleiner als 14 
und damit die Folge {[9°" BO" (k=1, 2,...) beschrankt (s. Abb. 1): 


: <<a) <a—— 1) 
; i P (4.2) 
2, O<a<4, Dee ed, Pees Oi. 
Auf den Randern 
Ch< (Op — — 
a=0, eto, 
% ==10; 2303 (4.3) 
<< Se ees 
Oa oS 4e 5° 
a4, = — 2 


dieser Gebiete ist mindestens einer der beiden Eigenwerte A, und A, dem Betrage 
nach gleich 1, wahrend nur in der restlichen «-6-Ebene mindestens ein Eigenwert 
dem Betrage nach gréBer als 1 ist. Es bleibt also noch das Verhalten der Folge 
auf den Randern (4.3) zu untersuchen. 
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Wie man mit Hilfe der Jordanschen Normalform einer Matrix erkennt, kann 
im Falle | A,| eagaipesd 5.2) die Beschranktheit unserer Matrizenfolge auf dem 
Rande héchstens dann gestort sein, wenn mindestens zwei Eigenwerte einander 
gleich und dem Betrage nach gleich 1 sind. Das ist nur méglich auf der Geraden 
a—=0 (dann ist 4, =A,=1) und im Punkte «=4, B=— 2 (dann ist A, =/;=—1). 
In allen anderen Randpunkten ist unsere Matrizenfolge daher beschrankt. Da 


der Fall «=4, B=—2 aus weiter unten genannten Griinden nicht interessiert, 
FEAL vo: geben wir noch das Ergebnis 

~ (wey <g im Falle «=0. Hier besteht 
Begrenzung des moglichen fe Beschranktheit der Folge 
ok (a, Bereiohs /m Fallem=0 2 fpactor™ Bly nur, wenn b=0 
ys Lrmittelter Instabitatsbererch cedaleR B>0 ote B <_ 216m 


des Verfahrens m=O 


Auf Grund der Konver- 
genzbedingung fiir die Itera- 
tion der fortlaufenden Rech- 
nung kénnen nun aber die 
Werte « und f a priori ge- 
wisse Schranken nicht tiber- 
schreiten. Der Durchschnitt 
des dadurch  definierten 
Wertebereichs mit der Kom- 

Abb. 1 plementarmenge des oben an- 

gegebenen Beschranktheits- 

bereichs der Matrizenfolge liefert somit den Instabilitatsbereich des Verfahrens 
fiir die eingangs beschriebene Anfangsfehlerverteilung. 


Die Konvergenzbedingung lautet 


Te 
n<2(1+s% -1), (4.4) 
wobei M die Lipschitz-Konstante fiir die Bedingung 


WAG? 20) Pada) — 1 (Xs Ms 2,21, %)| < M(|z— Z| = |p2— p,| =e | qo — ai|) 


darstellt. Man wird M = Max {|A|, |B, |C |} wahlen und mit diesem M die Schritt- 
weite h aus (4.4) bestimmen. Dann folgt 


jal sie <2—sM (1+ 54,1) <2, 


2M 
ae 40) 
ue 2am <aM (1+ 54, —4] Step 

Somit ergibt sich folgender Instabilitatsbereich (s. Abb=4): 

a<0O, 6 beliebig; 

a=0, B=+0, _ £8 beliebig; 

¢=05) pBiss Ooh «PS 0; 

on 


0, a 
o% po a 


Im Falle der Berechnung der Werte unterhalb y= x ist h durch —h zu ersetzen. 
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5. Der Fall m=1 
Die in (3.3) auftretenden Matrizen 2!) und Bf! lauten: 


2 2 2 
eG fee A Me hone G 
| 6 3 3 5 
gi) — Metiae| Pema eulRee , S2— A fe BR aes 
2 2 2 2 

h h h h h h 

lige guateie BSB oy oth toed be 

2 2 a 2 “ Z) are) Z . 
Die zu untersuchenden Eigenwerte von 77 B14 sind: 

2 eS Ee. 2 
— A | 4— : A 0) 
A=1; dgg=— Rol = E 1 ja()] jue] (5.1) 
mit 
2 
| | = 14 te eae Lae 


Rad Seschrankthertsbereich von 
{a a B 
____ Begrenzung des moglichen 
(«,f)-Bereichs im Falle m=1 
yy, Ermittelter Instabibiatsbe ~ 
reich des Verfahrens m= 1 


LQQQ erovm 
Wd , 


Abb. 2 


Die Betrachtungen zur Ermittlung des Instabilitatsbereichs verlaufen zu denen 
in 4. analog, so daB wir uns hier kiirzer fassen k6nnen. Die Konvergenzbedingung 
fiir die Iteration der fortlaufenden Rechnung lautet jetzt 


n<s(\ir 2-1) 


Mit M=Max{|A], | B|,|C|} und den gleichen Abkiirzungen wie in 4. folgt 


|| = 6, |B | SD; 
Innerhalb dieses «-f-Bereichs besteht Instabilitat, wenn (s. Abb. 2) 
C= Owe p, belichig; 
—0N =O.ep DelleDigs 
Cia Oho 3) 2105 et (37105 
Dee ep. 


Im Falle der Berechnung der Werte unterhalb y= x ist / durch — h zu ersetzen. 


yes 
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6. Der Fall m=2 


Wir betrachten noch den Fall m2 und bilden zunachst ahnlich wie in 4. 
und 5. die Matrizen 2"), B?! und BI]?! und mit ihnen die Matrix 


[232 92] [21-2 9a i2) 
spl] — U Bo » Bi ; 
G © 


Zu bestimmen sind die 6 Eigenwerte dieser Matrix. Man errechnet 4,=/,=0, 
Az=1. Die restlichen 3 Eigenwerte sind die Wurzeln der Gleichung 


1+73a+3 6 , 2a+B yy ate 
(A—41)3 + ay (A —1)24 at) (A —1) + joel] —3() (6.1) 
mit 


|2]| =1+ yeB + eae. 


Dabei wurden die gleichen Abkiirzungen wie friiher benutzt. Das Verfahren ist 
sicher dann instabil, wenn eine der Lésungen von (6.1) dem Betrage nach gréBer 
als 1 ist, was sich fiir jede spezielle Differentialgleichung leicht nachpriifen laBt. 
Im Falle mehrfacher Eigenwerte vom Betrag 1, z.B. fiir «.=0, kann wiederum 
die Transformation der Matrix $S!@] auf die Jordansche Normalform Auskunft 
iiber die Beschranktheit der Folge {$2!"} (k=1, 2, ...) geben. 


Die Verfasser danken Herrn Dr.-Ing. R. ZURMUHL fiir einen matrizen-theoretischen 
Hinweis. 
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Feblerabschatzungen bei hyperbolischen 
Differentialgleichungen 


WOLFGANG WALTER 


Vorgelegt von H. COLLATZ 


1. Einleitung 


Das Hauptziel der vorliegenden Arbeit besteht darin, fiir die hyperbolische 
Differentialgleichung 


(1) Uny=](%, ¥,U,U,,u,)  [bzw. =f(x, y, 4)] 


Abschatzungssatze aufzustellen. Daneben werden auch Systeme von solchen 
Differentialgleichungen und einige verwandte Probleme betrachtet. Von den ver- 
schiedenen Aufgaben fiir die Gleichung (1) behandeln wir die charakteristische 
und die Cauchysche Anfangswertaufgabe; genauer gesagt legen wir eine von 
KisyNskI [5] angegebene allgemeinere Formulierung zugrunde. Sie bietet den 
Vorteil der gleichzeitigen Behandlung von charakteristischem und Cauchyschem 
Problem. Mit Existenzfragen beschaftigen wir uns hier nicht; Existenzsatze unter 
sehr schwachen Voraussetzungen tiber f/ finden sich in [8], [9]. 

In Arbeiten von Diaz [6], TOrNiIG [70] und Moore [//] wurden Verfahren 
zur numerischen Behandlung der Gleichung (1) untersucht. Es handelt sich 
dabei um die Ubertragung von Verfahren, die bei gewdhnlichen Differential- 
gleichungen bekannt sind, etwa des Euler-Cauchy-Polygonzugverfahrens, der 
Extrapolations- und Interpolationsverfahren, auf hyperbolische Differential- 
gleichungen (in [/7] wird eine Ubertragung des Runge-Kutta-Verfahrens ange- 
kiindigt). Hier soll nicht ein weiteres Verfahren zur naherungsweisen Berechnung 
einer Losung besprochen werden. Vielmehr werden allgemeine Satze bewilesen, 
welche gestatten, von einer vorliegenden Naherungslésung auf das Verhalten 
der Lésung zu schlieBen. Es werden dazu zwei verschiedene Arten von Satzen 
angegeben. Bei den Satzen der ersten Art — sie werden mit A bezeichnet — 
handelt es sich im wesentlichen darum, die bei gewohnlichen Differentialglei- 
chungen bekannten Begriffe Unterfunktion, Oberfunktion auf zweidimensionale 
Probleme auszudehnen. So besteht z.B. der folgende Sachverhalt: Liegen bei 
einem charakteristischen Anfangswertproblem die Anfangswerte der Funktion 
v(x, y) unterhalb, die Anfangswerte von w(x, y) oberhalb der fiir die Lésung 
vorgeschriebenen Anfangswerte und ist v,,</(%, ¥, v), Wyy2f(%, y, w), so liegt 
die Lésung zwischen v(x, y) und w(x, y). Wir nennen dann v eine Unterfunktion, 
w eine Oberfunktion. Aussagen dieser Art sind jedoch nur mdglich, wenn die 
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Funktion f gewisse Monotonieeigenschaften besitzt. Bei den Satzen der zweiten 
Art — wir bezeichnen sie mit B — werden Abschatzungen fiir die Differenz 
,Losung minus Naherungslésung” angegeben, bei welchen der , Defekt“ der 
Naherungslésung, das ist die GroBe v,,—f(%, y, v, -.-) in die Abschatzungsformel 
eingeht. Sie sind fiir beliebige rechte Seiten anwendbar. Der Spezialfall der 
Lipschitz-Abschatzung von f wird naher dargelegt. Er fiihrt zu einfachen Fehler- 
abschatzungen, die sich numerisch leicht auswerten lassen. 

Mit Hilfe dieser Siatze ist es méglich, sehr allgemeine Eindeutigkeitssatze 
fiir die einzelnen behandelten Probleme aufzustellen. So kann man ohne Miihe 
jene Sitze, welche Sao [4] fiir nichtlineare Volterra-Integralgleichungen in einer 
unabhangigen Variablen erzielt hat, auf Integralgleichungen in zwei unabhangigen 
Variablen, wie sie in Nr.4 betrachtet werden, iibertragen. Fiir die hyperbolische 
Differentialgleichung (1) wurde in den letzten Jahren das Eindeutigkeitsproblem 
von verschiedenen Autoren behandelt. Dabei gelang es, eine Reihe neuer Ein- 
deutigkeitskriterien aufzustellen. Alle diese Kriterien (und noch allgemeinere) 
sind aus den Satzen der Nr.5 und 8 ableitbar. Da wir in einer vor kurzem er- 
schienen Arbeit [12] auf das Eindeutigkeitsproblem ausfiihrlich eingegangen sind, 
verzichten wir hier auf eine nahere Diskussion. In [72] findet man auch Literatur- 
angaben. 

Was die Methode angeht, so stellen wir in Nr.3 zwei allgemeine Abschatzungs- 
satze fiir Operatoren an die Spitze. Sie werden in Nr.4 auf Integralgleichungen, 
in Nr.5 auf Differentialgleichungen u,,—/(x, y, «) spezialisiert. In Nr.7 werden 
die beiden allgemeinen Satze auf Vektorfunktionen ausgedehnt. Damit lassen 
sich dann Systeme von Integralgleichungen sowie Differentialgleichungen der 
Gestalt u,,=/(%, Vy, u, U,, u,) (und Systeme solcher Differentialgleichungen) be- 
handeln. 

Weitreichende Analogien zum Anfangswertproblem bei gewéhnlichen Diffe- 
rentialgleichungen, wie sie vom Existenz- und Eindeutigkeitsproblem her bekannt 
sind, treten auch hier sehr deutlich zu Tage. An gegebener Stelle wird darauf 
hingewiesen. 

Herrn Doz. Dr. K. NicKEL bin ich fiir manche Anregung zu Dank verpflichtet. 
Monotone Operatoren von der hier betrachteten Art wurden im eindimensionalen 
Fall erstmals von ihm eingefiihrt. Er wird iiber seine Ergebnisse in einer dem- 
nachst erscheinenden Arbeit berichten. 


2. Bezeichnungen und Problemstellung 


Unsere Betrachtungen beziehen sich auf einen in der (x, y)-Ebene gelegenen 
abgeschlossenen Grundbereich R. Dieser wird erhalten, indem man von einem 
Rechteck 0S x<a,0S yb, (a>0, b>0) durch einen Schnitt langs einer stetigen, 
vom Punkt (0, d9) (0S by 6) zum Punkt (a9, 0) (0S a) a) verlaufenden, monoton 
fallenden Kurve C, ein den Nullpunkt enthaltendes, krummliniges Dreieck ab- 
trennt (s. Abb. 1). Die beiden Extremfalle sind besonders wichtig im Hinblick auf 
hyperbolische Differentialgleichungen. Der erste hegt, dann vor, wenn a,=b,=0 
ist. Es wird dann gar nichts abgeschnitten, d.h. R ist gleich dem ganzen Rechteck. 
Mit einem solchen rechteckigen Grundbereich hat man es beim charakteristischen 
Anfangswertproblem fiir die Gleichung (1) zu tun. Beim zweiten Fall ist ao= a; 
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by=6. Solche Grundbereiche, die selbst die Form eines krummlinigen Dreiecks 
haben, treten beim Cauchy-Problem (oder nichtcharakteristischen Anfangswert- 
problem) auf. 


Die Kurve C, sei durch eine Funktion y=¥(x) definiert, welche im Intervall 
OS ~<a, stetig und im engeren Sinne monoton fallend ist; sie besitzt demnach 
eine in 0<y<0,y erklarte Umkehrfunktion x= %(y), welche ebenfalls die beiden 
genannten Eigenschaften hat. Wir erweitern nun den Definitionsbereich der 
Funktionen y (x) und x(y), indem wir y (x) =0 fiir a, x<a, X(y)=0firh<y<b 
setzen. Der Grundbereich R kann dann durch jeden der beiden Ausdriicke 


R={(x,y)|OSxSa, ¥(xs)SySd} 


={(x, y)|OSySb, Xy)SxSa} 


beschrieben werden. Unter 7(% 9, V9), (ely)y) 


(%9, Yo) € RK, verstehen wir den Durch- 

schnitt von R mit dem Rechteck 

OS *=%, OX ySy. bo 
Die Menge derjenigen Punkte von 

R, welche auf der Geraden x—0 oder 

auf der Geraden y=O oder auf der 

Kurve C, gelegen sind, wird mit C 0 

bezeichnet. Wir sprechen auch von Abb. 1, Das Grundgebiet R 

der ,,Anfangskurve“ C, weil bei den 

Anfangswertaufgaben fiir hyperbolische Differentialgleichungen auf C Anfangs- 

werte vorgegeben werden. Weiter sei C,,, der Durchschnitt von C und r(x, y). 


Die Klasse der in R stetigen Funktion p(x, y) wird mit C°(f), die Klasse 
der in R stetigen und mit stetigen Ableitungen ¢,, ,, p,, versehenen Funktion 
mit C*(R) bezeichnet. Die iiber einer Menge G summierbaren (d.h. im Lebesgue- 
schen Sinne integrierbaren) Funktionen bilden die Klasse L(G). 

Weiter werden Operatoren K betrachtet, welche den Raum C®(&) in sich 
abbilden; R= (Rk) sei die Klasse dieser Operatoren. Diese Klasse enthalt eine 
fiir alles folgende wichtige Unterklasse %* der monoton wachsenden Operatoren. 
Der Operator!) QER gehért zur Klasse ®*, wenn er die folgende Eigenschaft 
besitzt: 

Sind die Funktionen p(x, y) und @(x, y) aus C°(R) und geniigen ste fiir eonen 
in R gelegenen Punkt (%», Yo) der Unglerchung 


G(% ¥) SPX ¥) mM 7%, Vo) 
so tst 
QpsQGY an der Stelle (xo, Vo). 
Zur Klasse ® gehort z.B. der Operator 


Ko= ff k(x, 9, 5,4, p(s, 8) dsdt, 
r(x, y) 


1 Meist werden Operatoren der Klasse ® mit A’, Operatoren der Klasse &* mit 2 
bezeichnet. 
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wenn die Funktion k(x, y, s, ¢, 2) stetig ist; ist k auBerdem schwach monoton 
wachsend in der letzten Veranderlichen, so ist K€ 8*. Ein anderes Beispiel ftir 
einen monoton wachsenden Operator ist 
v 
Qo= J v(%,t) dt. 


y (%) 


3. Die beiden Hauptsatze 


Es werden nun zwei allgemeine Satze A und B iiber Operatoren bewiesen. 
Sie bilden die Grundlage fiir alle im weiteren Verlauf abzuleitenden Ungleichungen 
und Abschitzungen, welche Volterrasche Integralgleichungen und hyperbolische 
Differentialgleichungen in zwei Veranderlichen betreffen. 


Satz A. Es seien v(x, y) und w(x, y) Funktionen aus C°(R), 4(x, y) eime 
in R erklérte Funktion und 2 € &* ein monoton wachsender Operator. Fiir diese gelte 


(2) vSy +24, w>yt+Qw inR 
(3) v<w aufC. 
Dann ist 

vw mR. 


Hat der Operator Q die Eigenschaft, daf (fiir alle stetigen Funktionen y) Qh=0 
auf C ist, so ist die Voraussetzung (3) entbehrlich (sie folgt dann aus (2)). 

Nehmen wir fiir den Beweis an, die Behauptung sei falsch. Dann gibt es, 
da v<w auf C ist, einen Punkt (%, yp) CR derart, daB v=w in Punkt (%p, Vo) 
und v<w in 7(%, Yo) ist. Da 2 ein monoton wachsender Operator ist, ist dann 
Qv<Qw an der Stelle (%, yo). Fiir diese Stelle (x9, yo) gilt nach (2) 


vSnt+Qvsyn+QNw<w. 


Damit ist ein Widerspruch zur obigen Annahme v(xy, Vp) =w (xg, Vo) erreicht 
worden, der die Richtigkeit von Satz A erhartet. 


Aus dem Beweis geht tibrigens hervor, da die Voraussetzung (2) von Satz A 
auch in 


(27) v<ntQv, w2=n+Q2w inR 
abgeandert werden kann. 


Satz B. Zwischen den Operatoren K ECR, QER*, den Funktionen u, v, 0 C C°(R) 


und den in R erklarten Funktionen (x, y) und d(x, y) mégen die folgenden Bezie- 
hungen 


(4) u=H+Ku, |v—yn—Kv| Sd inR 

(5) |ju—v|<o aufC 

(6) |Kp — K@| <2(|p —G]) fir beliedige g, FE C°(R), 
(7) e>d+Q2o imR 


bestehen. Dann ist 


|u—v|<o mR. 
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Die Voraussetzung (5) ist wieder iiberfliissig, wenn K p und Q@ (fiir beliebige stetige 
Funktion p) auf C verschwinden; sie folgt dann aus (4) und (7). 

Dieser Satz laBt sich sehr einfach auf Satz A zuriickfiihren, indem wir dort v 
durch |w—v|, w durch 9 und 7 durch d ersetzen. Bei dieser Ersetzung wird die 
Voraussetzung (3) mit (5), der zweite Teil der Voraussetzung (2) mit (7) identisch. 
Es bleibt also lediglich noch zu zeigen, daB |v—w|<d+Q(|v—w]) ist. Der 
Nachweis wird mit (4) und (6) gefiihrt: 


|u—u|=|v—n—Ko+Kv—Ku|Sd4+|Kv—Ku|<d+Q(lv—u)). 


Bei den Anwendungen von Satz A liegt meist folgende Situation vor. Es ist 
eine Operatorengleichung u=7+Qu, etwa ein Anfangswertproblem fiir eine 
hyperbolische Differentialgleichung, gegeben, und die Aufgabe besteht darin, obere 
und untere Schranken fiir die Losung wu anzugeben. Man hat dann eine Funktion w 
zu bestimmen, fiir die w>7+2w (und eventuell w> wu auf C, vgl. die Bemerkung 
am SchluB von Satz A) ist. Nun ist es méglich, Satz A auf die Funktionen w 
(statt v) und w anzuwenden; so ergibt sich w< w. Ist man an unteren Schranken 
interessiert, so hat man entsprechend vorzugehen. 


Hat man also zwei Funktionen v, we C(R) bestimmt, fiir die 
(8) v<n+Qu, w>n+tQw mR 


(sowte, falls Q nicht die Eigenschajt ,,.Qp=0 auf C“ besitzt, vcu<w auf C) ist, 
so gilt 


(9) V=—U=wW MER 


fiir jede Losung der Operatorengleichung u=yn+Qu. 

In Ubertragung einer von PERRoN [/] stammenden, bei gewohnlichen und 
elliptischen Differentialgleichungen wohlbekannten Begriffsbildung nennen wir v 
eine Unterfunktion, w eine Oberfunktion. Es sei daran erinnert, daf diesen Be- 
trachtungen die einschneidende Voraussetzung zugrunde liegt, da 2 ein monoton 
wachsender Operator ist. Fiir die Anfangswertprobleme der hyperbolischen 
Differentialgleichung w,,—/(x, y, w) bedeutet das: Der Begriff der Ober- und 
Unterfunktion kann nur fiir solche rechte Seiten f(x, y, z) eingefiihrt werden, die 
in z schwach monoton wachsend sind. 

Gehen wir, um diese Situation zu beleuchten, kurz auf die analogen Ver- 
haltnisse bei gewdhnlichen Differentialgleichungen ein. Wir betrachten das 
Anfangswertproblem y’=f(x, y) fiir OS xa, y(0)=y,. Ist fiir zwei differen- 
zierbare Funktionen v, w 


em al eae) he OS l4ew), tir =< x= 2 


und v(0)< y)<w(0), so ist v< y<w fiir jede Losung y(x) des Problems. Dabei 
sind keinerlei Monotonievoraussetzungen iiber f notwendig. Anders liegen die Dinge, 


wenn man von der entsprechenden Integralgleichung y(*)=yo+/J f(é,y(t)) dé 
0 


ausgeht. Wir betrachten jetzt Lésungen des Anfangswertproblems im Sinne 
von CARATHEODORY, d.h. totalstetige Funktionen y(x), welche dieser Integral- 
gleichung (oder, was dasselbe ist, der Anfangsbedingung y (0) =, und fast tiberall 
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der Gleichung y’=/(x, y)) geniigen. Es gilt dann: Sind die Funktionen v, w 
totalstetig und geniigen sie der Voraussetzung 


v(x) <yo+ S(t») dt, w(x) > vot Sit) de inj (0S *4 
0 0 


oder der (starkeren) Voraussetzung 
v(0)<¥y<w(0), vw Sfx»), 21, OY 1." In OS Fa 


so ist v<y<w, falls die Funktion f(x, 2) im der Variablen z monoton wdachst. 
Anhand einfacher Beispiele laBt sich zeigen, daB diese Monotonievoraussetzung 
notwendig ist. Bei dem in der vorliegenden Arbeit eingeschlagenen Weg zur 
Behandlung hyperbolischer Differentialgleichungen wird ebenfalls von der ent- 
sprechenden Integralgleichung ausgegangen. Eine Monotonievoraussetzung er- 
scheint deshalb unumganglich. 

Der Satz B hingegen ist fiir Operatorengleichungen «= -+ Ku mit beliebigem 
Operator K €§ giiltig. Der monotone Operator 2 tritt hier in der Voraussetzung 
(6) als Schranke fiir den Operator K auf. Meist wird Satz B in folgender Weise 
angewandt. Man versucht, um den Verlauf der unbekannten Lésung w kennenzu- 
lernen, eine ,,Naherungslésung‘‘ v zu bestimmen, die einen mdglichst kleinen 
,Defekt'’ v—1 —Kv besitzt (vgl. (4)), und die sich den auf C vorgegebenen 
Anfangswerten moglichst gut anpaBt (vgl. (5)). Die Schranke g fiir die Differenz 
v—u bestimmt sich aus der GréBe des Defektes und der Abschatzung von K 
durch den monotonen Operator Q. 


4. Anwendung auf Integralgleichungen 
Durch Spezialisierung der beiden allgemeinen Satze der vorigen Nummer auf 
gewisse Integraloperatoren ergeben sich eine Reihe von Resultaten tuber nicht- 
lineare Volterrasche Integralgleichungen in zwei unabhangigen Veranderlichen. 
Die in Rede stehenden Integraloperatoren K und Q2 sind gegeben durch 


(10) Ky [fray sels, t)) ds dt, 
7 (x,y) 

(11) Qy = ff w(x, y, x, 4, y(s, t)) ds dt. 
¥ (x, Vv) 


Die ,,Kerne“ k(x, y, s, f, 2) undw(x, y, s, t, 2) seien auf dem topologischen Produkt 
RXRX{—0o<2< 00} erklart, @ sei auBerdem schwach monoton wachsend in 
der letzten Variablen. Ferner sollen sie solche Regularitatseigenschaften besitzen, 
daB fiir eine beliebige Funktion y€ C°(R) die Integrale (10) und (41) existieren 
und da8 die Funktionen Ky, Qy in R stetig sind und auf C verschwinden?, Diese 
Voraussetzung fassen wir zu REL, wet zusammen; aus ihnen folet KER, 
QER". Es ist nicht schwierig, hinreichende Kriterien dafiir anzugeben, daB 
ke ist. Das ist z.B. der Fall, wenn k(x, y, s, t, z) auf dem ganzen Definitions- 
bereich stetig ist. Eine schwachere Voraussetzung lautet: k sei fiir festes (s, t)C R 
in (x, y) und 2 stetig, fiir festes (x, y)€ R und festes z in (s, 2) meBbar, und zu 


fu, = fast uberall. 


3 Letzteres ist nicht selbstverstandlich; man denke etwa an den durch k (S562) — 
z/y(%, y) definierten Mittelwert-Operator. 
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jedem M> 0 gebe es eine Funktion H(s, t)€ L(R), so daB 
eyes et) | Sore tires ay) CR, |z]}< M 


ist. Offenbar gilt dann ebenfalls kE&. Nun gibt es eine wichtige Klasse von 
Kernen, die sog. Differenzkerne, fiir die auch dieses Kriterium nicht sehr brauch- 
bar ist. Ist k(x, y, s, t, 2) =k, (x —s, y —#, z) ein Differenzkern und wird h, (x, y, z) 
an einer Stelle (%), ¥o) € R unendlich, so versagt das oben angefiihrte Kriterium. 
Bei Volterra-Integralgleichungen in einer Variablen ist die klassische Abelsche 
Integralgleichung von dieser Art; der zugehérige Kern lautet k(x, s, z)=2/]/x—s. 
Das folgende Kriterium fiir Differenzkerne wird wohl in den meisten Fallen aus- 
reichen. Ein Differenzkern k(x, y, s, t, z)=k,(x—s, y—t, z) ist aus &%, wenn 
k, (x, y, 2) stetig in z fiir festes (x, y)€ R und meBbar in (x, y) fiir festes z ist 
und wenn zu jeder positiven Konstanten M eine Funktion H(x, y) € L(R) existiert, 
derart daB 


|k, (x, y,2)| S A(x, y) fiir allez mit |z|<M 
ist. 
Die beiden durch (10) und (11) gegebenen Integraloperatoren AK und 2 haben 
die am Ende von Satz A und Satz B genannte Eigenschaft. Deshalb werden 
diese beiden Satze hier besonders einfach. 


IN, JES: 
v(x, ¥) S(x,y) + Sf ols, 3,60) dsdt 
(12) 1 (x, 9) 
w(x, y)>n(x,y) + Sf w(x, 9, 5,4, w) dsdt 
7 (x, ) 
in R, so gilt 


v<w mR. 


Hierbet ist Voraussetzung, daB v, wEC°(R), n in R erklart und we X* ist. 
B. Es sei w€C°(R) eine Lésung, v€ C°(R) eine Naherungslisung der zu (10) 
gehérigen Integralgleichung; es gelte also in R 
UW (%, V) =a) Sf R(x, y,s,t,u)dsdt 
7 (x, ) 


| v(x, y)— (x, y) Ee, y,S,t,v)dsdt d(x, y). 
r(x, ¥ 


(13) 


Ferner sei REX, wCX*, und es set 


(14) |A(x, y, s,4,2) —R(x, y, s, t,2)|= (x, y, s,t, |z —2)) 


fiir beliebige z,Z und (x, y,s,t)E RXR. Die Funktion @€ C°(R) gentige der Un- 
gleachung 
(15) o(x,y)>d(x,y)+ffo(xy,s,te)dsdt mR. 
v(x, y) 
Dann ist in ganz R 
|u—v|<e. 


Diese beiden Satze folgen, wie gesagt, unmittelbar aus Nr. 4. Sie sind insofern 
von spezieller Art, als wir vorausgesetzt haben, da® die Integralgleichung in Kw 
linear ist. Man kann mit unserer Methode auch allgemeinere Probleme behandeln, 
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etwa Integralgleichungen der Gestalt 
u=h(x,y,u, Ku), 


wobei /# von vier reellen Veranderlichen abhangt. Aufgaben dieser Art werden 
wir indessen nicht weiter verfolgen. 


5. Hyperbolische Differentialgleichungen 
Wir wenden uns zundchst der speziellen hyperbolischen Differentialgleichung 


(16) Uy = T(x, ¥, U) 


zu, bei der auf der rechten Seite die partiellen Ableitungen w,, wu, nicht vor- 
kommen. Von den verschiedenen Problemen fiir die Gleichung (16) lassen sich 
das charakteristische und das nichtcharakteristische oder Cauchysche Anfangs- 
wertproblem in unsere Betrachtungen einbeziehen. Beide Probleme werden in 
der Literatur an vielen Stellen behandelt (etwa KAMKE [3)). 

Beim charakteristischen Anfangswertproblem ist der Grundbereich R das 
Rechteck 0< x<a, OX y<b; es werden auf den Charakteristiken x=0 und y=0 


Anfangswerte 


(17) u(x,0)=o(x) (Sx<a), wO,y)=t(y) (OS ySd) 


vorgeschrieben (o(0)=1(0)). Beim Cauchy-Problem ist das Grundgebiet ein 
krummliniges Dreieck. Die Kurve C, lauft vom Punkt (0, 6) zum Punkt (a, 0). 
Als Anfangswerte werden auf C(=(C,) die Funktionswerte der Lésung sowie 
die Werte der partiellen Ableitungen u, und wu, vorgegeben. Diese diirfen jedoch 
nicht willkiirlich vorgeschrieben werden; vielmehr miissen sie vertraglich sein 
mit der Tatsache, daB fiir jede stetige und mit stetigen ersten Ableitungen ver- 
sehene Funktion (x, y) 


(18) J (peds+9ydt)=I 9.(8,F(s)) ds— J (2), #) dt=o (x9 (9) — 9 FO), 9) 


Cay a 
ist (C,, ist der in r(x, y) gelegene Teil der Kurve C, das erste Integral ist also 
ein tiber dieses Kurvenstiick* erstrecktes Kurvenintegral). 

Was die Differenzierbarkeitsvoraussetzungen anbelangt, so werden der Ein- 
fachheit halber nur Losungen im klassischen Sinne betrachtet. Es wird also 
angenommen, da die Lésung aus der Klasse C*(R) der mit Einschlu8 ihrer 
Ableitungen ,, p, und y,, in R stetigen Funktion @ ist. 

Da die beiden Anfangswertprobleme, wie wir sehen werden, in der Form einer 
nichtlinearen Volterra-Integralgleichung gestellt werden k6énnen, ist es nahe- 
legend, den Begriff der Lésung zu erweitern. Auf eine Méglichkeit der Erwei- 
terung gehen wir kurz ein. Sie schafft einen Lésungsbegriff, der weitgehend 
analog zur Lésung im Sinne von CARATHEODORY bei gewohnlichen Differential- 
gleichungen ist. Sie benutzt den Begriff der Totalstetigkeit fiir Funktionen von 
zwei Verdnderlichen (Hopson [2], S. 346). Fiir totalstetige Funktionen (x, y) 
existiert fast tiberall die Ableitung », y» welche durch 


,(%, y) = lim elath, y+k)+o(*, y)—9(*+h, y)—o(%, y +h) 
ED h,k—>0 ideas oe. ae 


4 Die Kurve C ist ,,von links oben nach rechts unten“ orientiert. 
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definiert ist. Unter einer Lésung etwa der charakteristischen Anfangswertaufgabe 
versteht man dann eine in R stetige und totalstetige Funktion w, welche fast 
tiberall in R der Gleichung (16) sowie der Anfangsbedingung (17) geniigt. Es 
ist nicht schwierig, hierfiir einen allgemeinen Existenzsatz, entsprechend dem 
Existenzsatz von CARATHEODORY fiir gewdhnliche Differentialgleichungen, zu 
beweisen (f(x, y, 2) ist dabei als meBbar in (x, y) und stetig in z vorausgesetzt, 
und es muB eine Abschatzung | f(x, y, z)| S H(x, y) € L(R) bestehen). Eine Theorie 
der Abschatzung, insbesondere der Eindeutigkeit, kann mit den hier entwickelten 
Methoden angegriffen werden; doch sehen wir hier von einer Durchfiihrung ab 
und begniigen uns mit diesem Hinweis. 

Kehren wir zu unserem Problem zuriick. Ist g eine Funktion aus der Klasse 
C*(R), so gilt fiir (x, y)ER 


(x;y) = p(x, 9 (x) +S oy (x \dt= p (x, y (x)) a ‘p, FON+S Pals i)ds| dt. 


Jede Funktion »€C*(R) gestattet also eine eindeutige Zerlegung 


(19) p(x, v) = 9° (x, y) t LD Pays, t) ds dt, 
1 (x,y 
wobei 
(20) P(x, v) =9(% V(x) — fv, dt=o9(F(y), vy) +S yas 


ist (die Ubereinstimmung dieser beiden Ausdriicke fiir y® folgt aus (18); sie ergibt 
sich auch, wenn man bei der obigen Herleitung zuerst in x-Richtung und dann 
erst in y-Richtung integriert). Im Falle eines Rechtecks R hat man einfach 


(21) p(x, Y) = p(x, 0) + (0, v) — ~(0, 0); 


es ist in jedem Falle y°€ C*(R) und g?,=0, also g® von der Gestalt o, (x) + 1, (y). 
Die Formel (19) ist das zweidimensionale Analogon zu der Formel 


y (x) = y+ fy dt, y®°=y(0). 


Jedes der beiden genannten Anfangswertprobleme fiir die Gleichung (16) ist 
also auf eine Integralgleichung 
(22) u(x, y)=n(x, y) +S f(s, 6 w(s, t)) ds dt 

r(x, Vv) 

zuriickfiihrbar, in der 4 (x, vy) =0,(x) +7 (y) € C*(R) ist. Die Funktion 7 ist durch 
die vorgegebenen Anfangswerte eindeutig bestimmt; sind diese beim charakte- 
ristischen Anfangswertproblem durch (17) gegeben, so ist 7=o(x)+7(y) —a(0) 
nach (21). Da die Integralgleichung (22) von der in Nr.4 behandelten Art ist, 
lassen sich die beiden Satze anwenden. 


A. Es sev® fER*, v, wEC°(R) und 


(23) v<ntf{[ft(s,trv)dsdt, w>nt+ff[ftls,t,w)dsdt mR, 
7 (x,y) r(x, ¥) 


5 f€ N+ bedeutet natiirlich, daB k(x, y, s, t, z)=f(s, 4, z)€X* ist. 
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wobei n eine in R erklérte Funktion ist. Dann st in R 
(24) vou<w 
fiir jede Lésung uC C°(R) der Integralgleichung (22). 
Ist insbesondere u,v, wEC*(R), 
(25) OL, se] (%, y; v), Uny =f (%, y; u) , aes y; w) mR, 
(26) YP<w<w mR 
(die Funktion u°, v°, w° reprasentieren die Anfangswerte von u, v, w; sie sind gemap 
(20) definiert), so gilt ebenfalls (24). 

Man beachte, daB die Ungleichungen (26) nicht nur auf der Anfangskurve C, 
sondern im ganzen Bereich R gelten miissen. Liegt etwa das charakteristische 
Anfangswertproblem vor, und sind o(x) bzw. 6(x) die Anfangswerte von u bzw. w 
auf der Charakteristik y=0, t(y) bzw. T(y) die Anfangswerte auf der Charakte- 
ristik «=O, so ist u°=o(x)+7(y) —o(0) und w=a(x)+7(y) —G(0). Es geniigt 
dann fiir das Bestehen der Ungleichung u°< w® nicht, daB o<o und t<T ist. 
Gegenbeispiel : 


—x, T=3—y, w(1,1)=0> w(t, 1) =—F. 


toleo 


C= G=7 =O, .C= 


Der erste Teil der Behauptung A ist ein Spezialfall von Nr.5, der zweite 
Teil laBt sich sofort auf den ersten Teil zuriickfithren (y= 1°). 


B. Die Funktion fC & lasse sich gemaB 
(27) f(x, y.2) —f(%y,2)|\S@(x%y,|z—2Z|) fir (x, y)ER, z, % beliebig 


durch eine Funktion wEC&* abschatzen. Dann besteht, falls die Funktionen 
u,v, @€C°(R) wm Bereich R den Relationen 


(28) u=n+fffls,t.u)dsat, lv—n — [ft(s.t,v) asad < d(x, y) 
(29) o>d+ffa(s,t,o)dsdt 


gentigen, die Abschatzung 


(30) ja—vj<o mR. 
Sind insbesondere die Funktionen u,v, oC C*(R) und ist (mit ener Funktion 
d(x, y) €L(R)) 


(31) yy =1(%, 9,4), — |Yxy — Tf (%, ¥,.0)| SO(x,y) im R, 
(32) Ory 20(%,¥,0)+4(x%,y) mR, 
(33) |uw9—v°]<6 mR, 


so gult wieder die Abschdtzung (30). 


Der erste Teil von B ist wieder in der entsprechenden Aussage von Nr.4 
enthalten. Der zweite Teil ergibt sich so: es ist 


v—u?— Jf f(s, t,v) dsdt| < |W v°| + |v — v°— ff f(s, t, v) ds dt 


S |v — v°| + ff |v,,— f(s, 4,v)|dsdt< |u° — v°| + ff d(s, 2) dsdt. 
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Setzen wir u°=n und |w°—v|+ ff dds dt=d, so gilt (28). Weiterhin ist 


e=0°+ If ox, dsdt>|uo — v| + ff (d+ a(s, 4, 0)) ds dt, 


womit auch (29) nachgewiesen ist. Also gilt in der Tat die Behauptung (30). 


Beziiglich der Voraussetzung (33) ist eine ahnliche Bemerkung wie bei Satz A 
am Platze. Hat beim charakteristischen Anfangswertproblem die Funktion « 
die Anfangswerte a(x), t(y) (auf y=0 bzw. x=0), die Funktion v die Anfangs- 
werte G(x), T(y) und ist |o —G|<e, |r —7|Se, so kann man lediglich behaupten, 
daB | w° —v°|<3¢ ist. Weiter beachte man, daB fiir die Schrankenfunktion @ in 
(29) ein echtes >-Zeichen gefordert wird, wahrend in (32), jedoch nicht in (33), 
das =-Zeichen zugelassen ist. Man kénnte den Satz auch so formulieren, daB 
an allen drei Stellen das =-Zeichen erlaubt ist. Allerdings wiirde man dann 
Maximalintegrale benétigen, ahnlich wie bei gewohnlichen Differentialgleichungen. 
Insbesondere kann man zeigen: 

Ist f(x, y, 2) €X* stetig in Rx{—co<z< co} und so beschaffen, dak héchstens 
eine Losung der Gleichung (22) existiert, so bleibt A richtig, wenn man itiberall die 
Zeichen <, > durch <, = ersetzt. Ebenso kann in B das Gleichheitszeichen in 
(29), (30), (33) hinzugesetzt werden, wenn wCX* die beiden genannten Voraus- 
setzungen der Stetigkeit und Eindeutigkert erfiillt. 


6. Die Lipschitz-Abschatzung 


Von groBem Interesse im Hinblick auf die Anwendungen ist der Fall, daB 
die rechte Seite f(x, y, z) der Differentialgleichung (16) eine Lipschitz-Abschatzung 


(34) |7(%, y, 2) —F(%, y,2)|S k(x, y) |z —2| 


fiir (x, y)€ R und beliebige z, Z zulaBt. Man kann dann, wie aus Nr.5 B folgt, 
Schrankenfunktionen @ als Lésungen der linearen Differentialgleichung 


(35) Oxy =*(x, 9) @ + O(%, ¥) 

bestimmen. Wir betrachten etwas allgemeiner die Integralgleichung 
(36) o(%, y) =a(x, 9) als (s, 4) o(s,“)dsdt 

und setzen dabei d€ C°(R), RE L(R) voraus. Die Lésung lautet 

(37) o(x, y) =4(x, 9) thie, t) d(s, t) w(x, y; 8, t) ds dt, 


wobei y die zur Gleichung (35) gehérende Riemannsche Funktion ist. Sie wird, 
wie man am einfachsten mit Hilfe sukzessiver Approximation beweist, gegeben 
durch 

(38) y(%, 95 8,4) = DR, (%, 9; 8,8) 


n=1 


mit 


xy 
ky=1,) Raley st=JS SRE) 4k, (EN: 8,0) dé dy 
st 
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und hat die Eigenschaften p(x, y; x, t)=y(*, y; 8, y)=1 und 
(39) We1(%, V3 S, 4) =R(s, 4) p(%, ¥5 5,8). 


Durch zweimalige partielle Integration erhalt man aus (37) — allerdings nur 
unter geeigneten Differenzierbarkeitsvoraussetzungen beziiglich d, etwa fiir 
d€C*(R) — die wohlbekannte Formel 


20(x, y) =4(%(y), y) +4(x,9(x)) + 
a) + f (wd, —y.4@) ds — (pd,—y,d) dt] + 2 Jf d.pdsdt 
Cay 


r(x, 9) 
(die beiden ersten Argumente von y lauten tiberall x, y). Fiir das charakteristische 
Anfangswertproblem laBt sich g in die einfachere Gestalt 


o(x, y) =(x, ¥; 0,0) d(0,0) + f p(x, y; s, 0) d,(s, 0) ds + 

(41) y . xy 

+ f(x, y; 0,2) 4,(0, t)dt+f f(x; 5,0) d,,(s,t)dsdt 
0 00 

bringen. 


Es sei uC C*(R) eine Lésung und v€C*(R) eine Naherungslosung der Diffe- 
rentialglerchung (16). Fiir die letztere gelte 


(42) lvry—f(x, ¥,2)|S d(x, y) und |v°—u°|Se(x,y) mR, 


wobet 6(x, y)EL(R), e(x, y)€ CXR) tst. Die Funktion f(x, y, z) gestatte die Ab- 
schatzung (34). Dann wird durch (37) (und damit auch durch (40) bzw. (41), falls 
d€C*(R) ist), wenn man dort fiir d die Funktion 


(43) d(x, y)=e(x, yy) + ff 6(s,t) dsdt 
7 (x,y) 


ernsetzt, ene Schranke o fiir die Differenz v —u geliefert: 
ju—ul|So mR. 


Diese Behauptung beweisen wir, indem wir sie auf Nr.5 B mit w(x, y, z)= 
k(x, y)z zuriickfiihren. Dabei tritt eine kleine Schwierigkeit auf, weil die Funktion 
9 in (29) einer Ungleichung, hier aber der entsprechenden Gleichung geniigt. 
Sie macht folgende zusitzliche Uberlegung notwendig. Bezeichnen wir mit 0» 
die Funktion, die sich gemaB (37) ergibt, wenn dort statt d die Funktion d+1/n 
eingesetzt wird, so haben wir 


g,=d+i4 [[enkasde>a+ ff o,hasdt, 


also nach Nr.5 B, da die Voraussetzung (29) von o, erfiillt wird, u—v|<o,. 
Daraus ergibt sich fiir »—+co die behauptete Abschatzung |w—v|<o. 


In einigen Fallen lat sich die Riemannsche Funktion geschlossen angeben 
(vgl. etwa Saver [7], S. 148ff.). Wir notieren hier nur den Fall, daB 


(44) k(x, v) = h(x) L(y) 
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ist. Die Funktionen /(x) und /(y) seien nichtnegativ und iiber (0, a) bzw. (0, 0) 
summierbar; ihre Stammfunktionen werden mit 


(45) Hy =fh(s)ds, LO) =f) ae 


bezeichnet. Aus (38) leitet man ohne Schwierigkeit ab, daB die Riemannsche 
Funktion durch die Formel 


(46) v(x, 9; 8, t) = E ([A(x) — H(s)| [L(y) — L())) 

gegeben wird, worin EF die modifizierte Bessel-Funktion 0-ter Ordnung 
co {” 

(47) E(t) = Dis = 10(2 Ve) 
n=0 

bedeutet. 


Besonders einfach wird die Abschaétzung, wenn wir k, 6 und « als konstant 
annehmen. Dann ist p(x, y; s, t) =E(k(x—s)(y—d2)), und die Integrale in (40) 
lassen sich in Sonderfallen, von denen wir zwei wichtige angeben, auswerten. 


Charakteristisches Anfangswertproblem. Ist R ein Rechteck, uc C*(R) eine 
Lésung der Differentialgleichung (16), bei der die rechte Seite f einer Lipschitz- 
Bedingung mit der Lipschitz-Konstante k 


(48) [7 (%, 2) — F(%, ¥,2)| Sk z—2| 
geniigt, und gilt fiir die Funktion vE€C*(R) mit geeigneten Konstanten 6 und e 
(49) lu —f(*,y,2)|S6 und |v—wl|<Se mR, 


so besteht im ganzen Rechteck R die Abschatzung 
0) 0) 
(50) ju —v| SB(kxy) (e+ 4)— 4. 


Die Ungleichung |w°—v°|<e gilt nach (21) sicher dann, wenn die Werte 
von uw und v auf den Charakteristiken x=0 und y=O sich um héchstens ¢/3 
unterscheiden. — Im vorliegenden Fall ist d(s, t)=e-+06 st. Die beiden ersten 
Integrale in (41) verschwinden also, und das Gebietsintegral laBt sich auf Grund 
der Gleichungen (46), (47) leicht auswerten. Man erhalt auf diese Weise genau 
die rechte Seite von (50). 

Als zweites Beispiel nehmen wir das Cauchy-Problem, bei dem die Kurve C 
die Gerade x + y=0 ist (die Kurve C verlauft nicht mehr im ersten Quadranten, 
was aber offenbar fiir die ganze Theorie unwesentlich ist). 


Cauchy-Problem mit Anfangskurve x+y=0. Die Funktionen u, v seen 
aus der Klasse® C*(R), und u sei eine Lésung von (16). Fiir geergnete Konstanten 
k, 6, &, &, mégen die Abschitzung (48), die erste der Ungleichungen (49) sowte die 
Ungleichungen’ 

(51) |ju—v|Se aufC und |u,—v,|Se aufC 


6 Unter R verstehen wir jetzt das durch die Ungleichungen ¥+y 20, *Sa, ySb 
definierte Dreieck. 

7 Die zweite Bedingung (51) kann auch durch | u,—v,|<¢, oder allgemeiner durch 
«| 4,—v,|+(1—«) | wy—v,| Se, OSa<1, ersetzt werden. 
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gelten. Dann besteht in R die Abschatzung 


(52) |u —v| <A ltt) 4 Beet”) — : 
mit 
(53) eden sss amet alisha 


Der Beweis kann ahnlich wie oben gefiihrt werden, indem man die Integrale 

n (40) auswertet. Im vorliegenden Fall ist nach (20) |w°—v°|Se+e(*+4), 
also 3 

d(x, y)=e+a(xty) + 


Aus der Formel 


) ener 


NS Sf fos (y—2) "dsdt= er 


—% —t 


sowie ahnlichen Formeln fiir die Kurvenintegrale iiber y, y, und y, ergibt sich 
dann die Behauptung. Eine andere und vielleicht einfachere Art der Herleitung 
besteht darin, fiir die Schranke o den Ansatz 0 (x, y)=g(x-+y) zu machen. Man 
wird dann, mit t=x+ , auf ein Anfangswertproblem fiir eine gewohnliche Diffe- 
rentialgleichung zweiter Ordnung 


p’=kp+d (0StSa+d), g()=e, @PO=4 


gefiihrt, dessen Lésung y(#)=q(x+y) gleich der rechten Seite von (52) ist. 

Aufgaben der eben behandelten Art werden meist in einem um 45° gedrehten 
Koordinatensystem gestellt. Wir schreiben deshalb die Formel (52) auf die 
Koordinaten *+ y=é, x —-y=yn um. Die Differentialgleichung (16) transformiert 
sich dann in 


een a Enki Rs act 
Ll] = thes — Uqy (72 Nae i) 0 


fiir die Funktion u(&,7)=«u(x, y). Die Voraussetzungen lauten jetzt 
L[jv4jJ=0, |LP]|S6 mR 
ja—v| Se und |4%—v,|<Se fir €=0. 
Aus ihnen folgt die Abschatzung 


Ja —d|< Acie Behe 

Die Betrachtungen dieser und der vorangehenden Nummer gelten selbstver- 
standlich auch fiir den Fall, daB u und v zwei Lésungen ein und derselben Glei- 
chung sind. Sie enthalten dann quantitative Aussagen tiber die stetige Abhangig- 
keit der Lésung von den Anfangswerten. La8t man auch die Anfangswerte von 
u und v tibereinstimmen, so hat man zwei Lésungen eines Anfangswertproblems 
und kann aus den obigen Abschatzungen Eindeutigkeitsaussagen ableiten. Auf 
das Eindeutigkeitsproblem sind wir an anderer Stelle [72] ausfithrlich eingegangen, 
weshalb wir uns hier mit der Bemerkung begniigen, daB sich alle bekannten 
Eindeutigkeitskriterien fiir die Gleichung (16) aus dem Abschatzungssatz der 
Nr.5 herleiten lassen. Darin sind enthalten Ubertragungen der Kriterien von 
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LipscHITz, OsGoop, MonTEL, NAGUMO, KAMKE u.a. von gewohnlichen auf hyper- 
bolische Differentialgleichungen. 

Wir wollen zum SchluB noch erwahnen, da8 die in dieser Nummer entwickelten 
Formeln groBe Ahnlichkeit mit den entsprechenden Abschatzungsformeln fiir 
gewohnliche Differentialgleichungen aufweisen. Besonders in die Augen fallend 
ist diese Ubereinstimmung bei der Abschatzung (50). Ihr entspricht bei gewohn- 
lichen Differentialgleichungen die Abschatzung 


x t) 6 
| u(x) —v(x)|<e (e+) Se 
Die Voraussetzungen hierzu lauten u'=f (x, w), |v’ —f (x, v)| <6 fir O40) 
|~(0) —v(0)|Se sowie |f(x, z) —f(x, 2)|<k|z—2|. 


7. Systeme von Integralgleichungen 
Die Ubertragung unserer Uberlegungen auf Systeme von Differential- bzw. 
Integralgleichungen bietet keinerlei Schwierigkeiten. Wir beginnen mit einigen 
mezeichmungen, Essel #24 .. Fiir Vektoren u— (a, ,...:, Ui), U= (Uj 5 »-29-Uy)y 25 
definieren wir |w| durch 
un) 


und Ungleichungen w<v durch die Festsetzung, daB u<v gleichbedeutend mit 
u,<v; fiir i=1,..., ist (entsprechendes gilt fiir <, >, =). Die GréBe |u| ist 
also kein Skalar, sondern ein Vektor. Sind die Komponenten von wu Funktionen, 
so sollen Ausdriicke wie u€C9(R), uC L,(R) besagen, daB alle Komponenten 
u, in der entsprechenden Funktionenklasse liegen. Die Klasse §,, wird von den- 
jenigen Operatoren K gebildet, welche den Raum C?(R) in sich abbilden. Die 
Unterklasse %* der monoton wachsenden Operatoren Q ist genau wie friiher 
charakterisiert durch die Eigenschaft: 


|e] = (| a 


ypeery 


Aus PX in 1(Xo, Vo) folgt A PSA@ an der Stelle (xq, Yo). 
Hierbei ist natiirlich p, PEC), (R) vorausgesetzt. 
Sind z.B. die Komponenten k; des Vektors k Funktionen k; (x, y,s,¢, 2, ..-,2,) = 


an 


k;(x, y, s, t, 2), so wird unter geeigneten Regularitatsvoraussetzungen durch 
(54) Kep = ff k(x, y, 5,1.) ds dt 

1 (%,9) 
ein Operator KES, erklart. Sind alle k; in allen Variablen z; (17, 7) schwach 
monoton wachsend, so ist HE R;. : 

Die beiden Hauptsatze von Nr.3 bleiben ohne irgendwelche Anderungen auch 
fiir Systeme erhalten. Es seien u,v,w,e Funktionen aus der Klasse C}(R), 
y und d in R erklarte Funktionen sowie K ein Operator aus der Klasse &,,, 
Q €&* ein monoton wachsender Operator. 


A, Aus 
(2’) ~syntQv, w>yjntQw nk, 
(3’) v<w aufC 
folgt 


C2.w ink. 
Arch. Rational Mech. Anal., Vol. 7 18a 
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B. Aus 
(4’) u=y+ Ku, |jo—n—Kov|Sd mR, 
(5’) ju—v|<e afc, 
(6’) |Kep —K@|<Q(|\p— |) fir beliebige ~, PcCiR), t 
(7’) e>diQp inR 
folgt 


ju—vl<e inR. 


Gelten auf der Anfangskurve C die Beziehungen® Kp=0, Qp=0 fiir beliebige 
stetige Funktionen ep, so sind die beiden Voraussetzungen (3'), (5") tiberflissig. 

Die Ubereinstimmung mit Nr.3 erstreckt sich auch auf die Beweise. Nur 
am Anfang bei Satz A wird eine kleine Anderung notwendig. Der Punkt (x9, V9) 
ist jetzt so zu wahlen, daB v<w in 7(%, Yo) ist und daB im Punkt (%9, yo) fiir 
mindestens eine Komponente Gleichheit v;(%9, Vp) =w; (9, Yo) herrscht. 

Auf Grund dieser Ubereinstimmung mit dem eindimensionalen Fall bereitet 
es keinerlei Schwierigkeiten, die Betrachtungen der Nr.4 und 5 auf Systeme von 


Integralgleichungen 
u=—n+ {f kG, 1, w) ds dt 
r(x, y) 
(bzw. entsprechende Systeme mit den monoton wachsenden Kernen w (x, y, s, ¢, z)) 
und auf Systeme von hyperbolischen Differentialgleichungen der speziellen Gestalt 


WU, = (%, ¥, UW) 


zu iibertragen. Die einzigen Anderungen an den dort formulierten Ergebnissen 
betreffen die Ersetzung von Skalaren durch Vektoren. Da solche Probleme in 
den Anwendungen selten auftreten, lassen wir es bei dieser Bemerkung bewenden 
und gehen in der nachfolgenden Nummer auf einen anderen weit wichtigeren 
Sonderfall ein. 


8. Die hyperbolische Differentialgleichung Uy =f (x, y, U, Up, Uy) 
Aus den beiden allgemeinen Satzen der vorigen Nummer lassen sich auch 
Aussagen tiber die Lésungen der Differentialgleichung 


(55) Uny =f (%, y,U, Uy, Uy) 


(bei der alle vorkommenden GréBen wieder Skalare sind) ableiten. Beide Anfangs- 
wertprobleme, sowohl das charakteristische wie auch das Cauchysche, lassen sich 
in ein System von drei Integralgleichungen der Form 


U(x, y) =n (x, y) + IS Fs, t, Uy, Ug, Us) ds dt 


r(x,y 


ug 
(56) ty (X,Y) =m (% ¥) + J F(X, b ty, Uy, Ug) at 
¥ (x) 


( 


U3 (Xx, y) = Ny (%, y) ot i AS y, Uy; Us, Us) ds 
x (¥) 


8 0 ist der Nullvektor. 
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iiberfithren (es ist w4=u, uy=u,, Us=u,). Die Funktion 7 ist wie friiher von 
der Gestalt n=o,(x)+7,(y). Anstelle von (56) schreiben wir kurz 


(57) u=y+Fu 


und benutzen dabei die Bezeichnungen w= (1, u2, Us), N=(4, Ny» Ny) Sowie 


y x 
(58) Fp=( ff f(ste)dsat, fi(ete)dt, ff(s,y,.@) ds). 
1 (%, 9) y (*) x (Y) 
Um Schwierigkeiten, welche mit dem eigentlichen Problem wenig zu tun haben, 
von vornherein aus dem Wege zu raumen, setzen wir voraus, daB die Funktion 
T(%, Y, 2%, 4g, %3) auf der Menge Rx{—c< 4%, %, %< 00} stetig ist. Dann ist 
FESR;. Ist die Funktion f in jeder der Variablen z,, 2), 23 monoton wachsend, 
so ist FERS. 

Besitzt die rechte Seite f der Gleichung (55) die zuletzt genannte Monotonie- 
eigenschaft, so gilt fiir F also der Satz A von Nr.7. Wichtig ist nun die folgende 
Bemerkung. Obwohl der Operator F die am Ende von Nr.7 B genannte Eigen- 
schaft, auf C zu verschwinden, nicht besitzt, kann die Voraussetzung (3’) weg- 
gelassen werden, da sie aus (2’) folgt. Es wird also behauptet, daB aus 


(59) oSyntFv, w>yntiFw inR 


sich v<w auf C ergibt; dabei ist v=(v,, vg, v3), W= (W,, Wy, W3) C CZ(R) voraus- 
gesetzt. 

Fiir die ersten Komponenten ist diese Ungleichung v,< w, auf C sicher richtig, 
da die erste Komponente von Fe auf C verschwindet. Aus demselben Grunde 
gilt v»< wz, auf jenem Teil von C, auf dem die zweite Komponente von F@ ver- 
schwindet, das ist also auf der Punktmenge {(x,y(x))| OS <a}. Ebenso ist 
U3<W3 auf {(x (y) y)| OS ys by. Nehmen wir nun an, es sei nicht v»<w, auf 
ganz C. Dann me es auf dem Geradenstiick {( (0, ¥)| Sys b} eine erste Stelle 
Yo> by, an der v2=w, ist (man beachte v(0, b9)<w (0, b9)). Auf dem zwischen 
b, und y, gelegenen Stiick der y-Achse ist dann v<w und deshalb [mit d)=¥ (0) ] 


io Pedra) MOltanldt. 

by bo 
Hieraus und aus (59) folgt aber v.(0, Vo)<w2(0, Yo) im Widerspruch zur oben 
gemachten Annahme 7, (0, Yo) =W2(0, Yo). Es ist also in der Tat v,.<w, auf C 
und, wie man in gleicher Weise zeigt, v3< w, auf C. 

A. Die Funktion } (x, y, 21, 2, 2%) set auf der Menge RX{—0<%, 2%, %< 00} 

stetig und monoton wachsend in den drei letzten Verdnderlichen. Fir die Vektor- 
funktionen u, v, wC C3(R) gelte? 


(60) o<ynt+Fv, u=yn+Fu, w>yntkw mk. 
Dann ist 
(61) vo<u<w mR. 


9 Man beachte, daB es sich im folgenden jeweils um drei skalare Ungleichungen 
handelt. 
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Fiir die Anwendung auf Differentialgleichungen ist folgende Fassung oft bequemer. 
Die skalaren Funktionen u,v, w seien aus der Klasse C*(R), und es ser 


(62) Uygy = f(%, y; v, Vy» Vy) , Uyy = t(%, y, us; Uz, Uy) i Oxy = f(x, ys eo ee Wy). 


Fiir die Anfangswerte dieser Funktionen gelte 
v<u<w auf der Menge {(x,y(x))| OS *S ay} 
(63) V,<U,<w, auf der Menge {(x,¥(x))| OS xSa} 
Vy<Uy<w, auf der Menge {(X(y),¥)| OS ¥S b}. 
Dann ist in R 
(64) V<U<w, Op Shy Vy < Uy < Wy. 
Die Anfangsbedingung (63) lautet im Falle des charakteristischen Anfangs- 
wertproblems 
v(0, 0) << “(0,0)<w(0,0), v,(*,0)< u,(%,0)<w,(*%,0) fir OS%Sa 
vy (0, ¥) <u, (0, y)<w,(0,y) fir OSySd, 
im Falle des Cauchy-Problems 
viu<w, Uy by We, Vy<Uy<w, aufC. 


Der erste Teil von A folgt direkt aus Nr.7 A, da die Voraussetzung (3’), 
wie oben gezeigt wurde, iiberfliissig ist. Der zweite Teil laBt sich mit Hilfe der 
Darstellung (19), (20) leicht auf den ersten zuriickfitihren. Dabei mu8 nachge- 
wiesen werden, da8 im ganzen Bereich R 


0 0 0 0 0 0 0 0 0 
PO WO 1 Saye Mh oO aac lg ity abe 


ist. Es ist nicht schwierig und soll dem Leser iiberlassen werden, diese drei 
Ungleichungen aus der Voraussetzung (63) abzuleiten. 

Die bisherigen, die Differentialgleichung (55) betreffenden Uberlegungen gelten 
nur unter einer einschneidenden Monotonievoraussetzung beziiglich f. Fiir solche 
rechten Seiten # kann man also von Unter- und Oberfunktionen beziiglich eines 
Anfangswertproblems fiir die Gleichung (55) sprechen. Im allgemeinen wird / 
diese Monotonieeigenschaft nicht besitzen, und man muB sich dann auf den 
folgenden Satz B stiitzen. 


B. Die Funktionen f(x, y, 2, 2g, %3) und w(x, Y, 21, 2, 23) Seien in 
Rx{—2<2z,<0o} 


stetig, die Funktion wm auperdem monoton wachsend in %, 2, 23 und es gelte die 
Abschiatzung 


(65) [4(%, 9, 1, 255.25) — f(%, ¥, 2, 2p, %)| S(x,y, 


2 = %|, |Z, — Z|, 


23 — 2s|) 
fiir (x, y)ER und beliebige z;, %;. Fiir die Funktionen u, », 0 €C*(R) gelte 
(66) Uny =](%, Y, U, Uy, My) , | Wry — £(%, y, v, v V,)| S d(x, y) 


eA ose 


(67) Ovy 2 O(%, VY, O, Ox» Oy) + 6(%, y) 
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in R sowte 
|u—v|<e@ auf der Menge {(x,¥(x))| 0S xSa,} 


(68) |u,—v,| <0, auf der Menge {(x,¥(x))| OLE 
|uy—v|<@, auf der Menge {(X(y),y)| OS y Sd}. 


Dann hat man im ganzen Bereich R die Abschatzung 
(69) |u —v|<e@, | uw, — Vy] < @,, |, — vy | < @,. 


Auch dieser Satz kann, genau wie der vorangehende, auf den entsprechenden 
Satz von Nr.7 zuriickgefiihrt werden. Dabei hat man sich zu tiberlegen, daB 
die Voraussetzung (5’) aus (4’) und (7’) folgt, ebenso wie friiher (3’) aus (2’) 
gefolgert wurde. Die Einzelheiten sollen dem Leser iiberlassen bleiben. 

Die in (68) stehenden <-Zeichen sind wohl noch unbequemer als beim ent- 
sprechenden Satz von Nr.5. Stimmen namlich die Ableitungen von w und v auf 
C iiberein, so darf man nicht eine Funktion o wahlen, fiir die 9° konstant ist. 
Doch gilt auch hier die friihere Bemerkung, da8 man durch einen einfachen 
Grenziibergang, der im wesentlichen darauf hinauslauft, fiir die Differential- 
gleichung 0,,=@(, Y, 0, @,, Qy) Maximal- und Minimalintegrale zu definieren, 
zu Satzen gelangen kann, bei denen Gleichheitszeichen zugelassen sind. Im 
besonderen gilt: 

Ist die Funktion f in A bzw. die Funktion w in B so beschaffen, daB ein E xistenz- 
und Eindeutigkertssatz besteht (etwa Lipschitz-Bedingung beztiglich %, %, 23), $0 
kann man ber allen Ungleichungen in A bzw. B das Gleichhettszeichen hinzusetzen, 
ohne die Giiltigkett dieser Satze zu beeintrachtigen. 


9. Die Lipschitz-Bedingung 
Die rechte Seite f der Gleichung (55) geniige einer Lipschitz-Bedingung 


(70) | f (4, Y. 21» 22, 23) — Ff (%, ¥, 21,22, %)| Shlt z,| + pz, zZ,| + g |Z FAR 


bei der k, # und g von x und y abhangende Funktionen sind. Dann 1aBt sich 
eine Schranke 9 fiir die Differenz zweier Funktionen uv und v, welche der Bedin- 
gung (66) geniigen, als Lésung der linearen Differentialgleichung 


(71) Ory =ko+t Por: + 70,19 


gewinnen. 

Gentigt 9 dieser Differentialgleichung sowie der Anfangsbedingung (68), so gilt 
die Abschatzung (69). Durch einen Grenziibergang, ahnlich wie in Nr.6, zergt man, 
da es hinreichend ist, wenn @ die Ungleichungen (68) mit dem <-Zetchen erfiilit. 
In der Behauptung (69) muB dann ebenfalls das <-Zeichen zugelassen werden. 

Da die Lésungen der Gleichung (71) mit Hilfe der Riemannschen Funktion 
geschlossen angegeben werden kénnen, ist es nicht schwierig, hieraus explizite 
Abschitzungen abzuleiten. Wir begniigen uns mit zwei Beispielen, welche auf 
besonders einfache Formeln fiihren. 
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Charakteristisches Anfangswertproblem. Fir geeignete Konstanten k, p, 9, 
6, ¢ mogen die Abschatzungen (70), 


(72) | apy = fey ynent;, %)| Sd PL” son 


lee oo Bes 
|w(0, 0) — (0, 0)| Se 
(73) |“, —v,|Segel fiir y=0 
|u,—v|Seper” fir x=0 
gelten. Dann besteht, wenn u eine Losung der Differentialgleichung (55) ist, die 


Abschatzung 
(74) |\u—v| Se? B (m x y) (e 4 <| = =a mt m=k+pq 


mw Ww 
(sowie entsprechende Abschitzungen fiir die ersten Ableitungen) in R. 
Man rechnet leicht nach, daB die auf der rechten Seite von (74) stehende 
Funktion @ allen an sie gestellten Anforderungen geniigt. 


Cauchy-Problem mit Anfangskurve x+y=0. Fiir geeignete Konstanten 
k, p, 7, 0, € & mogen die Abschdtzungen (70), 


(75) Oe he, y, VU, Vz, Vy)| =O mn fis 


(76) |w—v|Se, |u,—v,| Se, und |u,—v,|Se fir x+y=0 
gelten. Dann gilt in R 
(77) | Na v| oy eh le +y) a A, et) thy 2 


mit 


24=b+at+VGO+tO?+4k, 24,=p+¢—-VO+92+4k 
V@+Pr4RA=a—hlets), VE+O +44 =A(e+5)— a. 


Macht man namlich, wie in Nr.6, fiir @ den Ansatz 0 (x, y)=q@(x+-y), so ergeben 
sich fiir w(t) die Bedingungen 

gy’ =kyo+(p+q)9'+6, gp()=c, y'(0)=4. 
Die Lésung fiihrt gerade auf die Schranke (77). 


10. Ein Beispiel und Schlu8Bbemerkungen 
Wir betrachten die quasilineare Differentialgleichung 


(78) Un =E+thu+ku(u,+4,), 


worin g, h und k stetige Funktionen des Punktes (x, y)€ R sind. Die bekannten 
Existenz- und Eindeutigkeitssatze lassen sich auf diese Gleichung anwenden, 
wenn man sich dabei des folgenden wohlbekannten Kunstgriffes bedient. Sind 
p<gq zwei Zahlen, so bezeichnen wir mit [w]f die Funktion von uw, welche =p 
fiir usp, =u fiir p<u<q und =q fiir w=q ist. Unter Verwendung dieser 
Bezeichnung ersetzt man die Gleichung (78) z.B. durch 
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Man sieht sofort, daB die Funktion f, einer Lipschitz-Bedingung in u, w,, U, 
gentigt und da eine Lésung wu der Differentialgleichung (79) auch Lésung der 
urspriinglichen Gleichung (78) ist, wenn # hinreichend klein und gq hinreichend 
groB ist (genauer: solange <u, u,, u,<q ist). Es gilt also ein Eindeutigkeitssatz 
und ein Existenzsatz fiir die Gleichung (78) (beim letzteren mu8 man sich eventuell 
auf ,,kleine“’ Grundgebiete beschranken, da die Lésungen sehr stark anwachsen 
k6nnen). 

Im folgenden seien g, # und k nichtnegativ, und die durch die Anfangswerte 
eindeutig definierte Funktion (x, y)=0,(*)+1,(y) (vgl. (56)) habe ebenfalls die 
Eigenschaften 720, 7,20, 7,20. Setzen wir in (79) p=0, so ist die Funktion 
v= 0 eine Unterfunktion im Sinne von Nr.8 A. [In (63) sind, da ein Existenzsatz 
besteht, Gleichheitszeichen zugelassen; man kann jedoch auf diese Bemerkung 
verzichten und statt dessen die Unterfunktionen v= —e(1+%+y+ xy), e>0 
betrachten.| Die Gleichung (79) hat demnach genau eine Lésung w, und._ fiir 
diese gilt w, w,,u,20. Dasselbe ist dann auch fiir die urspriingliche Gleichung 
richtig. Unter den am Anfang dieses Abschnittes eingefiihrten Voraussetzungen 
iiber die Koeffizienten g,h, k und die Anfangswerte kann also der Satz A von 
Nr.8 angewandt werden (auch unter Zulassung von Gleichheitszeichen in (63) 
und (64)). 

Bei der Suche nach moglichst guten Ober- und Unterfunktionen kénnen ver- 
schiedene Wege eingeschlagen werden. Man kann etwa einfache und naheliegende 
Funktionsansatze, welche noch freie Parameter enthalten, machen und versuchen, 
die letzteren so zu bestimmen, daB den Bedingungen (62), (63) Geniige getan wird. 
Eine zweite Méglichkeit besteht darin, mit Ansatzen der Form v(x, y)=y(x+y) 
oder v(x, vy) =y(xy) [entsprechend w=@(x-+ y) oder w= (x y)] Ober- und Unter- 
funktionen als Lésungen gewohnlicher Differentialgleichungen zu bestimmen. Ist 
0=2,=2=2., 0S mShNSh,, 0ShSkSk, und t=x+y, so ergeben sich in 
unserem Beispiel fiir g(t), P(é) aus (62) die Ungleichungen 


GP’ SHntmp+2neP"’, P= Sot heP + 2h. 


Dazu kommen dann noch Anfangsbedingungen gema8 (63) (hierbei ist zugelassen, 
daB g,,...,4, Funktionen von t=x-+y sind). In ahnlicher Weise fithrt der 
Ansatz v=y(xy), w= (xy) mit s=xy auf 

yptsy’Sathypt2yskhipy, 


80 ar. —07 a SSeS 
i y+syp’ 28+ hey (a4 =) ape". 


Hierbei wurde vorausgesetzt, daB fiir das Grundgebiet a= gilt; es wurden die 
Ungleichungen 


(81) at SxtyS2yxy 


benutzt. Eine dritte Methode empfiehlt sich besonders dann, wenn ,,gute“ Unter- 
funktionen leicht, Oberfunktionen dagegen schwierig zu bestimmen sind (oder 
umgekehrt). Man setzt dabei, ausgehend von einer Unterfunktion v, w=v-+@ 
und bestimmt o mit Hilfe des Defektes von v nach Satz B von Nr.8. Wir geben 
unten auch dafiir ein Beispiel an. 
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Fiihren wir einige dieser Ansdtze an der Differentialgleichung 
(82) Uy = 1+ U(uU, + Uy) 


durch (h=0, g=k=1)! Es sei fiir das Quadrat OS 4, ySa das charakteristische 
Anfangswertproblem mit Anfangswerten o(x)=1(y)=0 vorgegeben. Der wohl 
einfachste Ansatz v, w=C xy erfiillt die Anfangsbedingungen (63) mit dem Gleich- 
heitszeichen, und aus (62) wird 


C—12,Cx (x+y). 


Fiir C=1 erhilt man eine Unterfunktion, fiir C>1 eine Oberfunktion, das letztere 
jedoch nur, solange das quadratische Polynom P(C)=C*st—C+1S0 ist (es 
war s=xy, t=x-+-y). Nun ist, wie man durch Lésen der entsprechenden qua- 
dratischen Gleichung leicht nachrechnet, P(C))<0 fiir OSstS«, wenn Co= 
(1—/1—4«)/2« ist. Die Funktion w=Cyxy ist also eine Oberfunktion, solange 
xy(x-+ y)<a ist, insbesondere fiir x y(x-+ y)=a. Wir haben also aus dem Ansatz 
C xy zwei Schranken fiir die Lésung 


Vi-—4xy(x+y))=xy+ 2 y(xt+y)+--- 


(83) xySus 


1 
1 
2(4+y) ( 
(x, y=0, 4xy(x-+)<1) gewonnen, welche fiir kleine Werte von x oder y schon 
recht brauchbar ist (0,0625< (4, 3) <0,0646; 0,25<u(%, 3)<0,5). Man itiber- 
legt sich leicht, da8 fiir x y> 0 in (83) an beiden Stellen das <-Zeichen stehen kann. 
Bessere Unter- und Oberfunktionen lassen sich gemaB (80) aus 


yp +sp"S14+2ysypy’, y(0) =0 
yp +se21+(a+Z\9p, (0) =0 


gewinnen. Man erhalt auf diese Weise mit wenig Miihe die Funktionen y=s+ 0,32 3 
und, wenn etwa a= b=} ist, p=s+0,32s?, dh. 


(84) xy +0,32(xy)i<u< xy+0,32(xy)2. 


Offenbar wird beim vorliegenden Beispiel ein Ansatz von der Form w(xy) 
dem Sachverhalt nicht ganz gerecht, da die dann auf der rechten Seite auf- 
tretende GréBe x+y nicht durch xy ausgedriickt werden kann und gemaB (81) 
abgeschatzt werden mu8. Ein Versuch mit dem Ansatz s y(r), wobei r=st= 
xy(x-+y) ist, fiihrt auf 


(85) UASry + ary”! Slt+retr yy! +s%(2x7/— 7) 


fiir y=y/(r). Setzt man y als Polynom 2. Grades an, so ergibt sich mit ys 
1+-7/6+ 77/30 eine Unterfunktion, mit x4=1+7/6+0,0424r? eine Oberfunktion, 
letzteres fiir a=b=}. Damit sind zwei Schranken gefunden, welche sich im 
Quadrat OX x, ySt um weniger als 0,0006 unterscheiden. Beim Versuch, obere 
Schranken fiir gréBere Quadrate aufzustellen, treten Schwierigkeiten auf. Sie 
haben ihren Grund im raschen Anwachsen der Lésung. Versuchen wir, um 


dartiber naheren AufschluB zu bekommen, einen Ansatz x (r)=c/(c—r). Aus (85) 
wird dann 


22(6 — ¢) — er +92 Z)20(c — 1) = (0 fe % <0) 
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Fiir c=6 ergibt sich eine Unterfunktion, fiir c=3 eine Oberfunktion (man beachte 
s?/r<7/4). Die Lésung wird also zwischen y=3 und y=6 unendlich, sie existiert 
sicher in einem Quadrat mit a< 73/2, jedoch nicht im Quadrat mit a= 3, und 
sie geniigt den Ungleichungen 


6xY 3xY 
4 ik ee 
6—x*y(*+4) 3—axy(%¥+y) 


Wir geben zum Schlu8 ein Beispiel fiir die Art und Weise, wie man aus einer 
Unterfunktion mit Satz 8 B eine Oberfunktion ableiten kann. Kennt man bereits 
eine Funktion w(x, y) mit der Eigenschaft ww, u,< Wy, U,Sw,, so darf man 
annehmen, da8 die rechte Seite unserer Differentialgleichung einer Abschatzung 
(65) mit w(x, y, 2, 2_, 23) =% (w,+w,) + w(z,+2,) geniigt (man darf dann, ahnlich 
wie in (79), « durch [w]o, w, durch [u]~, ... ersetzen). Wir gehen, um einen ganz 
einfachen Fall zu betrachten, aus von der Unterfunktion v=xy und setzen 
w—Cxy, C>1. Wir werden sehen, daB man die obere Schranke w nicht explizit 
anzugeben braucht, es mu8 nur sichergestellt sein, daB sie fiir kleine Werte von 
xy oberhalb der Lésung verlauft. Die Gleichung (67) lautet jetzt 


Ory 2 C(% + VOT CXY(O, + Oy) + XY(% + 9). 
Setzt man etwa 90=A x? y?, so hat man 


(86) 4A 23ACKy(x+y)+(%+). 


Es besteht also der folgende Sachverhalt: Solange wSw=Cxy, u,Sw,, usw, 
ist, gilt auch wSv+o=xy+Ax*y*, u,Sv,+0,, UySv,to,, wenn auBerdem 
die Beziehung (86) besteht. Aus dieser Formulierung wird klar, daB man sich 
nicht von vornherein auf eine Funktion w festlegen muB. Der Giiltigkeitsbereich 
der Ungleichungen u<w, u,Sw,, u,Sw, kann vielmehr zugleich mit @ mit- 
bestimmt werden. Wei man ndmlich, daB diese drei Ungleichungen fiir kleine 
xy gelten, so sind sie sicher richtig, solange noch v-+ 9<w, v,+ 0,<W,, ¥y+0,< wy 
ist. Beschranken wir uns etwa auf 0<%, y< 4, so sind die zuletzt genannten 
Ungleichungen fiir 1+-24=C erfiillt. Aus (86) berechnet man dann A zu A= 
4 (13 — 145) <0,32. Die auf diese Weise erhaltene obere Schranke xy+ A x? y2 
ist um ein weniges besser als (84). 

Wir fiigen noch einige allgemeine Bemerkungen an. Das am Ende von Nr.6 
iiber das Eindeutigkeitsproblem Gesagte gilt entsprechend auch fiir die Glei- 
chung (55). Aus dem Satz B von Nr.8 kénnen alle bekannten Eindeutigkeits- 
kriterien fiir die Gleichung (55) hergeleitet werden; vgl. [12]. 

Einige Probleme, auf die wir nicht eingegangen sind, die sich aber mit unserer 
Methode behandeln lassen, seien noch genannt. Das sind zunachst die beiden 
Anfangswertprobleme fiir ein System 


U,, =f (%, ¥, U, U,, Uy). 


Gewisse Differentialgleichungen héherer Ordnung in zwei unabhangigen Variablen, 
z.B. die Gleichung 


Ur vy = (%, Vs Uy Uys Uys Une, Uy) » 
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kénnen ebenfalls auf Integralgleichungen vom Volterra-Typ zuriickgefiihrt wer- 
den. SchlieBlich treten bei der Ausdehnung der ganzen Theorie auf mehr als 
zwei unabhingige Veranderliche keinerlei neue Schwierigkeiten auf. Es ist damit 
méglich, etwa Differentialgleichungen der Gestalt 


Ux y 2 (; y, a, u, Ux» Uy» U,, Ungy >» Ux 2, Uy s) » (a u(x, y; z) 


zu behandeln. 
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